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

Timestamp:

Sep 26, 2023, 3:57:38 PM (14 months ago)

Author:

jbclement

Message:

Mars PEM:
Minor changes concerning the form of the code in the PEM.
JBC

File:

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

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

@@ -80 +80 @@
       B=1/m_noco2
 
-  print *, "Opening ", fichnom, "..."
+  write(*,*) "Opening ", fichnom, "..."
 
 !  Open initial state NetCDF file
@@ -86 +86 @@
   CALL err(NF90_OPEN(var,NF90_NOWRITE,fID),"open",var)
 
-     print *, "Downloading data for vmr co2..."
+     write(*,*) "Downloading data for vmr co2..."
 
   CALL get_var3("co2_layer1"   ,q_co2_dyn)
 
-     print *, "Downloading data for vmr co2 done"
-     print *, "Downloading data for vmr h20..."
+     write(*,*) "Downloading data for vmr co2 done"
+     write(*,*) "Downloading data for vmr h20..."
 
   CALL get_var3("h2o_layer1"   ,q_h2o_dyn)
 
-     print *, "Downloading data for vmr h2o done"
-     print *, "Downloading data for surface pressure ..."
+     write(*,*) "Downloading data for vmr h2o done"
+     write(*,*) "Downloading data for surface pressure ..."
 
   CALL get_var3("ps"   ,ps_GCM)
 
-     print *, "Downloading data for surface pressure done"
-     print *, "nslope=", nslope
+     write(*,*) "Downloading data for surface pressure done"
+     write(*,*) "nslope=", nslope
 
 if(nslope.gt.1) then
 
-     print *, "Downloading data for co2ice_slope ..."
+     write(*,*) "Downloading data for co2ice_slope ..."
 
 DO islope=1,nslope
@@ -112 +112 @@
 ENDDO
 
-     print *, "Downloading data for co2ice_slope done"
-     print *, "Downloading data for h2o_ice_s_slope ..."
+     write(*,*) "Downloading data for co2ice_slope done"
+     write(*,*) "Downloading data for h2o_ice_s_slope ..."
 
 DO islope=1,nslope
@@ -120 +120 @@
 ENDDO
 
-     print *, "Downloading data for h2o_ice_s_slope done"
+     write(*,*) "Downloading data for h2o_ice_s_slope done"
 
 #ifndef CPP_STD
 
-     print *, "Downloading data for watercap_slope ..."
+     write(*,*) "Downloading data for watercap_slope ..."
 DO islope=1,nslope
        write(num,fmt='(i2.2)') islope
@@ -130 +130 @@
 !        watercap_slope(:,:,:,:)= 0.
 ENDDO            
-     print *, "Downloading data for watercap_slope done"
+     write(*,*) "Downloading data for watercap_slope done"
 
 #endif
     
- print *, "Downloading data for tsurf_slope ..."
+ write(*,*) "Downloading data for tsurf_slope ..."
 
 DO islope=1,nslope
@@ -141 +141 @@
 ENDDO
 
-     print *, "Downloading data for tsurf_slope done"
+     write(*,*) "Downloading data for tsurf_slope done"
 
 #ifndef CPP_STD
@@ -147 +147 @@
      if(soil_pem) then
 
-     print *, "Downloading data for soiltemp_slope ..."
+     write(*,*) "Downloading data for soiltemp_slope ..."
 
 DO islope=1,nslope
@@ -154 +154 @@
 ENDDO
 
-     print *, "Downloading data for soiltemp_slope done"
-
-     print *, "Downloading data for watersoil_density ..."
+     write(*,*) "Downloading data for soiltemp_slope done"
+
+     write(*,*) "Downloading data for watersoil_density ..."
 
 DO islope=1,nslope
@@ -163 +163 @@
 ENDDO
 
-     print *, "Downloading data for  watersoil_density  done"
-
-     print *, "Downloading data for  watersurf_density  ..."
+     write(*,*) "Downloading data for  watersoil_density  done"
+
+     write(*,*) "Downloading data for  watersurf_density  ..."
 
 DO islope=1,nslope
@@ -172 +172 @@
 ENDDO
 
-     print *, "Downloading data for  watersurf_density  done"
+     write(*,*) "Downloading data for  watersurf_density  done"
 
   endif !soil_pem
@@ -196 +196 @@
 
 ! Compute the minimum over the year for each point
-  print *, "Computing the min of h2o_ice_slope"
+  write(*,*) "Computing the min of h2o_ice_slope"
   min_h2o_ice_dyn(:,:,:)=minval(h2o_ice_s_dyn+watercap_slope,4)
-  print *, "Computing the min of co2_ice_slope"
+  write(*,*) "Computing the min of co2_ice_slope"
   min_co2_ice_dyn(:,:,:)=minval(co2_ice_slope_dyn,4)
 
 !Compute averages
 
-    print *, "Computing average of tsurf"
+    write(*,*) "Computing average of tsurf"
     tsurf_ave_dyn(:,:,:)=SUM(tsurf_gcm_dyn(:,:,:,:),4)/timelen
 
@@ -215 +215 @@
 
   if(soil_pem) then
-    print *, "Computing average of tsoil"
+    write(*,*) "Computing average of tsoil"
     tsoil_ave_dyn(:,:,:,:)=SUM(tsoil_gcm_dyn(:,:,:,:,:),5)/timelen
-    print *, "Computing average of watersurf_density"
+    write(*,*) "Computing average of watersurf_density"
     watersurf_density_ave(:,:) = SUM(watersurf_density(:,:,:),3)/timelen
   endif