Changeset 3928 for trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

Timestamp:

Oct 13, 2025, 4:47:04 PM (7 weeks ago)

Author:

jmauxion

Message:

Generic PCM:
Adding a slow_diagfi flag to the run.def/rcm1d.def file for 1D models only. When False, the netcdf
file is opened/closed once, thus saving significant computing time. When true,
the opening frequency is at output frequency (recommended in debug mode). Also
fixing a redundant loop on tracers when writing outputs in physiq_mod.
JM

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -2705 +2705 @@
             ! call write_output(trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
             !              'kg m^-2',qsurf(:,iq) )
-
-            if(watercond.or.CLFvarying)then
-               call write_output("rneb_man","H2O cloud fraction (conv)"," ",rneb_man)
-               call write_output("rneb_lsc","H2O cloud fraction (large scale)"," ",rneb_lsc)
-               call write_output("CLF","H2O cloud fraction"," ",cloudfrac)
-               call write_output("CLFt","H2O column cloud fraction"," ",totcloudfrac)
-               call write_output("RH","relative humidity"," ",RH)
-               call write_output("vteta","virtual potential temperature","K",zh * (1.+(1./epsi-1.) * &
-                       zq(1:ngrid,1:nlayer,1)))
-            endif
-
-            if(waterrain)then
-               call write_output("rain","rainfall","kg m-2 s-1",zdqsrain)
-               call write_output("snow","snowfall","kg m-2 s-1",zdqssnow)
-               call write_output("reevap","reevaporation of precipitation","kg m-2 s-1",reevap_precip)
-            endif
-
-            if(generic_condensation)then
-               call write_output("CLF","GCS cloud fraction"," ",cloudfrac)
-               !AF24: TODO fix rneb_generic and RH_generic failure with write_output()
-               ! call write_output("rneb_generic","GCS cloud fraction (generic condensation)"," ",rneb_generic)
-               ! call write_output("RH_generic","GCS relative humidity"," ",RH_generic)
-               call writediagfi(ngrid,"rneb_generic","GCS cloud fraction (generic condensation)"," ",3,rneb_generic)
-               call writediagfi(ngrid,"RH_generic","GCS relative humidity"," ",3,RH_generic)
-               call write_output("vteta","virtual potential temperature","K",zh * (1.+(1./epsi_generic-1.) * &
-                       zq(1:ngrid,1:nlayer,1) ))
-            endif
-
-            if(generic_rain)then
-               call write_output("rain","generic rainfall","kg m-2 s-1",zdqsrain_generic)
-               call write_output("snow","generic snowfall","kg m-2 s-1",zdqssnow_generic)
-               call write_output("reevap","generic reevaporation of precipitation","kg m-2 s-1",reevap_precip_generic)
-            endif
-
-            if((hydrology).and.(.not.ok_slab_ocean))then
-               call write_output("hice","oceanic ice height","m",hice)
-            endif
-
-            call write_output("tau_col","Total aerosol optical depth","[]",tau_col)
          enddo ! end of 'nq' loop
+
+         if(watercond.or.CLFvarying)then
+            call write_output("rneb_man","H2O cloud fraction (conv)"," ",rneb_man)
+            call write_output("rneb_lsc","H2O cloud fraction (large scale)"," ",rneb_lsc)
+            call write_output("CLF","H2O cloud fraction"," ",cloudfrac)
+            call write_output("CLFt","H2O column cloud fraction"," ",totcloudfrac)
+            call write_output("RH","relative humidity"," ",RH)
+            call write_output("vteta","virtual potential temperature","K",zh * (1.+(1./epsi-1.) * &
+                     zq(1:ngrid,1:nlayer,1)))
+         endif
+
+         if(waterrain)then
+            call write_output("rain","rainfall","kg m-2 s-1",zdqsrain)
+            call write_output("snow","snowfall","kg m-2 s-1",zdqssnow)
+            call write_output("reevap","reevaporation of precipitation","kg m-2 s-1",reevap_precip)
+         endif
+
+         if(generic_condensation)then
+            call write_output("CLF","GCS cloud fraction"," ",cloudfrac)
+            !AF24: TODO fix rneb_generic and RH_generic failure with write_output()
+            ! call write_output("rneb_generic","GCS cloud fraction (generic condensation)"," ",rneb_generic)
+            ! call write_output("RH_generic","GCS relative humidity"," ",RH_generic)
+            call writediagfi(ngrid,"rneb_generic","GCS cloud fraction (generic condensation)"," ",3,rneb_generic)
+            call writediagfi(ngrid,"RH_generic","GCS relative humidity"," ",3,RH_generic)
+            call write_output("vteta","virtual potential temperature","K",zh * (1.+(1./epsi_generic-1.) * &
+                     zq(1:ngrid,1:nlayer,1) ))
+         endif
+
+         if(generic_rain)then
+            call write_output("rain","generic rainfall","kg m-2 s-1",zdqsrain_generic)
+            call write_output("snow","generic snowfall","kg m-2 s-1",zdqssnow_generic)
+            call write_output("reevap","generic reevaporation of precipitation","kg m-2 s-1",reevap_precip_generic)
+         endif
+
+         if((hydrology).and.(.not.ok_slab_ocean))then
+            call write_output("hice","oceanic ice height","m",hice)
+         endif
+
+         call write_output("tau_col","Total aerosol optical depth","[]",tau_col)
        endif ! end of 'tracer'