NLM_LBC in ROMS 3.7

[salik]

I've been running a modified version of the upwelling case with (among other modifications), a uniform 1 m thick sediment bed in ROMS 3.7
When trying to change the sediment bed to have a non-uniform thickness in the entire domain, I added the following code to ana_sediment:

Code: Select all

> #  if defined unised
>       DO j=JstrT,JendT
> 	DO i=IstrT,IendT
> 	  DO k=1,Nbed
>             bed(i,j,k,iaged)=time(ng)
>             bed(i,j,k,ithck)=0.10_r8
>             bed(i,j,k,iporo)=0.90_r8 
>             bed_frac(i,j,k,1)=1.0_r8
>           END DO
> !  Set exposed sediment layer properties.
> 
>           bottom(i,j,irlen)=0.10_r8
>           bottom(i,j,irhgt)=0.01_r8
>           bottom(i,j,izdef)=Zob(ng)
> #  elif defined SHELFSED
>       DO j=JstrT,JendT
> 	DO i=IstrT,IendT
> 	  val4=REAL(Lm(ng)+1-i,r8)
> !  Set bed layer properties.
> 	  IF(val4.le.20.0_r8) THEN
>             DO k=1,Nbed
>               bed(i,j,k,iaged)=time(ng)
>               bed(i,j,k,ithck)=0.10_r8
>               bed(i,j,k,iporo)=0.90_r8 
>               bed_frac(i,j,k,1)=1.0_r8
>             END DO
> !  Set exposed sediment layer properties.
> 
>             bottom(i,j,irlen)=0.10_r8
>             bottom(i,j,irhgt)=0.01_r8
>             bottom(i,j,izdef)=Zob(ng)
>           ELSE IF(val4.le.25.0_r8) THEN  
>             DO k=1,Nbed
>               bed(i,j,k,iaged)=time(ng)
>               bed(i,j,k,ithck)=0.10_r8*cos((val4-21)*pi/8)
>               bed(i,j,k,iporo)=0.90_r8
>               bed_frac(i,j,k,1)=1.0_r8
>             END DO
> !  Set exposed sediment layer properties.
> 
>             bottom(i,j,irlen)=0.10_r8
>             bottom(i,j,irhgt)=0.01_r8
>             bottom(i,j,izdef)=Zob(ng)
>           ELSE
>             DO k=1,Nbed
>               bed(i,j,k,iaged)=time(ng)
>               bed(i,j,k,ithck)=0.0_r8
>               bed(i,j,k,iporo)=0.9_r8
>               bed_frac(i,j,k,1)=1.0_r8
>             END DO
> !  Set exposed sediment layer properties.
> 
>             bottom(i,j,irlen)=0.0_r8
>             bottom(i,j,irhgt)=0.0_r8
>             bottom(i,j,izdef)=Zob(ng)
>           END IF
>         END DO
>       END DO
> #  endif

My model runs in three phases with different forcing, so I'm restarting from the RST-file a lot. This has never been a problem until now.
Since making this change, I've been getting the LBC_GETATT error when trying to restart from a previous solution.
I have tried changing the ana_sediment.h file back, I've tried overwriting it with a different file that I know to work, but this does not solve the problem.

Does anyone have experience with this issue? I guess I could just activate the #define NO_LBC_ATT option, but it seems like this should not be necessary when running 3.7