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

Last change on this file since 5005 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
RevLine 
[1862]1!
[3065]2! $Id: calcul_divers.h 3065 2017-11-10 13:25:09Z dcugnet $
3!
4!
[1862]5! $Header$
6!
[1912]7! Initialisations diverses au tout debut
8      IF(itap.EQ.1) THEN
9         itapm1=0
[1454]10      ENDIF
11
[1912]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
[3065]19      IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN
[1912]20!        print*,'calcul nday_rain itap ',itap
[1454]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
[1912]26
27! Initialisation fin de mois
[3065]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)
[1912]30!       print*,'initialisation itapm1 ',itapm1
31      ENDIF
[2103]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
[3065]40  IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.1) THEN
[2103]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
[3065]50  IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN
[2103]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.