Changeset 3717 for trunk/LMDZ.GENERIC/libf/phystd

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

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (2 diffs)

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) )