source: LMDZ6/branches/DYNAMICO-conv/libf/phylmd/calcul_divers.h @ 3203

Last change on this file since 3203 was 3065, checked in by Laurent Fairhead, 7 years ago

Continuing LMDZ/DYNAMICO physiqs match-up. First step in replacing dtime in the physics
  • 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 Id
File size: 1.6 KB
Line 
1!
2! $Id: calcul_divers.h 3065 2017-11-10 13:25:09Z acozic $
3!
4!
5! $Header$
6!
7! Initialisations diverses au tout debut
8      IF(itap.EQ.1) THEN
9         itapm1=0
10      ENDIF
11
12! Initialisation debut de mois
13      IF(itap.EQ.itapm1+1) THEN
14        nday_rain(:)=0.
15!       print*,'initialisation mois suivants day_rain itap',itap
16      ENDIF
17
18! Calcul fin de journee : total_rain, nday_rain
19      IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN
20!        print*,'calcul nday_rain itap ',itap
21         DO i = 1, klon
22            total_rain(i)=rain_fall(i)+snow_fall(i) 
23            IF(total_rain(i).GT.0.) nday_rain(i)=nday_rain(i)+1.
24         ENDDO
25      ENDIF
26
27! Initialisation fin de mois
28      IF(MOD(itap-itapm1,NINT(mth_len*un_jour/phys_tstep)).EQ.0) THEN
29        itapm1=itapm1+NINT(mth_len*un_jour/phys_tstep)
30!       print*,'initialisation itapm1 ',itapm1
31      ENDIF
32!
33! calcul temperatures minimale et maximale moyennees sur le mois
34!
35!initialisation debut de mois ou de journee pour les fichiers mensuels
36  IF(itap.EQ.itapm1+1) THEN
37     t2m_min_mon=0.
38     t2m_max_mon=0.
39  ENDIF
40  IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.1) THEN
41     zt2m_min_mon=zt2m
42     zt2m_max_mon=zt2m
43  ENDIF
44!calcul a chaque pas de temps pour les fichiers mensuels
45     DO i = 1, klon
46        zt2m_min_mon(i)=MIN(zt2m(i),zt2m_min_mon(i))
47        zt2m_max_mon(i)=MAX(zt2m(i),zt2m_max_mon(i))
48     ENDDO
49!fin de journee
50  IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN
51   t2m_min_mon=t2m_min_mon+zt2m_min_mon
52   t2m_max_mon=t2m_max_mon+zt2m_max_mon
53  ENDIF
54!fin mois
55  IF(itap==itapm1) THEN
56   t2m_min_mon=t2m_min_mon/mth_len
57   t2m_max_mon=t2m_max_mon/mth_len
58  ENDIF
59!
Note: See TracBrowser for help on using the repository browser.