Changeset 4856

Timestamp:

Mar 19, 2024, 7:51:51 PM (7 weeks ago)

Author:

yann meurdesoif

Message:

• Phasing create_etat0_unstruct with regular ce0l
• Parameters for gravity wave parametrization or now computed by driver ICOSA_LMDZ

YM

File:

• LMDZ6/trunk/libf/phylmd/create_etat0_unstruct_mod.F90 (modified) (13 diffs)

LMDZ6/trunk/libf/phylmd/create_etat0_unstruct_mod.F90

@@ -1 +1 @@
 MODULE create_etat0_unstruct_mod
 
-
-
-
+  REAL, SAVE, ALLOCATABLE :: zmea_gw(:)
+  !$OMP THREADPRIVATE(zmea_gw)
+  REAL, SAVE, ALLOCATABLE :: zpic_gw(:)
+  !$OMP THREADPRIVATE(zpic_gw)
+  REAL, SAVE, ALLOCATABLE :: zval_gw(:)
+  !$OMP THREADPRIVATE(zval_gw)
+  REAL, SAVE, ALLOCATABLE :: zstd_gw(:)
+  !$OMP THREADPRIVATE(zstd_gw)
+  REAL, SAVE, ALLOCATABLE :: zsig_gw(:)
+  !$OMP THREADPRIVATE(zsig_gw)
+  REAL, SAVE, ALLOCATABLE :: zgam_gw(:)
+  !$OMP THREADPRIVATE(zgam_gw)
+  REAL, SAVE, ALLOCATABLE :: zthe_gw(:)
+  !$OMP THREADPRIVATE(zthe_gw)
+
+  PRIVATE zmea_gw, zpic_gw, zval_gw, zstd_gw, zsig_gw, zgam_gw, zthe_gw
 
 
@@ -41 +54 @@
   
  
+  SUBROUTINE init_param_gw(zmea, zpic, zval, zstd, zsig, zgam, zthe)
+  USE dimphy
+    REAL, INTENT(IN) :: zmea(klon)  
+    REAL, INTENT(IN) :: zpic(klon)
+    REAL, INTENT(IN) :: zval(klon) 
+    REAL, INTENT(IN) :: zstd(klon) 
+    REAL, INTENT(IN) :: zsig(klon) 
+    REAL, INTENT(IN) :: zgam(klon) 
+    REAL, INTENT(IN) :: zthe(klon)
+
+    ALLOCATE(zmea_gw(klon), zpic_gw(klon), zval_gw(klon), zstd_gw(klon), zsig_gw(klon), zgam_gw(klon), zthe_gw(klon)) 
+    
+    zmea_gw(:)=zmea(:)
+    zpic_gw(:)=zpic(:)
+    zval_gw(:)=zval(:)
+    zstd_gw(:)=zstd(:)
+    zsig_gw(:)=zsig(:)
+    zgam_gw(:)=zgam(:)
+    zthe_gw(:)=zthe(:)
+
+  END SUBROUTINE init_param_gw
+
+
+
+
   SUBROUTINE create_etat0_unstruct
   USE dimphy
@@ -57 +95 @@
   IMPLICIT NONE
   INCLUDE 'dimsoil.h'
+  include "clesphys.h"
 
     LOGICAL :: no_ter_antartique   ! If true, no land points are allowed at Antartic
@@ -66 +105 @@
     REAL,    DIMENSION(klon)                 :: fder
 
-    REAL,    DIMENSION(klon,nbsrf)           :: qsolsrf, snsrf
+    REAL,    DIMENSION(klon,nbsrf)           :: qsurf, snsrf
     REAL,    DIMENSION(klon,nsoilmx,nbsrf)   :: tsoil
     
@@ -74 +113 @@
     REAL,    DIMENSION(klon_mpi,nbsrf)       :: pctsrf_mpi
 
+    INCLUDE "compbl.h"
+    INCLUDE "alpale.h"
+   
     INTEGER :: ji,j,i
  
+
+!--- Initial atmospheric CO2 conc. from .def file
+    co2_ppm0 = co2_ppm
+
     IF (is_omp_master) THEN
       CALL xios_recv_field("ts",tsol_mpi)
@@ -81 +127 @@
       CALL xios_recv_field("mask",zmasq_mpi)
       IF (landice_opt .LT. 2) CALL xios_recv_field("landice",lic_mpi)
-      CALL xios_recv_field("zmea",zmea_mpi)
-      CALL xios_recv_field("zstd",zstd_mpi)
-      CALL xios_recv_field("zsig",zsig_mpi)
-      CALL xios_recv_field("zgam",zgam_mpi)
-      CALL xios_recv_field("zthe",zthe_mpi)
     ENDIF
     CALL scatter_omp(tsol_mpi,tsol)
