Changeset 2885 for trunk/LMDZ.COMMON/libf/evolution/read_data_GCM.F90

Timestamp:

Jan 31, 2023, 9:52:30 PM (22 months ago)

Author:

romain.vande

Message:

Mars PCM:
Move a endif misplaced
RV

File:

• trunk/LMDZ.COMMON/libf/evolution/read_data_GCM.F90 (modified) (5 diffs)

trunk/LMDZ.COMMON/libf/evolution/read_data_GCM.F90

@@ -69 +69 @@
   REAL, ALLOCATABLE ::  co2_ice_s(:,:,:)                               ! co2 ice, mesh averaged, of the concatenated file [kg/m^2]
   REAL, ALLOCATABLE ::  h2o_ice_s_slope(:,:,:,:)                       ! h2o ice per slope of the concatenated file [kg/m^2]
+  REAL, ALLOCATABLE :: watercap_slope(:,:,:,:)
 
 !-----------------------------------------------------------------------
@@ -80 +81 @@
       allocate(co2_ice_s(iim+1,jjm+1,timelen))
       allocate(h2o_ice_s_slope(iim+1,jjm+1,nslope,timelen))
+      allocate(watercap_slope(iim_input+1,jjm_input+1,nslope,timelen))
       allocate(h2o_ice_s(iim+1,jjm+1,timelen))
 
@@ -133 +135 @@
 
      print *, "Downloading data for h2o_ice_s_slope done"
+     print *, "Downloading data for watercap_slope ..."
+
+DO islope=1,nslope
+  write(num,fmt='(i2.2)') islope
+  call get_var3("watercap_slope"//num,watercap_slope(:,:,islope,:))
+ENDDO
+
+     print *, "Downloading data for watercap_slope done"
      print *, "Downloading data for tsurf_slope ..."
 
@@ -198 +208 @@
 
   print *, "Computing the min of h2o_ice_slope"
-  min_h2o_ice_slope(:,:,:)=minval(h2o_ice_s_slope,4)
+  min_h2o_ice_slope(:,:,:)=minval(h2o_ice_s_slope+watercap_slope,4)
   print *, "Computing the min of co2_ice_slope"
   min_co2_ice_slope(:,:,:)=minval(co2_ice_slope,4)
@@ -225 +235 @@
             min_co2_ice_slope(i,j,islope)  = 0.
           endif
-          if (min_h2o_ice_slope(i,j,islope).LT.0) then
-            min_h2o_ice_slope(i,j,islope)  = 0.
-          endif
+!          if (min_h2o_ice_slope(i,j,islope).LT.0) then
+!            min_h2o_ice_slope(i,j,islope)  = 0.
+!          endif
        ENDDO
     ENDDO