Changeset 295 for LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

Timestamp:

Nov 26, 2001, 2:47:57 PM (23 years ago)

Author:

lmdzadmin

Message:

Tag version 0 qui marche en couple/force
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90 (modified) (11 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

@@ -188 +188 @@
   integer, save        :: error
   integer              :: ii, index
-  logical,save              :: check = .true.
+  logical,save              :: check = .false.
   real, dimension(klon):: cal, beta, dif_grnd, capsol
 !!$PB  real, parameter      :: calice=1.0/(5.1444e+06*0.15), tau_gl=86400.*5.
@@ -332 +332 @@
      call albsno(klon,knon,dtime,agesno(:),alb_neig(:), precip_snow(:))  
      where (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0.
-     zfra(1:knon) = max(0.0,min(1.0,snow/(snow+10.0)))
+     zfra(1:knon) = max(0.0,min(1.0,snow(1:knon)/(snow(1:knon)+10.0)))
      alb_new(1 : knon)  = alb_neig(1 : knon) *zfra(1:knon) + &
     &                     alb_new(1 : knon)*(1.0-zfra(1:knon))
@@ -519 +519 @@
       CALL albsno(klon,knon,dtime,agesno(:),alb_neig(:), precip_snow(:))  
       WHERE (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0.
-      zfra(1:knon) = MAX(0.0,MIN(1.0,snow/(snow+10.0)))
+      zfra(1:knon) = MAX(0.0,MIN(1.0,snow(1:knon)/(snow(1:knon)+10.0)))
       alb_new(1 : knon) = alb_neig(1 : knon) *zfra(1:knon) + & 
      &                    0.6 * (1.0-zfra(1:knon))
@@ -589 +589 @@
      CALL albsno(klon,knon,dtime,agesno(:),alb_neig(:), precip_snow(:))  
      WHERE (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0.
-     zfra(1:knon) = MAX(0.0,MIN(1.0,snow/(snow+10.0)))
+     zfra(1:knon) = MAX(0.0,MIN(1.0,snow(1:knon)/(snow(1:knon)+10.0)))
      alb_new(1 : knon)  = alb_neig(1 : knon)*zfra(1:knon) + &
     &                     0.6 * (1.0-zfra(1:knon))
@@ -744 +744 @@
   character (len = 20) :: modname = 'interfsol'
   character (len = 80) :: abort_message
-  logical,save              :: check = .TRUE.
+  logical,save              :: check = .FALSE.
   real, dimension(klon) :: cal, beta, dif_grnd, capsol
 ! type de couplage dans sechiba
@@ -1045 +1045 @@
   cdrag(1:knon) = tq_cdrag(1:knon)
 
-!  where(cdrag > 0.01) 
-!    cdrag = 0.01
-!  endwhere
+   where(cdrag > 0.01) 
+     cdrag = 0.01
+   endwhere
 !  write(*,*)'Cdrag = ',minval(cdrag),maxval(cdrag)
 
@@ -1186 +1186 @@
   character (len = 20) :: modname = 'interfoce_cpl'
   character (len = 80) :: abort_message
-  logical,save              :: check = .true.
+  logical,save              :: check = .FALSE.
 ! variables pour moyenner les variables de couplage
   real, allocatable, dimension(:,:),save :: cpl_sols, cpl_nsol, cpl_rain
@@ -1715 +1715 @@
   character (len = 20),save :: fich ='limit.nc'
   logical, save     :: newlmt = .TRUE.
-  logical, save     :: check = .true.
+  logical, save     :: check = .FALSE.
 ! Champs lus dans le fichier de CL
   real, allocatable , save, dimension(:) :: sst_lu, rug_lu, nat_lu
@@ -1959 +1959 @@
   character (len = 20),save :: fich ='limit.nc'
   logical,save     :: newlmt = .false.
-  logical,save     :: check = .true.
+  logical,save     :: check = .FALSE.
 ! Champs lus dans le fichier de CL
   real, allocatable , save, dimension(:) :: alb_lu, rug_lu
@@ -2138 +2138 @@
   REAL, parameter :: chasno = 3.334E+05/(2.3867E+06*0.15)
 !
-  logical, save         :: check = .true.
+  logical, save         :: check = .FALSE.
   character (len = 20)  :: modname = 'calcul_fluxs'
   logical, save         :: fonte_neige = .false.
@@ -2495 +2495 @@
   REAL, parameter :: chasno = 3.334E+05/(2.3867E+06*0.15)
 !
-  logical, save         :: check = .true.
+  logical, save         :: check = .FALSE.
   character (len = 20)  :: modname = 'fonte_neige'
   logical, save         :: neige_fond = .false.