Changeset 3150 for trunk/LMDZ.MARS

Timestamp:

Dec 6, 2023, 8:05:33 PM (12 months ago)

Author:

llange

Message:

Mars PCM
Switching vdif_cd to a module and convert it to F90.
Some cleaning have also been made (remove commented lines, outputs, etc.)
VDIFC: minor fix: zcdv_true zcdh_true are out and not in variables.
LL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/vdif_cd_mod.F90 (added)
• libf/phymars/vdifc_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4381 +4381 @@
 == 06/12/2023 == JBC
 Tiny fix from r3147 (Fortran fixed format).
- 
+
+== 06/12/2023 == LL
+Switching vdif_cd to a module and convert it to F90.
+Some cleaning have also been made (remove commented lines, outputs, etc.)
+VDIFC: minor fix: zcdv_true zcdh_true are out and not in variables.

trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F

@@ -35 +35 @@
       use paleoclimate_mod, only: d_coef,h2o_ice_depth,lag_layer
       use comsoil_h, only: layer, mlayer,adsorption_soil
-
+      use vdif_cd_mod, only: vdif_cd
       IMPLICIT NONE
 
@@ -113 +113 @@
       REAL zkq(ngrid,nlay+1)
       REAL zcdv(ngrid),zcdh(ngrid)
-      REAL, INTENT(IN) :: zcdv_true(ngrid),zcdh_true(ngrid)    ! drag coeff are used by the LES to recompute u* and hfx
+      REAL, INTENT(OUT) :: zcdv_true(ngrid),zcdh_true(ngrid)    ! drag coeff are used by the LES to recompute u* and hfx
       REAL zu(ngrid,nlay),zv(ngrid,nlay)
       REAL zh(ngrid,nlay)
@@ -1410 +1410 @@
      &                     "Ground ice latent heat flux",
      &                     "W.m-2",surf_h2o_lh(:,iflat))
-         call write_output('zdqsdif_ssi_frost',
-     &          'Flux between frost and subsurface','kg.m-2.s-1',
-     &                zdqsdif_ssi_frost(:,1))
-
-         call write_output('zdq_subtimestep',
-     &          'Actual flux zdqsdif_surf*subtimestep',
-     &          'kg.m-2',zdqsdif_tot(:)*subtimestep)
+!         call write_output('zdqsdif_ssi_frost',
+!     &          'Flux between frost and subsurface','kg.m-2.s-1',
+!     &                zdqsdif_ssi_frost(:,1))
+
+!         call write_output('zdq_subtimestep',
+!     &          'Actual flux zdqsdif_surf*subtimestep',
+!     &          'kg.m-2',zdqsdif_tot(:)*subtimestep)
 !         call write_output('zdq_end',
 !     &          'Flux after all contributions',