Changeset 3811

Timestamp:

Jun 24, 2025, 5:29:48 PM (5 days ago)

Author:

afalco

Message:

Pluto: compute cell_area_for_latlon_outputs even in 1D for outputs to work.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F

@@ -14 +14 @@
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat, volcapa
       use phyredem, only: physdem0,physdem1
-      use geometry_mod, only: init_geometry
+      use geometry_mod, only: init_geometry,
+     &                        init_geometry_cell_area_for_outputs
       use planete_mod, only: apoastr,periastr,year_day,peri_day,
      &         obliquit,z0,lmixmin,emin_turb,coefvis,coefir,
@@ -31 +32 @@
       use planete_mod, only: ini_planete_mod
       use physics_distribution_mod, only: init_physics_distribution
-      use regular_lonlat_mod, only: init_regular_lonlat
       use mod_interface_dyn_phys, only: init_interface_dyn_phys
       use inifis_mod, only: inifis
@@ -557 +557 @@
      &                   (/0.,0.,0.,0./),(/0.,0.,0.,0./),
      &                   cell_area,[1])
+      call init_geometry_cell_area_for_outputs(1,cell_area)
 ! Ehouarn: init_vertial_layers called later (because disvert not called yet)
 !      call init_vertical_layers(nlayer,preff,scaleheight,