Changeset 5458 for LMDZ6/trunk/libf/phylmd/lmdz_surf_wind.f90

Timestamp:

Dec 29, 2024, 12:41:42 AM (6 days ago)

Author:

fhourdin

Message:

Concering replay

File:

• LMDZ6/trunk/libf/phylmd/lmdz_surf_wind.f90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_surf_wind.f90

@@ -2 +2 @@
         CONTAINS 
 
-SUBROUTINE surf_wind(klon,nsrfwnd,zu10m,zv10m,sigmaw,cstar,ustar,wstar,wind10ms,probu)
+SUBROUTINE surf_wind(klon,nsurfwind,zu10m,zv10m,sigmaw,cstar,ustar,wstar,wind10ms,probu)
 
 USE lmdz_surf_wind_ini, ONLY : iflag_surf_wind
 
 IMPLICIT NONE
-INTEGER, INTENT(IN)                :: nsrfwnd, klon
+INTEGER, INTENT(IN)                :: nsurfwind, klon
 REAL, DIMENSION(klon), INTENT(IN)  :: zu10m, zv10m
 REAL, DIMENSION(klon), INTENT(IN)  :: cstar
 REAL, DIMENSION(klon), INTENT(IN)  :: sigmaw
 REAL, DIMENSION(klon), INTENT(IN)  :: ustar, wstar
-REAL, DIMENSION(klon,nsrfwnd), INTENT(OUT)         :: wind10ms, probu
+REAL, DIMENSION(klon,nsurfwind), INTENT(OUT)         :: wind10ms, probu
 
 
-REAL, DIMENSION(klon,nsrfwnd)         :: sigma_th, sigma_wk
-REAL, DIMENSION(klon,nsrfwnd)         :: xp, yp, zz
-REAL, DIMENSION(klon,nsrfwnd)         :: vwx, vwy, vw
-REAL, DIMENSION(klon,nsrfwnd)         :: vtx, vty
-REAL, DIMENSION(klon,nsrfwnd)         :: windx, windy, wind
+REAL, DIMENSION(klon,nsurfwind)         :: sigma_th, sigma_wk
+REAL, DIMENSION(klon,nsurfwind)         :: xp, yp, zz
+REAL, DIMENSION(klon,nsurfwind)         :: vwx, vwy, vw
+REAL, DIMENSION(klon,nsurfwind)         :: vtx, vty
+REAL, DIMENSION(klon,nsurfwind)         :: windx, windy, wind
 REAL, DIMENSION(klon)                 :: ubwk, vbwk      ! ubwk et vbwk sont les vitesses moyennes dans les poches
 REAL, DIMENSION(klon)                 :: weilambda, U10mMOD
@@ -30 +30 @@
 REAL    :: ktwk, ktth, kzth
 
-print*,'LLLLLLLLLLLLLLLLLLLLL nsrfwnd=',nsrfwnd
+print*,'LLLLLLLLLLLLLLLLLLLLL nsurfwind=',nsurfwind
 pi=2.*acos(0.)
 ray=7000.
@@ -37 +37 @@
 kzth=1.
 kref=3
-nwb=nsrfwnd
+nwb=nsurfwind
 
 ubwk(klon) = zu10m(klon)
@@ -53 +53 @@
 IF (iflag_surf_wind == 0) THEN
     !U10mMOD=sqrt(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
-    IF (nsrfwnd /= 1 ) THEN
-            STOP 'Si iflag_surf_wind=0, nsrfwnd=1'
+    IF (nsurfwind /= 1 ) THEN
+            STOP 'Si iflag_surf_wind=0, nsurfwind=1'
     ENDIF
     DO i=1,klon
@@ -66 +66 @@
 
     DO i=1, klon
-        DO nmc=1, nsrfwnd
+        DO nmc=1, nsurfwind
              ! Utilisation de la distribution de weibull
              !U10mMOD=sqrt(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
@@ -90 +90 @@
 
     DO i=1, klon
-        DO nmc=1, nsrfwnd
+        DO nmc=1, nsurfwind
             ! Utilisation de la distribution du vent a l interieur et a l exterieur des poches 
             call Random_number(zz)     ! tirage uniforme entre 0 et 1.
@@ -122 +122 @@
                   wind(i,nmc) = sqrt(windx(i,nmc)**2 + windy(i,nmc)**2)
                   wind10ms(i,nmc) = wind(i,nmc)
-                  probu(i,nmc) = wind(i,nmc)/nsrfwnd
+                  probu(i,nmc) = wind(i,nmc)/nsurfwind
 
             ELSE
@@ -143 +143 @@
                   wind(i,nmc) = sqrt(windx(i,nmc)**2 + windy(i,nmc)**2)
                   wind10ms(i,nmc) = wind(i,nmc)
-                  probu(i,nmc) = wind(i,nmc)/nsrfwnd
+                  probu(i,nmc) = wind(i,nmc)/nsurfwind
                   ! print*, 'wind10ms', wind10ms(i,nmc)         
             ENDIF