Changeset 443


Ignore:
Timestamp:
Jan 27, 2003, 11:47:56 AM (22 years ago)
Author:
lmdzadmin
Message:

Rajout qsol JLD
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/phyredem.F

    r433 r443  
    11c $Header$
    22c
     3cIM   SUBROUTINE phyredem (fichnom,dtime,radpas,co2_ppm,solaire,
    34      SUBROUTINE phyredem (fichnom,dtime,radpas,
    4      .           rlat,rlon, pctsrf,tsol,tsoil,deltat,qsol,snow,
     5     .           rlat,rlon, pctsrf,tsol,tsoil,deltat,qsurf,qsol,snow,
    56     .           albedo, evap, rain_fall, snow_fall,
    67     .           solsw, sollw,fder,
     
    2627      INTEGER radpas
    2728      REAL rlat(klon), rlon(klon)
     29cIM   REAL co2_ppm
     30cIM   REAL solaire
    2831      REAL tsol(klon,nbsrf)
    2932      REAL tsoil(klon,nsoilmx,nbsrf)
    3033      REAL deltat(klon)
    31       REAL qsol(klon,nbsrf)
     34      REAL qsurf(klon,nbsrf)
     35      REAL qsol(klon)
    3236      REAL snow(klon,nbsrf)
    3337      REAL albedo(klon,nbsrf)
     
    313317        ENDIF
    314318#ifdef NC_DOUBLE
    315       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,qsol(1,nsrf))
    316 #else
    317       ierr = NF_PUT_VAR_REAL (nid,nvarid,qsol(1,nsrf))
    318 #endif
    319       ENDDO
     319      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,qsurf(1,nsrf))
     320#else
     321      ierr = NF_PUT_VAR_REAL (nid,nvarid,qsurf(1,nsrf))
     322#endif
     323      END DO
     324C
     325      ierr = NF_REDEF (nid)
     326#ifdef NC_DOUBLE
     327      ierr = NF_DEF_VAR (nid,"QSOL",NF_DOUBLE,1,idim2,nvarid)
     328#else
     329      ierr = NF_DEF_VAR (nid,"QSOL",NF_FLOAT,1,idim2,nvarid)
     330#endif
     331      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20,
     332     .    "Eau dans le sol (mm)")
     333      ierr = NF_ENDDEF(nid)
     334#ifdef NC_DOUBLE
     335      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,qsol)
     336#else
     337      ierr = NF_PUT_VAR_REAL (nid,nvarid,qsol)
     338#endif
    320339c
    321340      DO nsrf = 1, nbsrf
Note: See TracChangeset for help on using the changeset viewer.