/[lmdze]/trunk/dyn3d/gcm.f

Diff of /trunk/dyn3d/gcm.f

Parent Directory | Revision Log | View Patch Patch

-trunk/libf/dyn3d/gcm.f90
revision 24 by guez,
Wed Mar  3 13:23:49 2010 UTC
+trunk/dyn3d/gcm.f
revision 108 by guez,
Tue Sep 16 14:00:41 2014 UTC
 Line 1
  PROGRAM gcm
-   ! General circulation model of LMD
-   ! From "gcm.F", version 1.4, 2006/04/04 15:05:16
-   ! Avec coordonnées verticales hybrides, avec nouveaux opérateurs de
-   ! dissipation * (gradiv2, divgrad2, nxgraro2)
    ! Authors: P. Le Van, L. Fairhead, F. Hourdin
+   ! From "gcm.F", version 1.4, 2006/04/04 15:05:16
-   ! Possibilité de choisir le schéma pour l'advection de "q", en
+   ! General circulation model of LMD. Avec coordonnée verticale
-   ! modifiant "iadv" dans "traceur.def".
+   ! hybride, avec nouveaux opérateurs de dissipation "*" ("gradiv2",
+   ! "divgrad2", "nxgraro2"). Possibilité de choisir le schéma pour
-   ! Pour Van-Leer plus vapeur d'eau saturée : iadv(1)=4
+   ! l'advection de "q", en modifiant "iadv" dans "traceur.def".
-   ! Pour Van-Leer : iadv=10
+   use comconst, only: daysec, dtvr, iniconst
-   USE IOIPSL, only: ioconf_calendar, histclo
+   use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom
+   use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond
+   use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &
+        nday, raz_date, periodav, conf_gcm, iflag_phys
    use dimens_m, only: iim, jjm, llm, nqmx
    use dimphy, only: klon
-   use paramet_m, only: ip1jm, ip1jmp1
+   USE disvert_m, ONLY : disvert
-   use comconst, only: daysec, cpp, dtvr, g, rad, r, initialize
+   use dynetat0_m, only: dynetat0, day_ini
-   use comdissnew, only: lstardis, nitergdiv, nitergrot, niterh, tetagdiv, &
+   use dynredem0_m, only: dynredem0
-        tetagrot, tetatemp
-   use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &
-        nday, raz_date, periodav, conf_gcm
-   use logic, only: iflag_phys
-   use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv
-   use temps, only: day_ref, annee_ref, day_ini, day_end, itau_dyn
-   use com_io_dyn, only: histid, histvid, histaveid
-   use tracstoke, only: istdyn, istphy
-   use abort_gcm_m, only: abort_gcm
-   use inithist_m, only: inithist
-   use initdynav_m, only: initdynav
-   use dynetat0_m, only: dynetat0
    use grid_change, only: dyn_phy, init_dyn_phy
+   use histclo_m, only: histclo
    use iniadvtrac_m, only: iniadvtrac
+   use inidissip_m, only: inidissip
+   use inifilr_m, only: inifilr
+   use initdynav_m, only: initdynav
+   use inithist_m, only: inithist
+   use init_dynzon_m, only: init_dynzon
+   USE ioconf_calendar_m, only: ioconf_calendar
+   use jumble, only: new_unit
    use leapfrog_m, only: leapfrog
-   use dynredem0_m, only: dynredem0
+   use netcdf95, only: nf95_close
-   use clesphys2, only: read_clesphys2
+   use suphec_m, only: suphec
+   use temps, only: day_ref, annee_ref, day_end, itau_dyn
+   use tracstoke, only: istdyn, istphy
+   use unit_nml_m, only: unit_nml
+   use yoethf_m, only: yoethf
+   use write_field, only: NbField, FieldId
    IMPLICIT NONE
-   REAL zdtvr ! time step for dynamics, in s
    ! Variables dynamiques :
