source: LMDZ5/branches/LMDZ5_AR5/libf/phylmd/inistats.F

Last change on this file was 1492, checked in by Laurent Fairhead, 14 years ago

Merge of development branch LMDZ5V2.0-dev r1455:r1491 into the trunk.
Validation made locally: restart files are strictly equal between the HEAD of the trunk
and r1491 of LMDZ5V2.0-dev


Synchro de la branche de développement LMDZ5V2.0-dev r1455:r1491 et de la trunk
Validation faite en local: les fichiers restart sont équivalents entre la HEAD de la trunk
et la révision r1491 de LMDZ5V2.0-dev

File size: 4.1 KB
Line 
1      subroutine inistats(ierr)
2
3      implicit none
4
5#include "dimensions.h"
6#include "paramet.h"
7#include "comgeom.h"
8#include "comvert.h"
9#include "comconst.h"
10#include "statto.h"
11#include "netcdf.inc"
12
13      integer,intent(out) :: ierr
14      integer :: nid
15      integer :: l,nsteppd
16      real, dimension(llm) ::  sig_s
17      integer :: idim_lat,idim_lon,idim_llm,idim_llmp1,idim_time
18      real, dimension(istime) :: lt
19      integer :: nvarid
20      real, dimension(llm) :: pseudoalt
21
22      write (*,*)
23      write (*,*) '                        || STATS ||'
24      write (*,*)
25      write (*,*) 'daysec',daysec
26      write (*,*) 'dtphys',dtphys
27      nsteppd=nint(daysec/dtphys)
28      write (*,*) 'nsteppd=',nsteppd
29      if (abs(float(nsteppd)-daysec/dtphys).gt.1.e-8*daysec)
30     &   stop'Dans Instat:  1jour .ne. n pas physiques'
31
32      if(mod(nsteppd,istime).ne.0)
33     &   stop'Dans Instat:  1jour .ne. n*istime pas physiques'
34
35      istats=nsteppd/istime
36      write (*,*) 'istats=',istats
37      write (*,*) 'Storing ',istime,'times per day'
38      write (*,*) 'thus every ',istats,'physical timestep '
39      write (*,*)
40
41      do l= 1, llm
42         sig_s(l)=((ap(l)+ap(l+1))/preff+bp(l)+bp(l+1))/2.
43         pseudoalt(l)=log(preff/presnivs(l))*8. 
44      enddo
45
46      ierr = NF_CREATE("stats.nc",NF_CLOBBER,nid)
47      if (ierr.ne.NF_NOERR) then
48         write (*,*) NF_STRERROR(ierr)
49         stop ""
50      endif
51
52      ierr = NF_DEF_DIM (nid, "latitude", jjp1, idim_lat)
53      ierr = NF_DEF_DIM (nid, "longitude", iip1, idim_lon)
54      ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
55      ierr = NF_DEF_DIM (nid, "llmp1", llm+1, idim_llmp1)
56      ierr = NF_DEF_DIM (nid, "Time", NF_UNLIMITED, idim_time)
57
58      ierr = NF_ENDDEF(nid)
59      call def_var_stats(nid,"Time","Time",
60     &            "hours since 0000-00-0 00:00:00",1,
61     &            idim_time,nvarid,ierr)
62! Time is initialised later by mkstats subroutine
63
64      call def_var_stats(nid,"latitude","latitude",
65     &            "degrees_north",1,idim_lat,nvarid,ierr)
66#ifdef NC_DOUBLE
67      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatu/pi*180)
68#else
69      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatu/pi*180)
70#endif
71      call def_var_stats(nid,"longitude","East longitude",
72     &            "degrees_east",1,idim_lon,nvarid,ierr)
73#ifdef NC_DOUBLE
74      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonv/pi*180)
75#else
76      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonv/pi*180)
77#endif
78
79! Niveaux verticaux, aps et bps
80      ierr = NF_REDEF (nid)
81! presnivs
82#ifdef NC_DOUBLE
83      ierr = NF_DEF_VAR (nid,"presnivs", NF_DOUBLE, 1,idim_llm,nvarid)
84#else
85      ierr = NF_DEF_VAR (nid,"presnivs", NF_FLOAT, 1,idim_llm,nvarid)
86#endif
87      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",15,
88     &                        "Vertical levels")
89      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"Pa")
90      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',4,"down")
91      ierr = NF_ENDDEF(nid)
92#ifdef NC_DOUBLE
93      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,
94     &                          presnivs(1:llm))
95#else
96      ierr = NF_PUT_VAR_REAL (nid,nvarid,
97     &                        presnivs(1:llm))
98#endif
99! Pseudo alts
100#ifdef NC_DOUBLE
101      ierr = NF_DEF_VAR (nid,"altitude", NF_DOUBLE, 1,idim_llm,nvarid)
102#else
103      ierr = NF_DEF_VAR (nid,"altitude", NF_FLOAT, 1,idim_llm,nvarid)
104#endif
105      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",8,"altitude")
106      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"km")
107      ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',2,"up")
108      ierr = NF_ENDDEF(nid)
109#ifdef NC_DOUBLE
110      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pseudoalt)
111#else
112      ierr = NF_PUT_VAR_REAL (nid,nvarid,pseudoalt)
113#endif
114!      call def_var_stats(nid,"aps","hybrid pressure at midlayers"," ",
115!     &            1,idim_llm,nvarid,ierr)
116!#ifdef NC_DOUBLE
117!      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aps)
118!#else
119!      ierr = NF_PUT_VAR_REAL (nid,nvarid,aps)
120!#endif
121
122!      call def_var_stats(nid,"bps","hybrid sigma at midlayers"," ",
123!     &            1,idim_llm,nvarid,ierr)
124!#ifdef NC_DOUBLE
125!      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bps)
126!#else
127!      ierr = NF_PUT_VAR_REAL (nid,nvarid,bps)
128!#endif
129
130      ierr=NF_CLOSE(nid)
131
132      end subroutine inistats
133
Note: See TracBrowser for help on using the repository browser.