Changeset 634 for LMDZ4/branches/LMDZ4_par_0/libf/phylmd/phyetat0.F

Timestamp:

May 4, 2005, 5:11:29 PM (19 years ago)

Author:

Laurent Fairhead

Message:

Modifications faites à la physique pour la rendre parallele YM
Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester
les modifs pleinement avant leurs inclusions dans le tronc principal
LF

File:

• LMDZ4/branches/LMDZ4_par_0/libf/phylmd/phyetat0.F (modified) (9 diffs)

LMDZ4/branches/LMDZ4_par_0/libf/phylmd/phyetat0.F

@@ -5 +5 @@
 c
       SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
-     .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsurf,qsol,snow,
-     .           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw,
-     .           fder,radsol,frugs,agesno,clesphy0,
-     .           zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,tabcntr0,
-     .           t_ancien,q_ancien,ancien_ok, rnebcon, ratqs,clwcon,
-     .           run_off_lic_0)
+     .            rlat_p,rlon_p, pctsrf_p, tsol_p,tsoil_p,deltat_p,
+     .           qsurf_p,qsol_p,snow_p,albe_p, alblw_p, evap_p,
+     .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,
+     .           fder_p,radsol_p,frugs_p,agesno_p,clesphy0,
+     .           zmea_p,zstd_p,zsig_p,zgam_p,zthe_p,zpic_p,zval_p,
+     .           rugsrel_p,tabcntr0,
+     .           t_ancien_p,q_ancien_p,ancien_ok, rnebcon_p, ratqs_p,
+     .           clwcon_p,run_off_lic_0_p)
+      USE dimphy,klon=>klon2,klon2=>klon,zmasq_p=>zmasq
+      USE iophy
       IMPLICIT none
 c======================================================================
@@ -17 +21 @@
 c======================================================================
 #include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimphy.h"
 #include "netcdf.inc"
 #include "indicesol.h"
@@ -27 +31 @@
       REAL dtime
       INTEGER radpas
-      REAL rlat(klon), rlon(klon)
+      REAL rlat_p(klon2), rlon_p(klon2)
       REAL co2_ppm_etat0
       REAL solaire_etat0
+      REAL tsol_p(klon2,nbsrf)
+      REAL tsoil_p(klon2,nsoilmx,nbsrf)
+      REAL deltat_p(klon2)
+      REAL qsurf_p(klon2,nbsrf)
+      REAL qsol_p(klon2)
+      REAL snow_p(klon2,nbsrf)
+      REAL albe_p(klon2,nbsrf)
+cIM BEG alblw
+      REAL alblw_p(klon2,nbsrf)
+cIM END alblw
+      REAL evap_p(klon2,nbsrf)
+      REAL radsol_p(klon2)
+      REAL rain_fall_p(klon2)
+      REAL snow_fall_p(klon2)
+      REAL sollw_p(klon2)
+      real solsw_p(klon2)
+      real fder_p(klon2)
+      REAL frugs_p(klon2,nbsrf)
+      REAL agesno_p(klon2,nbsrf)
+      REAL zmea_p(klon2)
+      REAL zstd_p(klon2)
+      REAL zsig_p(klon2)
+      REAL zgam_p(klon2)
+      REAL zthe_p(klon2)
+      REAL zpic_p(klon2)
+      REAL zval_p(klon2)
+      REAL rugsrel_p(klon2)
+      REAL pctsrf_p(klon2, nbsrf)
+      REAL run_off_lic_0_p(klon2)
+      
+      REAL t_ancien_p(klon2,klev), q_ancien_p(klon2,klev)
+      real rnebcon_p(klon2,klev),clwcon_p(klon2,klev)
+      real ratqs_p(klon2,klev)
+      LOGICAL ancien_ok
+      
+      REAL zmasq(klon)
+      REAL rlat(klon), rlon(klon)
       REAL tsol(klon,nbsrf)
       REAL tsoil(klon,nsoilmx,nbsrf)
@@ -60 +101 @@
       REAL fractint(klon)
       REAL run_off_lic_0(klon)
