Changeset 3535 for trunk/LMDZ.VENUS/libf

Timestamp:

Dec 6, 2024, 3:17:08 PM (2 weeks ago)

Author:

mlefevre

Message:

Addition of vertical eddy coeficient profile for the 1D model. ML

File:

• trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

@@ -225 +225 @@
       REAL,save,allocatable :: source(:,:)
       REAL ycoefh(klon,klev)    ! coef d'echange pour phytrac
+      REAL :: kzz_p(klev)       ! coef d'echange pour phytrac pour le 1D
       REAL yu1(klon)            ! vents dans la premiere couche U
       REAL yv1(klon)            ! vents dans la premiere couche V
@@ -1354 +1355 @@
          ENDDO
          ENDDO
+
+         if (klon.eq.1) then !For the 1D model
+            !Reading the prescribed profile from deftank
+            open(35,file='kzz_p.txt',form='formatted',status='old')
+            rewind(35)
+            DO k=1,klev
+               read(35,*) kzz_p(k)
+            ENDDO
+            close(35)
+
+            !Implementing the new profile in m2/s
+            DO k = 1, klev
+              ycoefh(1,k) = kzz_p(k)
+            ENDDO
+         endif
    
          DO iq=1, nqmax