- Timestamp:
- Jul 24, 2024, 4:23:34 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/exner_hyb_m.F90
r5116 r5117 45 45 INTEGER ngrid 46 46 REAL p(ngrid,llmp1),pk(ngrid,llm) 47 real, optional:: pkf(ngrid,llm)47 REAL, optional:: pkf(ngrid,llm) 48 48 REAL ps(ngrid),pks(ngrid), alpha(ngrid,llm),beta(ngrid,llm) 49 49 … … 57 57 58 58 ! Sanity check 59 if(firstcall) THEN59 IF (firstcall) THEN 60 60 ! sanity checks for Shallow Water case (1 vertical layer) 61 if(llm==1) THEN62 if(kappa/=1) THEN61 IF (llm==1) THEN 62 IF (kappa/=1) THEN 63 63 CALL abort_gcm(modname, & 64 64 "kappa!=1 , but running in Shallow Water mode!!",42) 65 65 endif 66 if(cpp/=r) THEN66 IF (cpp/=r) THEN 67 67 CALL abort_gcm(modname, & 68 68 "cpp!=r , but running in Shallow Water mode!!",42) 69 69 endif 70 endif ! of if (llm. eq.1)70 endif ! of if (llm.EQ.1) 71 71 72 72 firstcall=.FALSE. … … 74 74 75 75 ! Specific behaviour for Shallow Water (1 vertical layer) case: 76 if(llm==1) THEN76 IF (llm==1) THEN 77 77 ! Compute pks(:),pk(:),pkf(:) 78 78 … … 82 82 ENDDO 83 83 84 if(present(pkf)) THEN84 IF (present(pkf)) THEN 85 85 pkf = pk 86 86 CALL filtreg ( pkf, jmp1, llm, 2, 1, .TRUE., 1 ) … … 89 89 ! our work is done, exit routine 90 90 RETURN 91 endif ! of if (llm. eq.1)91 endif ! of if (llm.EQ.1) 92 92 93 93 ! General case: … … 137 137 ENDDO 138 138 139 if(present(pkf)) THEN139 IF (present(pkf)) THEN 140 140 ! calcul de pkf 141 141 pkf = pk
Note: See TracChangeset
for help on using the changeset viewer.