-
-      REAL t_ancien(klon,klev), q_ancien(klon,klev)
-      real rnebcon(klon,klev),clwcon(klon,klev),ratqs(klon,klev)
-      LOGICAL ancien_ok
+      REAL t_ancien(klon,klev)
+      REAL q_ancien(klon,klev)
+      real rnebcon(klon,klev)
+      real clwcon(klon,klev)
+      real ratqs(klon,klev)
+     
+
 
       INTEGER        longcles
@@ -78 +122 @@
       CHARACTER*7 str7
       CHARACTER*2 str2
+      real iolat(jjm+1)
 c
 c Ouvrir le fichier contenant l'etat initial:
@@ -1248 +1293 @@
       ENDIF
 c
+      clwcon=0.
       ierr = NF_INQ_VARID (nid, "CLWCON", nvarid)
       IF (ierr.NE.NF_NOERR) THEN
@@ -1270 +1316 @@
       PRINT*,'Eau liquide convective (ecart-type) clwcon:', xmin, xmax
 c
+      rnebcon=0.
       ierr = NF_INQ_VARID (nid, "RNEBCON", nvarid)
       IF (ierr.NE.NF_NOERR) THEN
@@ -1312 +1359 @@
 c Lecture ratqs
 c
+       ratqs=0.
       ierr = NF_INQ_VARID (nid, "RATQS", nvarid)
       IF (ierr.NE.NF_NOERR) THEN
@@ -1362 +1410 @@
       ierr = NF_CLOSE(nid)
 c
+
+cym  en attendant mieux
+        iolat(1)=rlat(1)
+        iolat(jjm+1)=rlat(klon)
+        do i=2,jjm
+          iolat(i)=rlat(2+(i-2)*iim)
+        enddo
+        
+        call init_iophy(iolat,rlon(2:iim+1))
+        
+      call ScatterField( rlat,rlat_p,1)
+      call ScatterField( rlon,rlon_p,1)
+      call ScatterField( tsol,tsol_p,nbsrf)
+      call ScatterField( tsoil,tsoil_p,nsoilmx*nbsrf)
+      call ScatterField( deltat,deltat_p,1)
+      call ScatterField( qsurf,qsurf_p,nbsrf)
+      call ScatterField( qsol,qsol_p,1)
+      call ScatterField( snow,snow_p,nbsrf)
+      call ScatterField( albe,albe_p,nbsrf)
+      call ScatterField( alblw,alblw_p,nbsrf)
+      call ScatterField( evap,evap_p,nbsrf)
+      call ScatterField( radsol,radsol_p,1)
+      call ScatterField( rain_fall,rain_fall_p,1)
+      call ScatterField( snow_fall,snow_fall_p,1)
+      call ScatterField( sollw,sollw_p,1)
+      call ScatterField( solsw,solsw_p,1)
+      call ScatterField( fder,fder_p,1)
+      call ScatterField( frugs,frugs_p,nbsrf)
+      call ScatterField( agesno,agesno_p,nbsrf)
+      call ScatterField( zmea,zmea_p,1)
+      call ScatterField( zstd,zstd_p,1)
+      call ScatterField( zsig,zsig_p,1)
+      call ScatterField( zgam,zgam_p,1)
+      call ScatterField( zthe,zthe_p,1)
+      call ScatterField( zpic,zpic_p,1)
+      call ScatterField( zval,zval_p,1)
+      call ScatterField( rugsrel,rugsrel_p,1)
+      call ScatterField( pctsrf,pctsrf_p, nbsrf)
+      call ScatterField( run_off_lic_0,run_off_lic_0_p,1)
+      call ScatterField( t_ancien,t_ancien_p,klev)
+      call ScatterField( q_ancien,q_ancien_p,klev)
+      call ScatterField( rnebcon,rnebcon_p,klev)
+      call ScatterField( clwcon,clwcon_p,klev)
+      call ScatterField( ratqs,ratqs_p,klev)
+      call ScatterField( zmasq,zmasq_p,1)
+      
       RETURN
       END