source: LMDZ6/branches/blowing_snow/libf/phylmd/ini_undefSTD.F90

Last change on this file was 3435, checked in by Laurent Fairhead, 6 years ago

"Historic" :-) commit merging the physics branch used for DYNAMICO with the LMDZ trunk.
The same physics branch can now be used seamlessly with the traditional lon-lat LMDZ
dynamical core and DYNAMICO.
Testing consisted in running a lon-lat LMDZ bucket simulation with the NPv6.1 physics package
with the original trunk sources and the merged sources. Tests were succesful in the sense that
numeric continuity was preserved in the restart files from both simulation. Further tests
included running both versions of the physics codes for one year in a LMDZOR setting in which
the restart files also came out identical.

Caution:

  • as the physics package now manages unstructured grids, grid information needs to be transmitted

to the surface scheme ORCHIDEE. This means that the interface defined in surf_land_orchidee_mod.F90
is only compatible with ORCHIDEE version orchidee2.1 and later versions. If previous versions of
ORCHIDEE need to be used, the CPP key ORCHIDEE_NOUNSTRUCT needs to be set at compilation time.
This is done automatically if makelmdz/makelmdz_fcm are called with the veget orchidee2.0 switch

  • due to a limitation in XIOS, the time at which limit conditions will be read in by DYNAMICO will be

delayed by one physic timestep with respect to the time it is read in by the lon-lat model. This is caused
by the line

IF (MOD(itime-1, lmt_pas) == 0 .OR. (jour_lu /= jour .AND. grid_type /= unstructured)) THEN ! time to read

in limit_read_mod.F90

Work still needed on COSP integration and XML files for DYNAMICO

EM, YM, LF
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.4 KB
Line 
1
2! $Id: ini_undefSTD.F90 3435 2019-01-22 15:21:59Z aborella $
3
4
5SUBROUTINE ini_undefstd(itap, itapm1)
6  USE dimphy
7  USE phys_state_var_mod ! Variables sauvegardees de la physique
8  USE phys_cal_mod, ONLY: mth_len
9  IMPLICIT NONE
10  include "clesphys.h"
11
12  ! ====================================================================
13
14  ! I. Musat : 09.2004
15
16  ! Initialisation - a des frequences differentes :
17
18  ! 1) des variables moyennees sur la journee "day" ou sur le mois "mth"
19  ! calculees a partir des valeurs "instantannees" de la physique
20
21  ! 2) des variables moyennes mensuelles "NMC" calculees a partir des val.
22  ! toutes les 6 heures
23
24  ! nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
25  ! ! de la physique
26  ! nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
27  ! ! de la physique
28  ! nout=3 !var. mensuelle "NMC" moyennee toutes les 6heures
29
30  ! NB: mettre "inst(X)" dans le write_hist*NMC.h !
31  ! ====================================================================
32
33  ! ym #include "dimphy.h"
34  ! variables Input/Output
35  ! INTEGER nlevSTD, klevSTD, itap
36  INTEGER itap, itapm1
37  ! PARAMETER(klevSTD=17)
38  ! REAL dtime
39
40  ! variables locales
41  ! INTEGER i, k, nout, n
42  INTEGER i, k, n
43  ! PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
44  ! REAL freq_outNMC(nout)
45  REAL un_jour
46  PARAMETER (un_jour=86400.)
47
48  ! variables Output
49  ! REAL tnondef(klon,klevSTD,nout)
50  ! REAL tsumSTD(klon,klevSTD,nout)
51
52  DO n = 1, nout
53
54    ! initialisation variables en debut du mois, de la journee ou des 6h
55
56    IF (n==1 .AND. itap-itapm1==1 .OR. n>1 .AND. mod(itap,nint( &
57        freq_outnmc(n)/phys_tstep))==1) THEN
58      ! print*,'ini_undefSTD n itap',n,itap
59      DO k = 1, nlevstd
60        DO i = 1, klon
61          tnondef(i, k, n) = 0.
62          tsumstd(i, k, n) = 0.
63          usumstd(i, k, n) = 0.
64          vsumstd(i, k, n) = 0.
65          wsumstd(i, k, n) = 0.
66          phisumstd(i, k, n) = 0.
67          qsumstd(i, k, n) = 0.
68          rhsumstd(i, k, n) = 0.
69          uvsumstd(i, k, n) = 0.
70          vqsumstd(i, k, n) = 0.
71          vtsumstd(i, k, n) = 0.
72          wqsumstd(i, k, n) = 0.
73          vphisumstd(i, k, n) = 0.
74          wtsumstd(i, k, n) = 0.
75          u2sumstd(i, k, n) = 0.
76          v2sumstd(i, k, n) = 0.
77          t2sumstd(i, k, n) = 0.
78          o3sumstd(i, k, n) = 0.
79          o3daysumstd(i, k, n) = 0.
80        END DO !i
81      END DO !k
82    END IF !
83  END DO !n
84  RETURN
85END SUBROUTINE ini_undefstd
Note: See TracBrowser for help on using the repository browser.