Changeset 437 for LMDZ.3.3

Timestamp:

Jan 24, 2003, 4:23:43 PM (22 years ago)

Author:

lmdzadmin

Message:

Introduction de qsurf dans l'etat de redemmarage JLD
IM

File:

• LMDZ.3.3/branches/rel-LF/libf/dyn3d/etat0_netcdf.F (modified) (4 diffs)

LMDZ.3.3/branches/rel-LF/libf/dyn3d/etat0_netcdf.F

@@ -35 +35 @@
       REAL :: q3d(iip1, jjp1, llm,nqmx), qsat(iip1, jjp1, llm)
       REAL :: tsol(klon), qsol(klon), sn(klon)
-      REAL :: tsolsrf(klon,nbsrf), qsolsrf(klon,nbsrf),snsrf(klon,nbsrf) 
+      REAL :: tsolsrf(klon,nbsrf), qsurf(klon,nbsrf),snsrf(klon,nbsrf) 
       REAL :: albe(klon,nbsrf), evap(klon,nbsrf)
       REAL :: tsoil(klon,nsoilmx,nbsrf) 
@@ -63 +63 @@
       CHARACTER*80 :: varname
       !
-      INTEGER :: i,j, ig, l, ji,ii1,ii2
+      INTEGER :: i,j, ig, l, ji,ii1,ii2, n
       REAL :: xpi
       !
@@ -635 +635 @@
       albe(:,is_sic) = 0.6
       evap(:,:) = 0.
-      qsolsrf(:,is_ter) = 150
-      qsolsrf(:,is_lic) = 150
-      qsolsrf(:,is_oce) = 150.
-      qsolsrf(:,is_sic) = 150.
+C     vapeur d'eau a la surface =- vapeur d'eau premier niveau
+      do n = 1, nbsrf 
+        DO j=2,jjm
+          DO i =1, iim
+            qsurf((j-2)*iim+1+i,n)=q3d(i,j,1,1)
+          END DO 
+        END DO
+        qsurf(1,n)=q3d(1,1,1,1)  ! pole nord 
+        qsurf(klon,n)=q3d(1,jjp1,1,1)  ! pole sud
+      END DO 
+      qsol(:) = 150.
       do i = 1, nbsrf
         do j = 1, nsoilmx
@@ -660 +667 @@
       ratqs = 0.0
 
-      call phyredem("startphy.nc",phystep,radpas, 
-     $    latfi, lonfi, pctsrf, tsolsrf, tsoil, deltat, qsolsrf, snsrf, 
-     $    albe, evap, rain_fall, snow_fall, solsw, sollw, fder,
+cIM   call phyredem("startphy.nc",phystep,radpas, co2_ppm, solaire,
+      call phyredem("startphy.nc",phystep,radpas,
+     $    latfi, lonfi, pctsrf, tsolsrf, tsoil, deltat, qsurf, qsol, 
+     $    snsrf, albe, evap, rain_fall, snow_fall, solsw, sollw, fder,
      $    radsol, frugs,  agesno, 
      $    zmea, zstd, zsig, zgam, zthe, zpic, zval, rugsrel,