-   REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants
+   REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm)  ! vent covariant
-   REAL teta(ip1jmp1, llm) ! température potentielle
+   REAL teta(iim + 1, jjm + 1, llm) ! température potentielle
-   REAL q(ip1jmp1, llm, nqmx) ! champs advectés
+   REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advectés
-   REAL ps(ip1jmp1) ! pression au sol (Pa)
+   REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)
+   REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
-   REAL masse(ip1jmp1, llm) ! masse d'air
    REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol
    ! Variables pour le fichier histoire :
-   REAL time_0
+   REAL time_0 ! time in day, as a fraction of day, in [0, 1[
-   !!INTEGER i
    ! Calendrier :
    LOGICAL:: true_calendar = .false. ! default value
-   ! Variables pour l'initialisation de la physique :
-   REAL zcufi(klon), zcvfi(klon) ! "cu" and "cv" values on the scalar grid
-   REAL latfi(klon), lonfi(klon)
-   REAL airefi(klon)
    logical mask_v(iim + 1, jjm)
    ! (mask for points in the "v" grid, first index is for longitude,
    ! second index is for latitude)
+   integer i
    namelist /main_nml/true_calendar
    !------------------------------------------------------------
+   call new_unit(unit_nml)
+   open(unit_nml, file="used_namelists.txt", status="replace", action="write")
+   CALL conf_gcm
    print *, "Enter namelist 'main_nml'."
    read (unit=*, nml=main_nml)
-   write(unit=*, nml=main_nml)
+   write(unit_nml, nml=main_nml)
-   ! Initialisations:
-   call initialize
    ! Choix du calendrier :
    if (true_calendar) then
-Line 87 
 PROGRAM gcm
+Line 79 
 PROGRAM gcm
       call ioconf_calendar('360d')
    endif
-   ! Lecture des fichiers "gcm.def" ou "run.def" :
-   call read_clesphys2
-   CALL conf_gcm
    ! Initialisation des traceurs
    ! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA
    call iniadvtrac
+   CALL iniconst
    ! Lecture du fichier "start.nc" :
    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)
-   ! Lecture des paramètres de contrôle pour la simulation :
+   ! On remet le calendrier à zéro si demandé :
-   ! on recalcule éventuellement le pas de temps
+   if (raz_date) then
-   IF (MOD(day_step, iperiod) /= 0) THEN
+      print *, 'On réinitialise à la date lue dans la namelist.'
-      call abort_gcm(modname = "gcm", message = &
+      annee_ref = anneeref
-           'Il faut choisir un nombre de pas par jour multiple de "iperiod".', &
+      day_ref = dayref
-           ierr = 1)
+      day_ini = dayref
-   ENDIF
+      itau_dyn = 0
+      time_0 = 0.
-   IF (MOD(day_step,iphysiq)/=0) THEN
+   else
-      call abort_gcm(modname = "gcm", message = &
+      print *, 'On garde les dates du fichier "start".'
-           'Il faut choisir un nombre de pas par jour multiple de "iphysiq".', &
-           ierr = 1)
-   ENDIF
-   ! On remet le calendrier à zero si demandé:
-   if (annee_ref /= anneeref .or. day_ref /= dayref) then
-      print *, 'Attention : les dates initiales lues dans le fichier ' // &
-           '"start" ne correspondent pas à celles lues dans "gcm.def".'
-      if (raz_date) then
-         print *, 'On réinitialise à la date lue dans "gcm.def".'
-         annee_ref = anneeref
-         day_ref = dayref
-         day_ini = dayref
-         itau_dyn = 0
-         time_0 = 0.
-      else
-         print *, 'On garde les dates du fichier "start".'
-      endif
-   ELSE
-      raz_date = .false.
    endif
-   ! Initialisation des constantes dynamiques :
+   CALL disvert
-   zdtvr = daysec / REAL(day_step)
+   CALL inigeom ! initialisation de la géometrie
-   IF (dtvr /= zdtvr) THEN
+   CALL inifilr ! initialisation du filtre
-      print *, 'Warning: the time steps in the ".def" file and in ' // &
+   CALL inidissip
-           '"start.nc" are different'
-      print *, 'dtvr (from "start.nc") = ', dtvr
-      print *, 'zdtvr (from ".def") = ', zdtvr
-      print *, 'Using the value from the ".def" file.'
-      dtvr = zdtvr
-   ENDIF
-   CALL iniconst
-   ! Initialisation de la géometrie :
-   CALL inigeom
-   ! Initialisation du filtre :
-   CALL inifilr
-   ! Initialisation de la dissipation :
-   CALL inidissip(lstardis, nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, &
-        tetatemp)
    call init_dyn_phy
    ! Initialisation de la physique :
    IF (iflag_phys == 1) THEN
-      latfi(1)=rlatu(1)
+      rlatd(1)=rlatu(1)
-      latfi(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
+      rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
-      latfi(klon)= rlatu(jjm + 1)
+      rlatd(klon)= rlatu(jjm + 1)
-      lonfi(1)=0.
+      rlond(1)=0.
-      lonfi(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
+      rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
-      lonfi(klon)= 0.
+      rlond(klon)= 0.
-      zcufi = pack(cu_2d, dyn_phy)
+      cuphy = pack(cu_2d, dyn_phy)
       ! Construct a mask for points in the "v" grid:
       mask_v = .true.
       mask_v(2:, 1) = .false.
       mask_v(iim + 1, 2:) = .false.
-      zcvfi(:klon - 1) = pack(cv_2d, mask_v)
+      cvphy(:klon - 1) = pack(cv_2d, mask_v)
-      zcvfi(klon) = cv_2d(1, jjm)
+      cvphy(klon) = cv_2d(1, jjm)
-      ! (that value of "cv_2d" is used twice in "zcvfi")
+      ! (that value of "cv_2d" is used twice in "cvphy")
-      airefi = pack(aire_2d, dyn_phy)
+      airephy = pack(aire_2d, dyn_phy)
-      print *, 'Attention : vitesse verticale nulle dans la physique.'
+      CALL suphec
-      CALL iniphysiq(klon, llm, latfi, lonfi, airefi, &
+      call yoethf
-           zcufi, zcvfi, rad, g, r, cpp)
    ENDIF
    ! Initialisation des entrées-sorties :
-Line 187 
 PROGRAM gcm
+Line 138 
 PROGRAM gcm
    print *, "day_end = ", day_end
    CALL dynredem0("restart.nc", day_end, phis)
-   CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, &
+   CALL inithist(day_ref, annee_ref, dtvr, nqmx, t_ops = iecri * daysec, &
-        infile="dyn_hist.nc", t_ops = iecri * daysec, t_wrt = iecri * daysec)
+        t_wrt = iecri * daysec)
-   CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, &
+   CALL initdynav(day_ref, annee_ref, dtvr, nqmx, t_ops = iperiod * dtvr, &
-        infile='dyn_hist_ave.nc', t_ops = iperiod * zdtvr, &
         t_wrt = periodav * daysec)
+   call init_dynzon(dt_app = dtvr * iperiod)
    ! Choix des fréquences de stockage pour le hors-ligne :
-   istdyn = day_step / 4     ! stockage toutes les 6 h = 1 jour / 4
+   istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
    istphy = istdyn / iphysiq
    ! Intégration temporelle du modèle :
    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)
+   close(unit_nml)
    call histclo
+   do i = 1, nbfield
+      call nf95_close(FieldId(i))
+   end do
    print *, 'Simulation finished'
    print *, 'Everything is cool'

 Legend:



Removed from v.24
 


changed lines


 
Added in v.108
 Legend:



Removed from v.24
 


changed lines


 
Added in v.108
-Removed from v.24
+Added in v.108

	ViewVC Help
Powered by ViewVC 1.1.21