@@ -91 +132 @@
     CALL scatter_omp(zmasq_mpi,zmasq)
     IF (landice_opt .LT. 2) CALL scatter_omp(lic_mpi,lic)
-    CALL scatter_omp(zmea_mpi,zmea)
-    CALL scatter_omp(zstd_mpi,zstd)
-    CALL scatter_omp(zsig_mpi,zsig)
-    CALL scatter_omp(zgam_mpi,zgam)
-    CALL scatter_omp(zthe_mpi,zthe)
 
     radsol(:)   = 0.0
@@ -166 +202 @@
     WHERE(pctsrf(:,is_oce)<EPSFRA) pctsrf(:,is_oce)=0.
     
-    zval(:)=max(0.,zmea-2*zstd(:))
-    zpic(:)=zmea+2*zstd(:)
+!    zval(:)=max(0.,zmea-2*zstd(:))
+!    zpic(:)=zmea+2*zstd(:)
     
 !! WARNING    DON'T FORGET FOR LATER
@@ -192 +228 @@
 !ym probably the uninitialized value was 0 for standard (regular grid) case
     falb_dif(:,:,:)=0
-
+    u10m(:,:)=0  
+    v10m(:,:)=0  
+    treedrg(:,:,:)=0
 !albedo SB <<<
     fevap(:,:) = 0.
-    DO i=1,nbsrf
-     qsolsrf(:,i)=150.
-    END DO
- 
+    qsurf = 0.
+    
     DO i=1,nbsrf 
       DO j=1,nsoilmx
@@ -207 +243 @@
     rain_fall = 0.; snow_fall = 0.
     solsw = 165.;   sollw = -53.
+    solswfdiff = 1. 
 !ym warning missing init for sollwdown => set to 0
   sollwdown  = 0.
@@ -215 +252 @@
     v_ancien=0
     q_ancien = 0.
+    ql_ancien = 0.
+    qs_ancien = 0.
+    prlw_ancien = 0.
+    prsw_ancien = 0.
+    prw_ancien = 0.
     agesno = 0.
 
+    wake_delta_pbl_TKE(:,:,:)=0
+    wake_dens(:)=0
+    awake_dens = 0.
+    cv_gen = 0.
+    ale_bl = 0.
+    ale_bl_trig =0.
+    alp_bl=0.
+    ale_wake=0.
+    ale_bl_stat=0.
+    
+    z0m(:,:)=0 ! ym missing 5th subsurface initialization
     z0m(:,is_oce) = rugmer(:)
-
-   z0m(:,is_ter) = 0.01 ! MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
-   z0m(:,is_lic) = 0.001 !MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
-
-   z0m(:,is_sic) = 0.001
-   z0h(:,:)=z0m(:,:)
+    z0m(:,is_ter) = 0.01 ! MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
+    z0m(:,is_lic) = 0.001 !MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
+    z0m(:,is_sic) = 0.001
+    z0h(:,:)=z0m(:,:)
 
     fder = 0.0
@@ -249 +300 @@
     entr_therm = 0.
     detr_therm = 0.
-    ale_bl = 0.
-    ale_bl_trig =0. 
-    alp_bl =0.
-    ale_wake = 0.
+    awake_s = 0.
+
     CALL fonte_neige_init(run_off_lic_0)
-    CALL pbl_surface_init( fder, snsrf, qsolsrf, tsoil )
+    CALL pbl_surface_init( fder, snsrf, qsurf, tsoil )
+
+    IF (iflag_pbl>1 .AND. iflag_wake>=1  .AND. iflag_pbl_split >=1) then
+     delta_tsurf = 0.
+     beta_aridity = 0.
+    END IF
+    ratqs_inter_ = 0.002
+    rneb_ancien = 0.
 
     CALL gather_omp(cell_area,cell_area_mpi)
@@ -264 +320 @@
     ENDIF
     
+    zmea(:) = zmea_gw(:)
+    zpic(:) = zpic_gw(:)
+    zval(:) = zval_gw(:)
+    zstd(:) = zstd_gw(:)
+    zsig(:) = zsig_gw(:)
+    zgam(:) = zgam_gw(:)
+    zthe(:) = zthe_gw(:)
+    DEALLOCATE(zmea_gw, zpic_gw, zval_gw, zstd_gw, zsig_gw, zgam_gw, zthe_gw)
+
     CALL phyredem( "startphy.nc" )