/[lmdze]/trunk/Sources/phylmd/getso4fromfile.f
ViewVC logotype

Annotation of /trunk/Sources/phylmd/getso4fromfile.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 68 - (hide annotations)
Wed Nov 14 16:59:30 2012 UTC (11 years, 6 months ago) by guez
Original Path: trunk/libf/phylmd/getso4fromfile.f90
File size: 2144 byte(s) 
Split "flincom.f90" into "flinclo.f90", "flinfindcood.f90",
"flininfo.f90" and "flinopen_nozoom.f90", in directory
"IOIPSL/Flincom".

Renamed "etat0_lim" to "ce0l", as in LMDZ.

Split "readsulfate.f" into "readsulfate.f90", "readsulfate_preind.f90"
and "getso4fromfile.f90".

In etat0, renamed variable q3d to q, as in "dynredem1". Replaced calls
to Flicom procedures by calls to NetCDF95.

In leapfrog, added call to writehist.

Extracted ASCII art from "grid_noro" into a file
"grid_noro.txt". Transformed explicit-shape local arrays into
automatic arrays, so that test on values of iim and jjm is no longer
needed. Test on weight:
          IF (weight(ii, jj) /= 0.) THEN
is useless. There is already a test before:
    if (any(weight == 0.)) stop "zero weight in grid_noro"

In "aeropt", replaced duplicated lines with different values of inu by
a loop on inu.

Removed arguments of "conf_phys". Corresponding variables are now
defined in "physiq", in a namelist. In "conf_phys", read a namelist
instead of using getin.
1 guez 68 SUBROUTINE getso4fromfile (cyr, so4)
2    
3     ! Routine for reading SO4 data from files
4    
5     use dimens_m
6     use dimphy
7     use netcdf
8     CHARACTER*15 fname
9     CHARACTER*4 cyr
10    
11     CHARACTER*6 cvar
12     INTEGER START(3), COUNT(3)
13     INTEGER STATUS, NCID, VARID
14     INTEGER imth, i, j, k, ny
15     PARAMETER (ny=jjm+1)
16    
17    
18     double precision so4mth(iim, ny, klev)
19     double precision so4(iim, ny, klev, 12)
20    
21    
22     fname = 'so4.run'//cyr//'.cdf'
23    
24     write (*,*) 'reading ', fname
25     STATUS = NF_OPEN (fname, NF_NOWRITE, NCID)
26     IF (STATUS .NE. NF_NOERR) write (*,*) 'err in open ',status
27    
28     DO imth=1, 12
29     IF (imth.eq.1) THEN
30     cvar='SO4JAN'
31     ELSEIF (imth.eq.2) THEN
32     cvar='SO4FEB'
33     ELSEIF (imth.eq.3) THEN
34     cvar='SO4MAR'
35     ELSEIF (imth.eq.4) THEN
36     cvar='SO4APR'
37     ELSEIF (imth.eq.5) THEN
38     cvar='SO4MAY'
39     ELSEIF (imth.eq.6) THEN
40     cvar='SO4JUN'
41     ELSEIF (imth.eq.7) THEN
42     cvar='SO4JUL'
43     ELSEIF (imth.eq.8) THEN
44     cvar='SO4AUG'
45     ELSEIF (imth.eq.9) THEN
46     cvar='SO4SEP'
47     ELSEIF (imth.eq.10) THEN
48     cvar='SO4OCT'
49     ELSEIF (imth.eq.11) THEN
50     cvar='SO4NOV'
51     ELSEIF (imth.eq.12) THEN
52     cvar='SO4DEC'
53     ENDIF
54     start(1)=1
55     start(2)=1
56     start(3)=1
57     count(1)=iim
58     count(2)=ny
59     count(3)=klev
60     ! write(*,*) 'here i am'
61     STATUS = NF_INQ_VARID (NCID, cvar, VARID)
62     write (*,*) ncid,imth,cvar, varid
63     ! STATUS = NF_INQ_VARID (NCID, VARMONTHS(i), VARID(i))
64     IF (STATUS .NE. NF_NOERR) write (*,*) 'err in read ',status
65     STATUS = NF_GET_VARA_DOUBLE &
66     (NCID, VARID, START,COUNT, so4mth)
67     IF (STATUS .NE. NF_NOERR) write (*,*) 'err in read data',status
68    
69     DO k=1,klev
70     DO j=1,jjm+1
71     DO i=1,iim
72     IF (so4mth(i,j,k).LT.0.) then
73     write(*,*) 'this is shit'
74     write(*,*) 'so4(',i,j,k,') =',so4mth(i,j,k)
75     endif
76     so4(i,j,k,imth)=so4mth(i,j,k)
77     ! so4(i,j,k,imth)=so4mth(k,j,i)
78     ENDDO
79     ENDDO
80     ENDDO
81     ENDDO
82    
83     STATUS = NF_CLOSE(NCID)
84     END SUBROUTINE getso4fromfile
  ViewVC Help
Powered by ViewVC 1.1.21