Changeset 4608 for LMDZ6/trunk/libf/phylmd

Timestamp:

Jul 5, 2023, 12:03:42 PM (18 months ago)

Author:

acozic

Message:

Add of an output file containing mass flow when offline parameter is to "yes"
this file will be on horizontal grid with vertical level klev
When LMDZ is coupled to Inca, we don't call anymore the routine phystoken if offline=y

Anne Cozic

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_output_write_mod.F90 (modified) (6 diffs)
• physiq_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -19 +19 @@
        ok_ade, ok_aie, ok_volcan, ivap, iliq, isol, ibs, ok_sync, &
        ptconv, read_climoz, clevSTD, ptconvth, &
-       d_u, d_t, qx, d_qx, zmasse, flag_aerosol, flag_aerosol_strat, ok_cdnc)
+       d_u, d_t, qx, d_qx, zmasse, flag_aerosol, flag_aerosol_strat, ok_cdnc, t, u1, v1)
 
     ! This subroutine does the actual writing of diagnostics that were
@@ -31 +31 @@
     USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy
     USE vertical_layers_mod, ONLY : ap, bp, aps, bps
+    USE phystokenc_mod, ONLY: offline
     USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &
          o_longitude, o_latitude, &
@@ -360 +361 @@
          ep, epmax_diag, &  ! epmax_cape
          p_tropopause, t_tropopause, z_tropopause, &
-         zxfluxt,zxfluxq
-!FC
+         zxfluxt,zxfluxq, &
+! offline 
+         da, mp, phi, wght_cvfd
+
 
 #ifdef CPP_StratAer 
@@ -458 +461 @@
     INTEGER :: flag_aerosol 
     LOGICAL :: ok_cdnc
+    REAL, DIMENSION(klon,klev) :: t   ! output for phystoken - offline flux 
+    REAL, DIMENSION(klon) :: u1, v1   ! output for phystoken - offline flux
+    
     REAL, DIMENSION(3) :: freq_moyNMC
 
@@ -482 +488 @@
 #endif
     REAL, PARAMETER :: un_jour=86400.
+    CHARACTER(len=12) :: nvar    
     INTEGER :: ISW, itr, ixt, it
     CHARACTER*1 ch1
     CHARACTER(LEN=maxlen) :: varname, dn
-
+    REAL, DIMENSION(klon,klev) :: coefh_stok 
+    
+    
 #ifdef CPP_StratAer 
     LOGICAL, PARAMETER :: debug_strataer=.FALSE.
@@ -874 +883 @@
        CALL histwrite_phy(o_topl0, toplw0)
 
+! offline 
+#ifdef CPP_XIOS
+       IF (offline) THEN
+
+          coefh_stok(:,1)      = cdragh(:)
+          coefh_stok(:,2:klev) = coefh(:,2:klev, is_ave)
+          
+          CALL histwrite_phy('upwd_stok', upwd)
+          CALL histwrite_phy('t_stok', t)
+          CALL histwrite_phy('fm_th_stok', fm_therm(:,1:klev))
+          CALL histwrite_phy('en_th_stok', entr_therm)
+          CALL histwrite_phy('da_stok',da )
+          CALL histwrite_phy('mp_stok',mp )
+          CALL histwrite_phy('dnwd_stok', dnwd)
+          CALL histwrite_phy('wght_stok', wght_cvfd)
+          CALL histwrite_phy('coefh_stok', coefh_stok)
+          CALL histwrite_phy('yu1_stok', u1) 
+          CALL histwrite_phy('yv1_stok', v1)
+
+          DO k=1,klev
+             IF (k<10) THEN
+                WRITE(nvar,'(i1)') k
+             ELSE IF (k<100) THEN
+                WRITE(nvar,'(i2)') k
+             ELSE
+                WRITE(nvar,'(i3)') k
+             END IF
+             nvar='phi_lev'//trim(nvar)
+             CALL histwrite_phy(nvar,phi(:,:,k))
+          END DO
+          
+       ENDIF
+#endif
+
+
+
+       
        IF (vars_defined) THEN
           zx_tmp_fi2d(:) = swup(:,klevp1)*swradcorr(:)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1834 +1834 @@
                               ptconv, read_climoz, clevSTD,                   &
                               ptconvth, d_u, d_t, qx, d_qx, zmasse,           &
-                              flag_aerosol, flag_aerosol_strat, ok_cdnc)
+                              flag_aerosol, flag_aerosol_strat, ok_cdnc, t, u1, v1)
 #endif
 
@@ -5213 +5213 @@
     !ENDIF    ! (iflag_phytrac=1)
 
+
+#ifndef INCA 
     IF (offline) THEN
 
@@ -5228 +5230 @@
 
     ENDIF
-
+#endif
     !
     ! Calculer le transport de l'eau et de l'energie (diagnostique)
@@ -5547 +5549 @@
          ok_sync, ptconv, read_climoz, clevSTD,          &
          ptconvth, d_u, d_t, qx, d_qx, zmasse,           &
-         flag_aerosol, flag_aerosol_strat, ok_cdnc)
+         flag_aerosol, flag_aerosol_strat, ok_cdnc, t, u1, v1 )
 #endif