source: LMDZ4/trunk/libf/phylmd/calcul_divers.h @ 5001

Last change on this file since 5001 was 1398, checked in by musat, 14 years ago

Last corrections for CMIP5:

  • Add O3 at standard level files histmthNMC.nc
  • Add positive attribute "down" for vertical axes for all output files
  • Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
  • Correct units for vertical axes : mb instead of hPa
  • Add mass flux at the bottom of clouds
  • Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
  • Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
  • Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
  • Fields albs, albslw become alb1, alb2 in output files
  • Correct title for rugs_* fields
  • Correct units for pbase and ptop are Pa (not mb)
  • Correct ndayrain field

FH/JLD/JYG/MAF/IM

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 655 bytes
Line 
1c
2c $Header$
3c
4c
5c initialisations diverses au "debut" du mois
6c
7      IF(MOD(itap,INT(ecrit_mth/dtime)).EQ.1) THEN
8         DO i=1, klon
9          nday_rain(i)=0.
10         ENDDO !i
11c
12c surface terre
13       DO i=1, klon
14         IF(pctsrf(i,is_ter).GT.0.) THEN
15            paire_ter(i)=airephy(i)*pctsrf(i,is_ter)
16         ENDIF
17       ENDDO
18c
19      ENDIF !MOD(itap,INT(ecrit_mth)).EQ.1
20c
21      IF(MOD(itap,INT(ecrit_day/dtime)).EQ.0) THEN
22c
23cIM calcul total_rain, nday_rain
24c
25       DO i = 1, klon
26        total_rain(i)=rain_fall(i)+snow_fall(i) 
27        IF(total_rain(i).GT.0.) nday_rain(i)=nday_rain(i)+1.
28       ENDDO
29      ENDIF !itap.EQ.ecrit_mth
Note: See TracBrowser for help on using the repository browser.