Changeset 4722 for LMDZ6/trunk/libf/phylmd/screenc_mod.F90

Timestamp:

Oct 9, 2023, 5:33:07 PM (12 months ago)

Author:

Laurent Fairhead

Message:

Modification by O. Torres to the cdrag routines to include different bulk formulae
to calculate cdrag coefficients over ocean as well as an iteration of that
calculation.
The iteration is controlled by flag ok_cdrag_iter which if set to FALSE by default
to converge with previous results.
The choice of bulk formulae is set with the choix_bulk parameter
The number of iterations to run is set with nit_bulk
OT, PB, CD, LF

Location:

LMDZ6/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/screenc_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/screenc_mod.F90

@@ -18 +18 @@
                          ts, qsurf, z0m, z0h, psol, &
                          ustar, testar, qstar, okri, ri1, &
-                         pref, delu, delte, delq)
+                         pref, delu, delte, delq, s_pblh, prain, tsol, pat1)
       IMPLICIT NONE
 !-----------------------------------------------------------------------
@@ -60 +60 @@
       REAL, dimension(klon), intent(in) :: speed, temp, q_zref
       REAL, intent(in) :: zref
-      REAL, dimension(klon), intent(in) :: ts, qsurf, z0m, z0h, psol
-      REAL, dimension(klon), intent(in) :: ustar, testar, qstar, ri1         
+      REAL, dimension(klon), intent(IN) :: ts
+      REAL, dimension(klon), intent(in) :: qsurf, psol
+      REAL, dimension(klon), intent(inout):: z0m, z0h
+      REAL, dimension(klon), intent(in) :: ustar, testar, qstar, ri1   
+
+      REAL, dimension(klon), intent(inout) :: s_pblh    
+      REAL, dimension(klon), intent(in) :: prain    
+      REAL, dimension(klon), intent(in) :: tsol    
+      REAL, DIMENSION(klon), INTENT(IN)    :: pat1 !pression premier lev      
 !
       REAL, dimension(klon), intent(out) :: pref, delu, delte, delq 
@@ -88 +95 @@
       CALL cdrag (knon, nsrf, &
                     speed, temp, q_zref, gref, &
-                    psol, ts, qsurf, z0m, z0h, &
+                    psol, s_pblh, ts, qsurf, z0m, z0h, &
                     zri_zero,0,                &
-                    cdram, cdrah, zri1, pref)
+                    cdram, cdrah, zri1, pref, prain, tsol, pat1)
       DO i = 1, knon
         IF(ok_prescr_ust) THEN
@@ -114 +121 @@
                          cdrm, cdrh,  okri, &
                          ri1, iri1, &
-                         pref, delm, delh, zri1)
+                         pref, delm, delh, zri1, s_pblh, prain, tsol, pat1)
       IMPLICIT NONE
 !-----------------------------------------------------------------------
@@ -156 +163 @@
       REAL, dimension(klon), intent(in) :: speed, temp, q_zref
       REAL, intent(in) :: zref
-      REAL, dimension(klon), intent(in) :: ts, qsurf, z0m, z0h, psol
+      REAL, dimension(klon), intent(in) :: ts, qsurf, psol
+      REAL, dimension(klon), intent(inout) :: z0m, z0h
       REAL, dimension(klon), intent(in) :: cdrm, cdrh, ri1         
+      REAL, dimension(klon), intent(inout) :: s_pblh    
+      REAL, dimension(klon), intent(in) :: prain    
+      REAL, dimension(klon), intent(in) :: tsol    
+      REAL, DIMENSION(klon), INTENT(IN) :: pat1 !pression premier lev      
       INTEGER, INTENT(IN)  :: iri1 ! Richardson de la 1ere couche
 !
@@ -180 +192 @@
       CALL cdrag(knon, nsrf, &
                     speed, temp, q_zref, gref, &
-                    psol, ts, qsurf, z0m, z0h, &
+                    psol, s_pblh, ts, qsurf, z0m, z0h, &
                     ri1, iri1, &
-                    cdram, cdrah, zri1, pref)
+                    cdram, cdrah, zri1, pref, prain, tsol, pat1)
       DO i = 1, knon
         delm(i) = sqrt(cdrm(i))/sqrt(cdram(i))