Ignore:
Timestamp:
Apr 13, 2001, 4:58:01 PM (24 years ago)
Author:
lmdzadmin
Message:

Rajout date0 dans physiq pour interface avec ORCHIDEE
Fermeture des fichiers restart par restclo dans abort_gcm
LF

Location:
LMDZ.3.3/branches/rel-LF
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/abort_gcm.F

    r2 r205  
     1c
     2c $Header
     3c
    14      SUBROUTINE abort_gcm(modname, message, ierr)
    25     
    … …  
    1619      write(*,*) 'in abort_gcm'
    1720      call histclo
     21      call restclo
    1822c     call histclo(2)
    1923c     call histclo(3)

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

    r177 r205  
    1       SUBROUTINE clmain(dtime,itap,pctsrf,
     1      SUBROUTINE clmain(dtime,itap,date0,pctsrf,
    22     .                  t,q,u,v,
    33     .                  jour, rmu0,
    … …  
    3636c dtime----input-R- interval du temps (secondes)
    3737c itap-----input-I- numero du pas de temps
     38c date0----input-R- jour initial
    3839c t--------input-R- temperature (K)
    3940c q--------input-R- vapeur d'eau (kg/kg)
    … …  
    7677c
    7778      REAL dtime
     79      real date0
    7880      integer itap
    7981      REAL t(klon,klev), q(klon,klev)
    … …  
    368370
    369371c calculer la diffusion de "q" et de "h"
    370       CALL clqh(dtime, itap, jour, debut,lafin,
     372      CALL clqh(dtime, itap, date0,jour, debut,lafin,
    371373     e          rlon, rlat, cufi, cvfi,
    372374     e          knon, nsrf, ni, pctsrf,
    … …  
    479481      RETURN
    480482      END
    481       SUBROUTINE clqh(dtime,itime, jour,debut,lafin,
     483      SUBROUTINE clqh(dtime,itime, date0,jour,debut,lafin,
    482484     e                rlon, rlat, cufi, cvfi,
    483485     e                knon, nisurf, knindex, pctsrf,
    … …  
    512514      INTEGER knon
    513515      REAL dtime              ! intervalle du temps (s)
     516      real date0
    514517      REAL u1lay(klon)        ! vitesse u de la 1ere couche (m/s)
    515518      REAL v1lay(klon)        ! vitesse v de la 1ere couche (m/s)
    … …  
    729732      ccanopy = 365.
    730733
    731       CALL interfsurf(itime, dtime, jour, rmu0,
     734      CALL interfsurf(itime, dtime, date0, jour, rmu0,
    732735     e klon, iim, jjm, nisurf, knon, knindex, pctsrf,
    733736     e rlon, rlat, cufi, cvfi,

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

    r201 r205  
    4949!############################################################################
    5050!
    51   SUBROUTINE interfsurf_hq(itime, dtime, jour, rmu0, &
     51  SUBROUTINE interfsurf_hq(itime, dtime, date0, jour, rmu0, &
    5252      & klon, iim, jjm, nisurf, knon, knindex, pctsrf, &
    5353      & rlon, rlat, cufi, cvfi,&
    … …  
    7777!   iim, jjm     nbres de pts de grille
    7878!   dtime        pas de temps de la physique (en s)
     79!   date0        jour initial
    7980!   jour         jour dans l'annee en cours,
    8081!   rmu0         cosinus de l'angle solaire zenithal
    … …  
    138139  integer, intent(IN) :: klon
    139140  real, intent(IN) :: dtime
     141  real, intent(IN) :: date0
    140142  integer, intent(IN) :: jour
    141143  real, intent(IN)    :: rmu0(klon)
    … …  
    345347      scatter= .true.
    346348      if (.not. scatter) then
    347         call interfsol(itime, klon, dtime, nisurf, knon, &
     349        call interfsol(itime, klon, dtime, date0, nisurf, knon, &
    348350     &  knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
    349351     &  debut, lafin, ok_veget, &
    … …  
    355357     &  tsol_rad, tsurf_new, alb_new, emis_new, z0_new, dflux_l, dflux_s)
    356358      else
    357         call interfsol_scat(itime, klon, dtime, nisurf, knon, &
     359        call interfsol_scat(itime, klon, dtime, date0, nisurf, knon, &
    358360     &  knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
    359361     &  debut, lafin, ok_veget, &
    … …  
    622624!#########################################################################
    623625!
    624   SUBROUTINE interfsol(itime, klon, dtime, nisurf, knon, &
     626  SUBROUTINE interfsol(itime, klon, dtime, date0, nisurf, knon, &
    625627     & knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
    626628     & debut, lafin, ok_veget, &
    … …  
    694696  integer, intent(IN) :: klon
    695697  real, intent(IN)    :: dtime
     698  real, intent(IN)    :: date0
    696699  integer, intent(IN) :: nisurf
    697700  integer, intent(IN) :: knon
    … …  
    754757  real, dimension(klon):: qsurf
    755758  real, dimension(klon):: snow, qsol
    756   real, save                 :: date0 = 0.
    757759  real, dimension(knon,2) :: albedo_out
    758760! Pb de nomenclature
    … …  
    963965!#########################################################################
    964966!
    965   SUBROUTINE interfsol_scat(itime, klon, dtime, nisurf, knon, &
     967  SUBROUTINE interfsol_scat(itime, klon, dtime, date0, nisurf, knon, &
    966968     & knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
    967969     & debut, lafin, ok_veget, &
    … …  
    10351037  integer, intent(IN) :: klon
    10361038  real, intent(IN)    :: dtime
     1039  real, intent(IN)    :: date0
    10371040  integer, intent(IN) :: nisurf
    10381041  integer, intent(IN) :: knon
    … …  
    10951098  real, dimension(klon):: qsurf
    10961099  real, dimension(klon):: snow, qsol
    1097   real, save                 :: date0 = 0.
    10981100  real, dimension(knon,2) :: albedo_out
    10991101! Pb de nomenclature
    … …  
    11921194!
    11931195  if (debut) then
     1196    if (check) write(*,*) 'debut orchidee itime - 1', itime-1,date0
    11941197    call intersurf_main (itime-1, iim, jjm+1 , knon, knindex, dtime, &
    11951198     & lrestart_read, lrestart_write, lon_sc, lat_sc, &

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

    r179 r205  
    9494ccc      PARAMETER (soil_model=.FALSE.)
    9595      logical ok_veget
    96       parameter (ok_veget = .false.)
    97 c      parameter (ok_veget = .true.)
     96      parameter (ok_veget = .true.)
     97c      parameter (ok_veget = .false.)
    9898c======================================================================
    9999c Dans les versions precedentes, l'eau liquide nuageuse utilisee dans
    … …  
    118118c
    119119      LOGICAL ok_instan ! sortir le fichier instantane
    120       PARAMETER (ok_instan=.true.)
     120      PARAMETER (ok_instan=.false.)
    121121c
    122122      LOGICAL ok_region ! sortir le fichier regional
    … …  
    499499      character*80 abort_message
    500500      logical ok_sync
    501 
     501      real date0
    502502
    503503c
    … …  
    15251525c Verifier les temperatures
    15261526c
    1527 
    15281527      CALL hgardfou(t_seri,ftsol,'debutphy')
    15291528c
    … …  
    15931592
    15941593      fder = 0.
    1595 
    1596       CALL clmain(dtime,itap,pctsrf,
     1594      date0 = day_ini
     1595
     1596      CALL clmain(dtime,itap,date0,pctsrf,
    15971597     e            t_seri,q_seri,u_seri,v_seri,
    15981598     e            julien, rmu0,
    … …  
    16071607     s            dsens, devap,
    16081608     s            ycoefh,yu1,yv1)
     1609
    16091610c
    16101611C§§§ PB
    … …  
    16331634c         evap(i) = - fluxq(i,1) ! flux d'evaporation au sol
    16341635         evap(i) = - zxfluxq(i,1) ! flux d'evaporation au sol
     1636C LF test signe flux
     1637         sens(i) = zxfluxt(i,1)
     1638         evap(i) = zxfluxq(i,1)
    16351639         fder(i) = dsens(i) + devap(i)
    16361640      ENDDO
    … …  
    17381742      CALL abort
    17391743      ENDIF
     1744
    17401745      CALL homogene(paprs, q_seri, d_q_con, u_seri,v_seri,
    17411746     .              d_u_con, d_v_con)
    17421747      DO k = 1, klev
    1743       DO i = 1, klon
     1748        DO i = 1, klon
    17441749         t_seri(i,k) = t_seri(i,k) + d_t_con(i,k)
    17451750         q_seri(i,k) = q_seri(i,k) + d_q_con(i,k)
    17461751         u_seri(i,k) = u_seri(i,k) + d_u_con(i,k)
    17471752         v_seri(i,k) = v_seri(i,k) + d_v_con(i,k)
    1748       ENDDO
     1753        ENDDO
    17491754      ENDDO
    17501755      IF (check) THEN

  • LMDZ.3.3/branches/rel-LF/makegcm

    r179 r205  
    7777#Force/couple
    7878set COUPLE=0
    79 set ocean=`grep '^      parameter (ocean' libf/phylmd/physiq.F | awk -F\' '{print $2}'`
     79set ocean=`grep '^      parameter (ocean' $LMDGCM/libf/phylmd/physiq.F | awk -F\' '{print $2}'`
    8080banner $ocean
    8181if ($ocean == couple) then
    8282  set COUPLE=1
    8383endif
    84 
    85 
     84#veget ou non
     85set veget=`grep '^      parameter (ok_veget' $LMDGCM/libf/phylmd/physiq.F | awk -F. '{print $2}'`
     86if ($veget == true) then
     87  banner orchidee
     88endif
    8689########################################################################
    8790#  Les differentes platformes reconnues
Note: See TracChangeset for help on using the changeset viewer.