- Timestamp:
- Jul 24, 2024, 2:54:37 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/exner_hyb_loc_m.F90
r5106 r5116 58 58 logical, save :: firstcall = .TRUE. 59 59 !$OMP THREADPRIVATE(firstcall) 60 character(len= *), parameter :: modname = "exner_hyb_loc"60 CHARACTER(LEN = *), parameter :: modname = "exner_hyb_loc" 61 61 62 62 !$OMP BARRIER 63 63 64 64 ! Sanity check 65 if (firstcall) then65 if (firstcall) THEN 66 66 ! sanity checks for Shallow Water case (1 vertical layer) 67 if (llm==1) then68 if (kappa/=1) then67 if (llm==1) THEN 68 if (kappa/=1) THEN 69 69 CALL abort_gcm(modname, & 70 70 "kappa!=1 , but running in Shallow Water mode!!", 42) 71 71 endif 72 if (cpp/=r) then72 if (cpp/=r) THEN 73 73 CALL abort_gcm(modname, & 74 74 "cpp!=r , but running in Shallow Water mode!!", 42) … … 82 82 83 83 ! Specific behaviour for Shallow Water (1 vertical layer) case: 84 if (llm==1) then 85 84 if (llm==1) THEN 86 85 ! Compute pks(:),pk(:),pkf(:) 87 86 ijb = ij_begin … … 96 95 97 96 !$OMP BARRIER 98 if (present(pkf)) then97 if (present(pkf)) THEN 99 98 jjb = jj_begin 100 99 jje = jj_end … … 104 103 105 104 ! our work is done, exit routine 106 return105 RETURN 107 106 endif ! of if (llm.eq.1) 108 107 … … 170 169 ENDDO 171 170 172 if (present(pkf)) then171 if (present(pkf)) THEN 173 172 ! calcul de pkf 174 173
Note: See TracChangeset
for help on using the changeset viewer.