Changeset 3717 for trunk/LMDZ.GENERIC

Timestamp:

Apr 9, 2025, 11:54:39 AM (3 months ago)

Author:

emillour

Message:

Generic PCM:
Fixing some outputs: when since "write_output" uses the shape of its argument
(scalar, 1D array for surface, 2D for volume) to infer which type of field to
output, one must provide it the sough array (and not just its first element).
MH+SB+JBM+EM

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• libf/phystd/physiq_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/changelog.txt ¶

@@ -2055 +2055 @@
 Follow-up of reindexing of columns for restartfi files: In 1D one should also
 initialize the "ind_cell_glo" array
+
+== 09/04/2025 == MH+SB+JBM+EM
+Fixing some outputs: when since "write_output" uses the shape of its argument
+(scalar, 1D array for surface, 2D for volume) to infer which type of field to
+output, one must provide it the sough array (and not just its first element).

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

@@ -2668 +2668 @@
       ! Output aerosols.
       if (igcm_co2_ice.ne.0.and.iaero_co2.ne.0) &
-         call write_output('CO2ice_reff','CO2ice_reff','m',reffrad(1,1,iaero_co2))
+         call write_output('CO2ice_reff','CO2ice_reff','m',reffrad(:,:,iaero_co2))
       if (igcm_h2o_ice.ne.0.and.iaero_h2o.ne.0) &
          call write_output('H2Oice_reff','H2Oice_reff','m',reffrad(:,:,iaero_h2o))
       if (igcm_co2_ice.ne.0.and.iaero_co2.ne.0) &
-         call write_output('CO2ice_reffcol','CO2ice_reffcol','um kg m^-2',reffcol(1,iaero_co2))
+         call write_output('CO2ice_reffcol','CO2ice_reffcol','um kg m^-2',reffcol(:,iaero_co2))
       if (igcm_h2o_ice.ne.0.and.iaero_h2o.ne.0) &
-         call write_output('H2Oice_reffcol','H2Oice_reffcol','um kg m^-2',reffcol(1,iaero_h2o))
+         call write_output('H2Oice_reffcol','H2Oice_reffcol','um kg m^-2',reffcol(:,iaero_h2o))
 
       ! Output tracers.
@@ -2682 +2682 @@
             call write_output(noms(iq),noms(iq),'kg/kg',zq(:,:,iq))
             call write_output(trim(noms(iq))//'_surf',trim(noms(iq))//'_surf',  &
-                            'kg m^-2',qsurf_hist(1,iq) )
+                            'kg m^-2',qsurf_hist(:,iq) )
             call write_output(trim(noms(iq))//'_col',trim(noms(iq))//'_col',    &
                            'kg m^-2',qcol(:,iq) )