- Timestamp:
- Jul 24, 2024, 4:23:34 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/exner_milieu_loc_m.F90
r5116 r5117 56 56 57 57 ! Sanity check 58 if(firstcall) THEN58 IF (firstcall) THEN 59 59 ! sanity checks for Shallow Water case (1 vertical layer) 60 if(llm==1) THEN61 if(kappa/=1) THEN60 IF (llm==1) THEN 61 IF (kappa/=1) THEN 62 62 CALL abort_gcm(modname, & 63 63 "kappa!=1 , but running in Shallow Water mode!!",42) 64 64 endif 65 if(cpp/=r) THEN65 IF (cpp/=r) THEN 66 66 CALL abort_gcm(modname, & 67 67 "cpp!=r , but running in Shallow Water mode!!",42) 68 68 endif 69 endif ! of if (llm. eq.1)69 endif ! of if (llm.EQ.1) 70 70 71 71 firstcall=.FALSE. … … 75 75 76 76 ! Specific behaviour for Shallow Water (1 vertical layer) case: 77 if(llm==1) THEN77 IF (llm==1) THEN 78 78 ! Compute pks(:),pk(:),pkf(:) 79 79 ijb=ij_begin … … 83 83 pks(ij) = (cpp/preff) * ps(ij) 84 84 pk(ij,1) = .5*pks(ij) 85 if(present(pkf)) pkf(ij,1)=pk(ij,1)85 IF (present(pkf)) pkf(ij,1)=pk(ij,1) 86 86 ENDDO 87 87 !$OMP ENDDO 88 88 89 89 !$OMP BARRIER 90 if(present(pkf)) THEN90 IF (present(pkf)) THEN 91 91 jjb=jj_begin 92 92 jje=jj_end … … 97 97 ! our work is done, exit routine 98 98 RETURN 99 endif ! of if (llm. eq.1)99 endif ! of if (llm.EQ.1) 100 100 101 101 ! General case: … … 140 140 !$OMP ENDDO NOWAIT 141 141 142 if(present(pkf)) THEN142 IF (present(pkf)) THEN 143 143 ! calcul de pkf 144 144
Note: See TracChangeset
for help on using the changeset viewer.