--- trunk/libf/phylmd/phyredem.f 2008/02/27 13:16:39 3 +++ trunk/phylmd/phyredem.f90 2013/11/15 18:45:49 76 @@ -1,474 +1,400 @@ -! -! $Header: /home/cvsroot/LMDZ4/libf/phylmd/phyredem.F,v 1.3 2005/05/25 13:10:09 fairhead Exp $ -! -c - SUBROUTINE phyredem (fichnom,dtime,radpas, - . rlat,rlon, pctsrf,tsol,tsoil, -cIM "slab" ocean - . tslab,seaice, - . qsurf,qsol,snow, - . albedo, alblw, evap, rain_fall, snow_fall, - . solsw, sollw,fder, - . radsol,frugs,agesno, - . zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel, - . t_ancien, q_ancien, rnebcon, ratqs, clwcon, - . run_off_lic_0) - use dimens_m - use indicesol - use dimphy - use conf_gcm_m - use dimsoil - use temps - use clesphys - IMPLICIT none -c====================================================================== -c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818 -c Objet: Ecriture de l'etat de redemarrage pour la physique -c====================================================================== - include "netcdf.inc" -c====================================================================== - CHARACTER*(*) fichnom - REAL dtime - INTEGER radpas - REAL, intent(in):: rlat(klon), rlon(klon) - REAL tsol(klon,nbsrf) - REAL tsoil(klon,nsoilmx,nbsrf) -cIM "slab" ocean - REAL tslab(klon), seaice(klon) - REAL qsurf(klon,nbsrf) - REAL qsol(klon) - REAL snow(klon,nbsrf) - REAL albedo(klon,nbsrf) -cIM BEG - REAL alblw(klon,nbsrf) -cIM END - REAL evap(klon,nbsrf) - REAL rain_fall(klon) - REAL snow_fall(klon) - real solsw(klon) - real sollw(klon) - real fder(klon) - REAL radsol(klon) - REAL frugs(klon,nbsrf) - REAL agesno(klon,nbsrf) - REAL zmea(klon) - REAL zstd(klon) - REAL zsig(klon) - REAL zgam(klon) - REAL zthe(klon) - REAL zpic(klon) - REAL zval(klon) - REAL rugsrel(klon) - REAL pctsrf(klon, nbsrf) - REAL t_ancien(klon,klev), q_ancien(klon,klev) - real clwcon(klon,klev),rnebcon(klon,klev),ratqs(klon,klev) - REAL run_off_lic_0(klon) -c - INTEGER nid, nvarid, idim1, idim2, idim3 - INTEGER ierr - INTEGER length - PARAMETER (length=100) - REAL tab_cntrl(length) -c - INTEGER isoil, nsrf - CHARACTER*7 str7 - CHARACTER*2 str2 -c - print *, "Call sequence information: phyredem" - ierr = NF_CREATE(fichnom, NF_CLOBBER, nid) - IF (ierr.NE.NF_NOERR) THEN - write(6,*)' Pb d''ouverture du fichier '//fichnom - write(6,*)' ierr = ', ierr - STOP 1 - ENDIF -c - ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 28, - . "Fichier redemmarage physique") -c - ierr = NF_DEF_DIM (nid, "index", length, idim1) - ierr = NF_DEF_DIM (nid, "points_physiques", klon, idim2) - ierr = NF_DEF_DIM (nid, "horizon_vertical", klon*klev, idim3) -c - ierr = NF_ENDDEF(nid) -c - DO ierr = 1, length - tab_cntrl(ierr) = 0.0 - ENDDO - tab_cntrl(1) = dtime - tab_cntrl(2) = radpas - tab_cntrl(3) = co2_ppm - tab_cntrl(4) = solaire - tab_cntrl(5) = iflag_con - tab_cntrl(6) = nbapp_rad - - IF( cycle_diurne ) tab_cntrl( 7 ) = 1. - IF( soil_model ) tab_cntrl( 8 ) = 1. - IF( new_oliq ) tab_cntrl( 9 ) = 1. - IF( ok_orodr ) tab_cntrl(10 ) = 1. - IF( ok_orolf ) tab_cntrl(11 ) = 1. - - tab_cntrl(13) = day_end - tab_cntrl(14) = annee_ref - tab_cntrl(15) = itau_phy -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "controle", NF_FLOAT, 1, idim1,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 22, - . "Parametres de controle") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "longitude", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 32, - . "Longitudes de la grille physique") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,rlon) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "latitude", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 31, - . "Latitudes de la grille physique") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,rlat) -c -C PB ajout du masque terre/mer -C - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "masque", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 16, - . "masque terre mer") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zmasq) -c BP ajout des fraction de chaque sous-surface -C -C 1. fraction de terre -C - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "FTER", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21, - . "fraction de continent") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon, is_ter)) -C -C 2. Fraction de glace de terre -C - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "FLIC", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 24, - . "fraction glace de terre") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon, is_lic)) -C -C 3. fraction ocean -C - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "FOCE", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 14, - . "fraction ocean") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon, is_oce)) -C -C 4. Fraction glace de mer -C - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "FSIC", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 18, - . "fraction glace mer") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon, is_sic)) -C -C -c - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "TS"//str2, NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28, - . "Temperature de surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,tsol(1,nsrf)) - ENDDO -c - DO nsrf = 1, nbsrf - DO isoil=1, nsoilmx - IF (isoil.LE.99 .AND. nsrf.LE.99) THEN - WRITE(str7,'(i2.2,"srf",i2.2)') isoil,nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "Tsoil"//str7,NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 29, - . "Temperature du sol No."//str7) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de couches" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,tsoil(1,isoil,nsrf)) - ENDDO - ENDDO -c -cIM "slab" ocean - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "TSLAB", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 33, - . "Ecart de la SST (pour slab-ocean)") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,tslab) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "SEAICE", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 33, - . "Glace de mer kg/m2 (pour slab-ocean)") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,seaice) -c - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"QS"//str2,NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 25, - . "Humidite de surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,qsurf(1,nsrf)) - END DO -C - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"QSOL",NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20, - . "Eau dans le sol (mm)") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,qsol) -c - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"ALBE"//str2,NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23, - . "albedo de surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,albedo(1,nsrf)) - ENDDO - -cIM BEG albedo LW - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"ALBLW"//str2,NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23, - . "albedo LW de surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,alblw(1,nsrf)) - ENDDO -cIM END albedo LW -c - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"EVAP"//str2,NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28, - . "Evaporation de surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,evap(1,nsrf)) - ENDDO - -c - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"SNOW"//str2,NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 22, - . "Neige de surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,snow(1,nsrf)) - ENDDO - -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "RADS", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28, - . "Rayonnement net a la surface") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,radsol) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "solsw", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 32, - . "Rayonnement solaire a la surface") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,solsw) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "sollw", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 27, - . "Rayonnement IF a la surface") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,sollw) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "fder", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 14, - . "Derive de flux") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,fder) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "rain_f", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21, - . "precipitation liquide") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,rain_fall) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "snow_f", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20, - . "precipitation solide") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,snow_fall) -c - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"RUG"//str2,NF_FLOAT,1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23, - . "rugosite de surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,frugs(1,nsrf)) - ENDDO -c - DO nsrf = 1, nbsrf - IF (nsrf.LE.99) THEN - WRITE(str2,'(i2.2)') nsrf - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid,"AGESNO"//str2,NF_FLOAT,1,idim2 - $ ,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 15, - . "Age de la neige surface No."//str2) - ierr = NF_ENDDEF(nid) - ELSE - PRINT*, "Trop de sous-mailles" - stop 1 - ENDIF - ierr = NF_PUT_VAR_REAL (nid,nvarid,agesno(1,nsrf)) - ENDDO -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "ZMEA", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zmea) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "ZSTD", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zstd) - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "ZSIG", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zsig) - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "ZGAM", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zgam) - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "ZTHE", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zthe) - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "ZPIC", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zpic) - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "ZVAL", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,zval) - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "RUGSREL", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,rugsrel) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "TANCIEN", NF_FLOAT, 1, idim3,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,t_ancien) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "QANCIEN", NF_FLOAT, 1, idim3,nvarid) - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,q_ancien) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "RUGMER", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28, - . "Longueur de rugosite sur mer") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,frugs(1,is_oce)) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "CLWCON", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28, - . "Eau liquide convective") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,clwcon) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "RNEBCON", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28, - . "Nebulosite convective") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,rnebcon) -c - ierr = NF_REDEF (nid) - ierr = NF_DEF_VAR (nid, "RATQS", NF_FLOAT, 1, idim2,nvarid) - ierr = NF_PUT_ATT_TEXT(nid,nvarid,"title", 5, - . "Ratqs") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,ratqs) -c -c run_off_lic_0 -c - ierr = NF_REDEF (nid) - ierr=NF_DEF_VAR(nid,"RUNOFFLIC0",NF_FLOAT, 1,idim2,nvarid) - ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 10, - . "Runofflic0") - ierr = NF_ENDDEF(nid) - ierr = NF_PUT_VAR_REAL (nid,nvarid,run_off_lic_0) -c -c - ierr = NF_CLOSE(nid) -c - RETURN - END +module phyredem_m + + IMPLICIT NONE + +contains + + SUBROUTINE phyredem(fichnom, rlat, rlon, pctsrf, tsol, tsoil, tslab, & + seaice, qsurf, qsol, snow, albedo, alblw, evap, rain_fall, snow_fall, & + solsw, sollw, fder, radsol, frugs, agesno, zmea, zstd, zsig, zgam, & + zthe, zpic, zval, t_ancien, q_ancien, rnebcon, ratqs, clwcon, & + run_off_lic_0, sig1, w01) + + ! From phylmd/phyredem.F, version 1.3 2005/05/25 13:10:09 + ! Author: Z. X. Li (LMD/CNRS) + ! Date: 1993/08/18 + ! Objet : écriture de l'état de démarrage ou redémarrage pour la physique + + USE dimphy, ONLY: klev, klon, zmasq + USE dimsoil, ONLY: nsoilmx + USE indicesol, ONLY: is_lic, is_oce, is_sic, is_ter, nbsrf + USE netcdf, ONLY: nf90_clobber, nf90_global, nf90_float + USE netcdf95, ONLY: nf95_create, nf95_put_att, nf95_def_dim, & + nf95_def_var, nf95_enddef, nf95_redef, nf95_put_var, nf95_close + USE temps, ONLY: itau_phy + + CHARACTER(len=*) fichnom + REAL, INTENT(IN):: rlat(klon), rlon(klon) + REAL, INTENT(IN):: pctsrf(klon, nbsrf) + REAL tsol(klon, nbsrf) + REAL tsoil(klon, nsoilmx, nbsrf) + REAL tslab(klon), seaice(klon) !IM "slab" ocean + REAL qsurf(klon, nbsrf) + REAL, intent(in):: qsol(klon) + REAL snow(klon, nbsrf) + REAL albedo(klon, nbsrf) + REAL alblw(klon, nbsrf) + REAL evap(klon, nbsrf) + REAL, INTENT(IN):: rain_fall(klon) + REAL snow_fall(klon) + REAL solsw(klon) + REAL, INTENT(IN):: sollw(klon) + REAL fder(klon) + REAL radsol(klon) + REAL frugs(klon, nbsrf) + REAL agesno(klon, nbsrf) + REAL zmea(klon) + REAL, intent(in):: zstd(klon) + REAL, intent(in):: zsig(klon) + REAL zgam(klon) + REAL zthe(klon) + REAL zpic(klon) + REAL zval(klon) + REAL t_ancien(klon, klev), q_ancien(klon, klev) + REAL rnebcon(klon, klev), ratqs(klon, klev), clwcon(klon, klev) + REAL run_off_lic_0(klon) + real, intent(in):: sig1(klon, klev) ! section adiabatic updraft + + real, intent(in):: w01(klon, klev) + ! vertical velocity within adiabatic updraft + + ! Local: + + INTEGER ncid, idim2, idim3 + integer varid, varid_run_off_lic_0, varid_sig1, varid_w01, varid_rlon + integer varid_rlat, varid_zmasq, varid_fter, varid_flic, varid_foce + integer varid_fsic + INTEGER isoil, nsrf + CHARACTER(len=7) str7 + CHARACTER(len=2) str2 + + !------------------------------------------------------------ + + PRINT *, 'Call sequence information: phyredem' + CALL nf95_create(fichnom, nf90_clobber, ncid) + + call nf95_put_att(ncid, nf90_global, 'title', & + 'Fichier redémarrage physique') + call nf95_put_att(ncid, nf90_global, "itau_phy", itau_phy) + + call nf95_def_dim(ncid, 'points_physiques', klon, idim2) + call nf95_def_dim(ncid, 'klev', klev, idim3) + + call nf95_def_var(ncid, 'longitude', nf90_float, idim2, varid_rlon) + call nf95_def_var(ncid, 'latitude', nf90_float, idim2, varid_rlat) + + call nf95_def_var(ncid, 'masque', nf90_float, idim2, varid_zmasq) + call nf95_put_att(ncid, varid_zmasq, 'title', 'masque terre mer') + + ! Fractions de chaque sous-surface + + call nf95_def_var(ncid, 'FTER', nf90_float, idim2, varid_fter) + call nf95_put_att(ncid, varid_fter, 'title', 'fraction de continent') + + call nf95_def_var(ncid, 'FLIC', nf90_float, idim2, varid_flic) + call nf95_put_att(ncid, varid_flic, 'title', 'fraction glace de terre') + + call nf95_def_var(ncid, 'FOCE', nf90_float, idim2, varid_foce) + call nf95_put_att(ncid, varid_foce, 'title', 'fraction ocean') + + call nf95_def_var(ncid, 'FSIC', nf90_float, idim2, varid_fsic) + call nf95_put_att(ncid, varid_fsic, 'title', 'fraction glace mer') + + call nf95_enddef(ncid) + + call nf95_put_var(ncid, varid_rlon, rlon) + call nf95_put_var(ncid, varid_rlat, rlat) + call nf95_put_var(ncid, varid_zmasq, zmasq) + call nf95_put_var(ncid, varid_fter, pctsrf(:, is_ter)) + call nf95_put_var(ncid, varid_flic, pctsrf(:, is_lic)) + call nf95_put_var(ncid, varid_foce, pctsrf(:, is_oce)) + call nf95_put_var(ncid, varid_fsic, pctsrf(:, is_sic)) + + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'TS'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Temperature de surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, tsol(:, nsrf)) + END DO + + DO nsrf = 1, nbsrf + DO isoil = 1, nsoilmx + IF (isoil<=99 .AND. nsrf<=99) THEN + WRITE (str7, '(i2.2, "srf", i2.2)') isoil, nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'Tsoil'//str7, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Temperature du sol No.'//str7) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de couches' + STOP 1 + END IF + call nf95_put_var(ncid, varid, tsoil(:, isoil, nsrf)) + END DO + END DO + + !IM "slab" ocean + call nf95_redef(ncid) + call nf95_def_var(ncid, 'TSLAB', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Ecart de la SST (pour slab-ocean)') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, tslab) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'SEAICE', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Glace de mer kg/m2 (pour slab-ocean)') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, seaice) + + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'QS'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Humidite de surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, qsurf(:, nsrf)) + END DO + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'QSOL', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', 'Eau dans le sol (mm)') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, qsol) + + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ALBE'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'albedo de surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, albedo(:, nsrf)) + END DO + + !IM BEG albedo LW + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ALBLW'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'albedo LW de surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, alblw(:, nsrf)) + END DO + !IM END albedo LW + + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'EVAP'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Evaporation de surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, evap(:, nsrf)) + END DO + + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'SNOW'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Neige de surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, snow(:, nsrf)) + END DO + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'RADS', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Rayonnement net a la surface') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, radsol) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'solsw', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Rayonnement solaire a la surface') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, solsw) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'sollw', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Rayonnement IF a la surface') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, sollw) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'fder', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', 'Derive de flux') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, fder) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'rain_f', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', 'precipitation liquide') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, rain_fall) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'snow_f', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', 'precipitation solide') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, snow_fall) + + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'RUG'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'rugosite de surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, frugs(:, nsrf)) + END DO + + DO nsrf = 1, nbsrf + IF (nsrf<=99) THEN + WRITE (str2, '(i2.2)') nsrf + call nf95_redef(ncid) + call nf95_def_var(ncid, 'AGESNO'//str2, nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Age de la neige surface No.'//str2) + call nf95_enddef(ncid) + ELSE + PRINT *, 'Trop de sous-mailles' + STOP 1 + END IF + call nf95_put_var(ncid, varid, agesno(:, nsrf)) + END DO + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ZMEA', nf90_float, idim2, varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, zmea) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ZSTD', nf90_float, idim2, varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, zstd) + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ZSIG', nf90_float, idim2, varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, zsig) + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ZGAM', nf90_float, idim2, varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, zgam) + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ZTHE', nf90_float, idim2, varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, zthe) + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ZPIC', nf90_float, idim2, varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, zpic) + call nf95_redef(ncid) + call nf95_def_var(ncid, 'ZVAL', nf90_float, idim2, varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, zval) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'TANCIEN', nf90_float, (/idim2, idim3/), varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, t_ancien) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'QANCIEN', nf90_float, (/idim2, idim3/), varid) + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, q_ancien) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'RUGMER', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', & + 'Longueur de rugosite sur mer') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, frugs(:, is_oce)) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'CLWCON', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', 'Eau liquide convective') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, clwcon(:, 1)) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'RNEBCON', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', 'Nebulosite convective') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, rnebcon(:, 1)) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'RATQS', nf90_float, idim2, varid) + call nf95_put_att(ncid, varid, 'title', 'Ratqs') + call nf95_enddef(ncid) + call nf95_put_var(ncid, varid, ratqs(:, 1)) + + call nf95_redef(ncid) + call nf95_def_var(ncid, 'RUNOFFLIC0', nf90_float, idim2, & + varid_run_off_lic_0) + call nf95_put_att(ncid, varid_run_off_lic_0, 'title', 'Runofflic0') + + call nf95_def_var(ncid, 'sig1', nf90_float, (/idim2, idim3/), varid_sig1) + call nf95_put_att(ncid, varid_sig1, 'long_name', & + 'section adiabatic updraft') + + call nf95_def_var(ncid, 'w01', nf90_float, (/idim2, idim3/), varid_w01) + call nf95_put_att(ncid, varid_w01, 'long_name', & + 'vertical velocity within adiabatic updraft') + + call nf95_enddef(ncid) + + call nf95_put_var(ncid, varid_run_off_lic_0, run_off_lic_0) + call nf95_put_var(ncid, varid_sig1, sig1) + call nf95_put_var(ncid, varid_w01, w01) + + call nf95_close(ncid) + + END SUBROUTINE phyredem + +end module phyredem_m