Changeset 4821 for LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_ini.F90

Timestamp:

Feb 14, 2024, 9:26:23 PM (5 months ago)

Author:

evignon

Message:

reecriture de la routine de sublimation de neige soufflee en vue
du debut d'un stage de M2 sur le sujet prochainement

File:

• LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_ini.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_ini.F90

@@ -3 +3 @@
 implicit none
 
-save
-   integer :: prt_level=0,lunout
-   real RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG
-   real coef_eva_bs, expo_eva_bs, fallv_bs, zeta_bs
-   real prt_bs, pbst_bs, qbst_bs
-   integer :: iflag_saltation_bs
+   real, save, protected :: RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG
+   real, save, protected :: coef_eva_bs, expo_eva_bs, fallv_bs, zeta_bs
+   real, save, protected :: prt_bs, pbst_bs, qbst_bs
+   
+   integer, save, protected :: iflag_saltation_bs
 
-   !$OMP THREADPRIVATE(prt_level, lunout)
    !$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG)
    !$OMP THREADPRIVATE(coef_eva_bs, expo_eva_bs, fallv_bs, zeta_bs)
@@ -16 +14 @@
    !$OMP THREADPRIVATE(iflag_saltation_bs)
 
-    real tbsmelt                  ! parameter to calculate melting fraction of BS sedimentation
-    parameter (tbsmelt=278.15)
+   real, save, protected :: tbsmelt=278.15     ! parameter to calculate melting fraction of BS sedimentation
+   real, save, protected :: taumeltbs0=1800.0  ! Melting time scale of blowing snow at 273.15K
+   real, save, protected :: qbsmin=1.E-10      ! Minimum blowing snow specific content
 
-    real taumeltbs0               ! Melting time scale of blowing snow at 273.15K
-    parameter (taumeltbs0=1800.0)
-
-    real qbsmin                   ! Minimum blowing snow specific content
-    parameter (qbsmin=1.E-10)
+   !$OMP THREADPRIVATE(tbsmelt, taumeltbs0, qbsmin)
 
 
-      contains
+    contains
 
-      subroutine blowing_snow_ini(prt_level_in,lunout_in, &
-                                  RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in,&
+   subroutine blowing_snow_ini(RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in,&
                                   RVTMP2_in, RTT_in,RD_in,RG_in)
 
          USE ioipsl_getin_p_mod, ONLY : getin_p
 
-         integer, intent(in) :: prt_level_in,lunout_in
          real, intent(in) :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in
          real, intent(in) ::  RVTMP2_in, RTT_in, RD_in, RG_in
 
 
-         prt_level=prt_level_in
          RG=RG_in
          RD=RD_in
@@ -49 +41 @@
          RG=RG_in
          RVTMP2=RVTMP2_in
-         lunout=lunout_in
 
 
@@ -60 +51 @@
          prt_bs= 0.0003
          CALL getin_p('prt_bs',prt_bs)
-
 
          zeta_bs= 3.