- 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_milieu_m.F90
r5116 r5117 42 42 INTEGER ngrid 43 43 REAL p(ngrid,llmp1),pk(ngrid,llm) 44 real, optional:: pkf(ngrid,llm)44 REAL, optional:: pkf(ngrid,llm) 45 45 REAL ps(ngrid),pks(ngrid) 46 46 … … 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) 66 66 endif 67 endif ! of if (llm. eq.1)67 endif ! of if (llm.EQ.1) 68 68 69 69 firstcall=.FALSE. … … 71 71 72 72 ! Specific behaviour for Shallow Water (1 vertical layer) case: 73 if(llm==1) THEN73 IF (llm==1) THEN 74 74 ! Compute pks(:),pk(:),pkf(:) 75 75 … … 79 79 ENDDO 80 80 81 if(present(pkf)) THEN81 IF (present(pkf)) THEN 82 82 pkf = pk 83 83 CALL filtreg ( pkf, jmp1, llm, 2, 1, .TRUE., 1 ) … … 86 86 ! our work is done, exit routine 87 87 RETURN 88 endif ! of if (llm. eq.1)88 endif ! of if (llm.EQ.1) 89 89 90 90 ! General case: … … 116 116 ENDDO 117 117 118 if(present(pkf)) THEN118 IF (present(pkf)) THEN 119 119 ! calcul de pkf 120 120 pkf = pk
Note: See TracChangeset
for help on using the changeset viewer.