Changeset 3780 for LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

Timestamp:

Oct 22, 2020, 2:50:18 PM (4 years ago)

Author:

evignon

Message:

Premiere comission Etienne: changements pour le 1D (forcage en Ts au dessus des continents) et inclusion drag arbres dans yamada4_num=6

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -233 +233 @@
        CALL getin('ok_flux_surf',ok_flux_surf)
 
+!Config  Key  = ok_forc_tsurf
+!Config  Desc = forcage ou non par la Ts
+!Config  Def  = false
+!Config  Help = forcage ou non par la Ts
+       ok_forc_tsurf=.false.
+       CALL getin('ok_forc_tsurf',ok_forc_tsurf)
+
 !Config  Key  = ok_prescr_ust
 !Config  Desc = ustar impose ou non
@@ -239 +246 @@
        ok_prescr_ust = .false.
        CALL getin('ok_prescr_ust',ok_prescr_ust)
+
+
+!Config  Key  = ok_prescr_beta
+!Config  Desc = betaevap impose ou non
+!Config  Def  = false
+!Config  Help = betaevap impose ou non
+       ok_prescr_beta = .false.
+       CALL getin('ok_prescr_beta',ok_prescr_beta)
 
 !Config  Key  = ok_old_disvert
@@ -280 +295 @@
 !Config  Desc = surface temperature
 !Config  Def  = 290.
-!Config  Help = not used if type_ts_forcing=1 in lmdz1d.F
+!Config  Help = surface temperature
        tsurf = 290.
        CALL getin('tsurf',tsurf)
@@ -297 +312 @@
        zsurf = 0.
        CALL getin('zsurf',zsurf)
+! EV pour accord avec format standard       
+       CALL getin('zorog',zsurf)
+
 
 !Config  Key  = rugos
@@ -359 +377 @@
        qsolinp = 1.
        CALL getin('qsolinp',qsolinp)
+
+
+
+!Config  Key  = betaevap
+!Config  Desc = beta for actual evaporation when prescribed
+!Config  Def  = 1.0
+!Config  Help = 
+       betaevap = 1.
+       CALL getin('betaevap',betaevap)      
 
 !Config  Key  = zpicinp
@@ -520 +547 @@
        CALL getin('forc_ustar',forc_ustar)
        IF (forc_ustar .EQ. 1) ok_prescr_ust=.true.
+
 
 !Config  Key  = nudging_u
@@ -1248 +1276 @@
       END
 
-!======================================================================
-       SUBROUTINE read_tsurf1d(knon,sst_out)
-
-! This subroutine specifies the surface temperature to be used in 1D simulations
-
-      USE dimphy, ONLY : klon
-
-      INTEGER, INTENT(IN)                  :: knon     ! nomber of points on compressed grid
-      REAL, DIMENSION(klon), INTENT(OUT)   :: sst_out  ! tsurf used to force the single-column model
-
-       INTEGER :: i
-! COMMON defined in lmdz1d.F:
-       real ts_cur
-       common /sst_forcing/ts_cur
-
-       DO i = 1, knon
-        sst_out(i) = ts_cur
-       ENDDO
-
-      END SUBROUTINE read_tsurf1d
-
+!!======================================================================
+!       SUBROUTINE read_tsurf1d(knon,sst_out)
+!
+!! This subroutine specifies the surface temperature to be used in 1D simulations
+!
+!      USE dimphy, ONLY : klon
+!
+!      INTEGER, INTENT(IN)                  :: knon     ! nomber of points on compressed grid
+!      REAL, DIMENSION(klon), INTENT(OUT)   :: sst_out  ! tsurf used to force the single-column model
+!
+!       INTEGER :: i
+!! COMMON defined in lmdz1d.F:
+!       real ts_cur
+!       common /sst_forcing/ts_cur
+
+!       DO i = 1, knon
+!        sst_out(i) = ts_cur
+!       ENDDO
+!
+!      END SUBROUTINE read_tsurf1d
+!
 !===============================================================
       subroutine advect_vert(llm,w,dt,q,plev)