Changeset 2979 for LMDZ5/trunk/libf/phylmd/phyetat0.F90

Timestamp:

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

Author:

Ehouarn Millour

Message:

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

File:

• LMDZ5/trunk/libf/phylmd/phyetat0.F90 (modified) (1 diff)

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: