Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/friction_loc.f90
- Timestamp:
- Jul 24, 2024, 2:54:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/friction_loc.f90
r5113 r5116 45 45 CHARACTER(len=80) :: abort_message 46 46 !$OMP THREADPRIVATE(firstcall,friction_type) 47 integer:: jjb,jje47 INTEGER :: jjb,jje 48 48 49 49 !$OMP SINGLE … … 51 51 ! set friction type 52 52 CALL getin("friction_type",friction_type) 53 if ((friction_type<0).or.(friction_type>1)) then53 if ((friction_type<0).or.(friction_type>1)) THEN 54 54 abort_message="wrong friction type" 55 write(lunout,*)'Friction: wrong friction type',friction_type55 WRITE(lunout,*)'Friction: wrong friction type',friction_type 56 56 CALL abort_gcm(modname,abort_message,42) 57 57 endif … … 99 99 ! les deux composantes du vent au pole sont obtenues comme 100 100 ! premiers modes de fourier de v pres du pole 101 if (pole_nord) then 102 101 if (pole_nord) THEN 103 102 upoln=0. 104 103 vpoln=0. … … 119 118 endif 120 119 121 if (pole_sud) then 122 120 if (pole_sud) THEN 123 121 upols=0. 124 122 vpols=0. … … 167 165 endif ! of if (friction_type.eq.0) 168 166 169 if (friction_type==1) then167 if (friction_type==1) THEN 170 168 ! for ucov() 171 169 jjb=jj_begin
Note: See TracChangeset
for help on using the changeset viewer.