52 |
_RL sm_xz (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
_RL sm_xz (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
53 |
_RL sm_yz (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
_RL sm_yz (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
54 |
|
|
|
#ifdef GAD_ALLOW_SOM_ADVECT |
|
55 |
C !LOCAL VARIABLES: ==================================================== |
C !LOCAL VARIABLES: ==================================================== |
56 |
C i,j,k :: loop indices |
C i,j,k :: loop indices |
57 |
_RL three |
_RL three |
74 |
& MAX(ABS(s1new)-slpmax,sm_zz(i,j,k)) ) |
& MAX(ABS(s1new)-slpmax,sm_zz(i,j,k)) ) |
75 |
sm_xz(i,j,k) = MIN( slpmax, MAX(-slpmax,sm_xz(i,j,k)) ) |
sm_xz(i,j,k) = MIN( slpmax, MAX(-slpmax,sm_xz(i,j,k)) ) |
76 |
sm_yz(i,j,k) = MIN( slpmax, MAX(-slpmax,sm_yz(i,j,k)) ) |
sm_yz(i,j,k) = MIN( slpmax, MAX(-slpmax,sm_yz(i,j,k)) ) |
77 |
sm_z (i,j,k) = s1new ; |
sm_z (i,j,k) = s1new |
78 |
sm_zz(i,j,k) = s2new ; |
sm_zz(i,j,k) = s2new |
79 |
ENDDO |
ENDDO |
80 |
ENDDO |
ENDDO |
81 |
ENDDO |
ENDDO |
82 |
ENDIF |
ENDIF |
83 |
|
|
|
#endif /* GAD_ALLOW_SOM_ADVECT */ |
|
|
|
|
84 |
RETURN |
RETURN |
85 |
END |
END |