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

Diff of /trunk/dyn3d/dynetat0.f

Parent Directory | Revision Log | View Patch Patch

-revision 130 by guez,
Tue Feb 24 15:43:51 2015 UTC
+revision 259 by guez,
Tue Mar  6 16:19:52 2018 UTC
 Line 1
  module dynetat0_m
+   use dimens_m, only: iim, jjm
    IMPLICIT NONE
+   private iim, jjm
    INTEGER day_ini
    ! day number at the beginning of the run, based at value 1 on
    ! January 1st of annee_ref
-   integer:: day_ref = 1 ! jour de l'année de l'état initial
+   integer:: day_ref = 1 ! jour de l'ann\'ee de l'\'etat initial
-   ! (= 350 si 20 décembre par exemple)
+   ! (= 350 si 20 d\'ecembre par exemple)
    integer:: annee_ref = 1998 ! Annee de l'etat initial (avec 4 chiffres)
+   REAL clon ! longitude of the center of the zoom, in rad
+   real clat ! latitude of the center of the zoom, in rad
+   real grossismx, grossismy
+   ! facteurs de grossissement du zoom, selon la longitude et la latitude
+   ! = 2 si 2 fois, = 3 si 3 fois, etc.
+   real dzoomx, dzoomy
+   ! extensions en longitude et latitude de la zone du zoom (fractions
+   ! de la zone totale)
+   real taux, tauy
+   ! raideur de la transition de l'int\'erieur \`a l'ext\'erieur du zoom
+   real rlatu(jjm + 1)
+   ! latitudes of points of the "scalar" and "u" grid, in rad
+   real rlatv(jjm)
+   ! latitudes of points of the "v" grid, in rad, in decreasing order
+   real rlonu(iim + 1) ! longitudes of points of the "u" grid, in rad
+   real rlonv(iim + 1)
+   ! longitudes of points of the "scalar" and "v" grid, in rad
+   real xprimu(iim + 1), xprimv(iim + 1)
+   ! 2 pi / iim * (derivative of the longitudinal zoom function)(rlon[uv])
+   REAL xprimm025(iim + 1), xprimp025(iim + 1)
+   REAL rlatu1(jjm), rlatu2(jjm), yprimu1(jjm), yprimu2(jjm)
+   REAL ang0, etot0, ptot0, ztot0, stot0
+   save
  contains
    SUBROUTINE dynetat0(vcov, ucov, teta, q, masse, ps, phis)
-Line 20 
 contains
+Line 58 
 contains
      ! This procedure reads the initial state of the atmosphere.
      use comconst, only: dtvr
-     use comgeom, only: rlonu, rlatu, rlonv, rlatv, cu_2d, cv_2d, aire_2d
      use conf_gcm_m, only: raz_date
      use dimens_m, only: iim, jjm, llm, nqmx
      use disvert_m, only: pa
-     use ener, only: etot0, ang0, ptot0, stot0, ztot0
      use iniadvtrac_m, only: tname
      use netcdf, only: NF90_NOWRITE, NF90_NOERR
      use netcdf95, only: NF95_GET_VAR, nf95_open, nf95_inq_varid, NF95_CLOSE, &
           NF95_Gw_VAR
      use nr_util, only: assert
-     use serre, only: clon, clat, grossismy, grossismx, dzoomx, dzoomy, taux, &
-          tauy
      use temps, only: itau_dyn
      use unit_nml_m, only: unit_nml
-Line 45 
 contains
+Line 79 
 contains
      ! Local variables:
      INTEGER iq
-     REAL, pointer:: tab_cntrl(:) ! tableau des paramètres du run
+     REAL, allocatable:: tab_cntrl(:) ! tableau des param\`etres du run
      INTEGER ierr, ncid, varid
      namelist /dynetat0_nml/ day_ref, annee_ref
-Line 63 
 contains
+Line 97 
 contains
           size(masse, 3)/) == llm, "dynetat0 llm")
      call assert(size(q, 4) == nqmx, "dynetat0 q nqmx")
-     ! Fichier état initial :
+     ! Fichier \'etat initial :
      call nf95_open("start.nc", NF90_NOWRITE, ncid)
      call nf95_inq_varid(ncid, "controle", varid)
-Line 102 
 contains
+Line 136 
 contains
      write(unit_nml, nml=dynetat0_nml)
      if (raz_date) then
-        print *, 'On réinitialise à la date lue dans la namelist.'
+        print *, 'Resetting the date, using the namelist.'
         day_ini = day_ref
         itau_dyn = 0
      else
-Line 114 
 contains
+Line 148 
 contains
      print *, "day_ini = ", day_ini
-     deallocate(tab_cntrl) ! pointer
      call NF95_INQ_VARID (ncid, "rlonu", varid)
      call NF95_GET_VAR(ncid, varid, rlonu)
-Line 128 
 contains
+Line 160 
 contains
      call NF95_INQ_VARID (ncid, "rlatv", varid)
      call NF95_GET_VAR(ncid, varid, rlatv)
-     call NF95_INQ_VARID (ncid, "cu", varid)
+     CALL nf95_inq_varid(ncid, 'xprimu', varid)
-     call NF95_GET_VAR(ncid, varid, cu_2d)
+     CALL nf95_get_var(ncid, varid, xprimu)
+     CALL nf95_inq_varid(ncid, 'xprimv', varid)
+     CALL nf95_get_var(ncid, varid, xprimv)
+     CALL nf95_inq_varid(ncid, 'xprimm025', varid)
+     CALL nf95_get_var(ncid, varid, xprimm025)
+     CALL nf95_inq_varid(ncid, 'xprimp025', varid)
+     CALL nf95_get_var(ncid, varid, xprimp025)
+     call NF95_INQ_VARID (ncid, "rlatu1", varid)
+     call NF95_GET_VAR(ncid, varid, rlatu1)
-     call NF95_INQ_VARID (ncid, "cv", varid)
+     call NF95_INQ_VARID (ncid, "rlatu2", varid)
-     call NF95_GET_VAR(ncid, varid, cv_2d)
+     call NF95_GET_VAR(ncid, varid, rlatu2)
-     call NF95_INQ_VARID (ncid, "aire", varid)
+     CALL nf95_inq_varid(ncid, 'yprimu1', varid)
-     call NF95_GET_VAR(ncid, varid, aire_2d)
+     CALL nf95_get_var(ncid, varid, yprimu1)
-     call NF95_INQ_VARID (ncid, "phisinit", varid)
+     CALL nf95_inq_varid(ncid, 'yprimu2', varid)
+     CALL nf95_get_var(ncid, varid, yprimu2)
+     call NF95_INQ_VARID (ncid, "phis", varid)
      call NF95_GET_VAR(ncid, varid, phis)
      call NF95_INQ_VARID (ncid, "ucov", varid)
-Line 151 
 contains
+Line 198 
 contains
      DO iq = 1, nqmx
         call NF95_INQ_VARID(ncid, tname(iq), varid, ierr)
-        IF (ierr /= NF90_NOERR) THEN
+        IF (ierr == NF90_NOERR) THEN
+           call NF95_GET_VAR(ncid, varid, q(:, :, :, iq))
+        ELSE
            PRINT *, 'dynetat0: "' // tname(iq) // '" not found, ' // &
                 "setting it to zero..."
            q(:, :, :, iq) = 0.
-        ELSE
-           call NF95_GET_VAR(ncid, varid, q(:, :, :, iq))
         ENDIF
      ENDDO

 Legend:



Removed from v.130
 


changed lines


 
Added in v.259
 Legend:



Removed from v.130
 


changed lines


 
Added in v.259
-Removed from v.130
+Added in v.259

	ViewVC Help
Powered by ViewVC 1.1.21