Changeset 2979 for LMDZ5/trunk

Timestamp:

Aug 3, 2017, 10:25:17 AM (7 years ago)

Author:

Ehouarn Millour

Message:

Missing initializations that caused problems when running in aquaplanet mode.
EM

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phyaqua_mod.F90 (modified) (2 diffs)
• phyetat0.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90

@@ -349 +349 @@
     wake_deltaq = 0.
     wake_s = 0.
+    wake_dens = 0. 
     wake_cstar = 0.
     wake_pe = 0.
@@ -358 +359 @@
     ale_bl_trig =0. 
     alp_bl =0.
-
+    treedrg(:,:,:)=0.
 
 

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -323 +323 @@
   ENDIF
 !FC
-     IF (ifl_pbltree>0) then
+  IF (ifl_pbltree>0) THEN
 !CALL get_field("FTER", pctsrf(:, is_ter), found)
-     treedrg(:,1:klev,1:nbsrf)= 0.0
-   CALL get_field("treedrg_ter", drg_ter(:,:), found)
+    treedrg(:,1:klev,1:nbsrf)= 0.0
+    CALL get_field("treedrg_ter", drg_ter(:,:), found)
 !  found=phyetat0_srf(1,treedrg,"treedrg","drag from vegetation" , 0.)
-  !lecture du profile de freinage des arbres
-   IF (.not. found ) then
-     treedrg(:,1:klev,1:nbsrf)= 0.0
-   else
-     treedrg(:,1:klev,is_ter)= drg_ter(:,:)
+    !lecture du profile de freinage des arbres
+    IF (.not. found ) THEN
+      treedrg(:,1:klev,1:nbsrf)= 0.0
+    ELSE
+      treedrg(:,1:klev,is_ter)= drg_ter(:,:)
 !     found=phyetat0_srf(klev,treedrg,"treedrg","freinage arbres",0.)
-   ENDIF
-      endif
-
+    ENDIF
+  ELSE
+    ! initialize treedrg(), because it will be written in restartphy.nc
+    treedrg(:,:,:) = 0.0 
+  ENDIF
 
   ! Lecture de l'age de la neige:

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -2082 +2082 @@
        CALL zenang_an(iflag_cycle_diurne.GE.1,jH_cur, &
             latitude_deg,longitude_deg,rmu0,fract)
-       JrNt = 1.0
+       swradcorr(:) = 1.0
+       JrNt(:) = 1.0
+       zrmu0(:) = rmu0(:)
     ELSE
        ! recode par Olivier Boucher en sept 2015