- Timestamp:
- Jul 24, 2024, 2:54:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/exner_milieu_m.F90
r5106 r5116 51 51 52 52 logical,save :: firstcall=.TRUE. 53 character(len=*),parameter :: modname="exner_milieu"53 CHARACTER(LEN=*),parameter :: modname="exner_milieu" 54 54 55 55 ! Sanity check 56 if (firstcall) then56 if (firstcall) THEN 57 57 ! sanity checks for Shallow Water case (1 vertical layer) 58 if (llm==1) then59 if (kappa/=1) then58 if (llm==1) THEN 59 if (kappa/=1) THEN 60 60 CALL abort_gcm(modname, & 61 61 "kappa!=1 , but running in Shallow Water mode!!",42) 62 62 endif 63 if (cpp/=r) then63 if (cpp/=r) THEN 64 64 CALL abort_gcm(modname, & 65 65 "cpp!=r , but running in Shallow Water mode!!",42) … … 71 71 72 72 ! Specific behaviour for Shallow Water (1 vertical layer) case: 73 if (llm==1) then 74 73 if (llm==1) THEN 75 74 ! Compute pks(:),pk(:),pkf(:) 76 75 … … 80 79 ENDDO 81 80 82 if (present(pkf)) then81 if (present(pkf)) THEN 83 82 pkf = pk 84 83 CALL filtreg ( pkf, jmp1, llm, 2, 1, .TRUE., 1 ) … … 86 85 87 86 ! our work is done, exit routine 88 return87 RETURN 89 88 endif ! of if (llm.eq.1) 90 89 … … 117 116 ENDDO 118 117 119 if (present(pkf)) then118 if (present(pkf)) THEN 120 119 ! calcul de pkf 121 120 pkf = pk
Note: See TracChangeset
for help on using the changeset viewer.