source: LMDZ6/branches/Ocean_skin/libf/phylmd/calcul_divers.h @ 5054

Last change on this file since 5054 was 3605, checked in by lguez, 5 years ago

Merge revisions 3427:3600 of trunk into branch Ocean_skin

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