Changeset 1064

Timestamp:

Dec 17, 2008, 1:27:20 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Correction: initialisation rh2m, qsat2m en dehors boucle nbsrf
IM

File:

• LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -1 @@
-!
-! $Header$
 !
 MODULE pbl_surface_mod
@@ -152 +150 @@
 
     IF (type_ocean /= 'slab  ' .AND. type_ocean /= 'force ' .AND. type_ocean /= 'couple') THEN
-      WRITE(lunout,*)' *** Warning ***'
-      WRITE(lunout,*)'Option couplage pour l''ocean = ', type_ocean
-      abort_message='option pour l''ocean non valable'
-      CALL abort_gcm(modname,abort_message,1)
+       WRITE(lunout,*)' *** Warning ***'
+       WRITE(lunout,*)'Option couplage pour l''ocean = ', type_ocean
+       abort_message='option pour l''ocean non valable'
+       CALL abort_gcm(modname,abort_message,1)
     ENDIF
 
@@ -255 +253 @@
     INCLUDE "temps.h"
     INCLUDE "control.h"
-
 ! Input variables
 !****************************************************************************************
@@ -545 +542 @@
     ytsoil = 999999. 
 
+    rh2m(:)        = 0.
+    qsat2m(:)      = 0.
 !****************************************************************************************
 ! 3) - Calculate pressure thickness of each layer
@@ -642 +641 @@
 ! 
 !****************************************************************************************
-    
+   
     loop_nbsrf: DO nsrf = 1, nbsrf
 
@@ -988 +987 @@
        trmb2(:,nsrf)  = 0.        ! inhibition 
        trmb3(:,nsrf)  = 0.        ! Point Omega
-       rh2m(:)        = 0.
-       qsat2m(:)      = 0.
 
 #undef T2m