Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3d/friction.F90
- Timestamp:
- Jul 24, 2024, 2:54:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d/friction.F90
r5113 r5116 47 47 ! set friction type 48 48 CALL getin("friction_type", friction_type) 49 if ((friction_type<0).or.(friction_type>1)) then49 if ((friction_type<0).or.(friction_type>1)) THEN 50 50 abort_message = "wrong friction type" 51 write(lunout, *)'Friction: wrong friction type', friction_type51 WRITE(lunout, *)'Friction: wrong friction type', friction_type 52 52 CALL abort_gcm(modname, abort_message, 42) 53 53 endif … … 55 55 ENDIF 56 56 57 if (friction_type==0) then57 if (friction_type==0) THEN 58 58 ! calcul des composantes au carre du vent naturel 59 59 do j = 1, jjp1 … … 118 118 endif ! of if (friction_type.eq.0) 119 119 120 if (friction_type==1) then120 if (friction_type==1) THEN 121 121 do l = 1, llm 122 122 ucov(:, :, l) = ucov(:, :, l) * (1. - pdt * kfrict(l))
Note: See TracChangeset
for help on using the changeset viewer.