/[MITgcm]/MITgcm/pkg/fizhi/fizhi_init_chem.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_init_chem.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.13 by molod, Thu Aug 5 20:11:11 2004 UTC revision 1.14 by molod, Wed Sep 1 18:35:33 2004 UTC
# Line 68  C*************************************** Line 68  C***************************************
68        _RL lats(nlat)        _RL lats(nlat)
69        _RL levs(nlev)        _RL levs(nlev)
70    
71        integer time        integer time0
72        integer lat        integer lat
73        integer lev        integer lev
74    
# Line 90  c ----------------------------- Line 90  c -----------------------------
90    
91  c Read Moisture Amounts by Month and Level  c Read Moisture Amounts by Month and Level
92  c ----------------------------------------  c ----------------------------------------
93        do time=1,ntime        do time0=1,ntime
94        read (ku,1001)        read (ku,1001)
95        do  lat=1,nlat        do  lat=1,nlat
96        read(ku,1000) (qz(lat,lev,time),lev=1,nlev)        read(ku,1000) (qz(lat,lev,time0),lev=1,nlev)
97        enddo        enddo
98        enddo        enddo
99    
100  c Convert from Volume Mixing Ratio to Mass Mixing Ratio  c Convert from Volume Mixing Ratio to Mass Mixing Ratio
101  c -----------------------------------------------------  c -----------------------------------------------------
102        do time = 1,ntime        do time0 = 1,ntime
103        do  lev = 1,nlev        do  lev = 1,nlev
104        do  lat = 1,nlat        do  lat = 1,nlat
105        qz(lat,lev,time) = qz(lat,lev,time)*voltomas        qz(lat,lev,time0) = qz(lat,lev,time0)*voltomas
106        enddo        enddo
107        enddo        enddo
108        enddo        enddo
# Line 136  C*************************************** Line 136  C***************************************
136        _RL lats(nlat)        _RL lats(nlat)
137        _RL levs(nlev)        _RL levs(nlev)
138    
139        integer time        integer time0
140        integer lat        integer lat
141        integer lev        integer lev
142        integer nrec        integer nrec
# Line 165  c ------------------------------------- Line 165  c -------------------------------------
165        open(ku,file='gcmo3.data',form='unformatted',access='direct',        open(ku,file='gcmo3.data',form='unformatted',access='direct',
166       .                                                 recl=nlat*4)       .                                                 recl=nlat*4)
167    
168        do time=1,ntime        do time0=1,ntime
169        do lev=1,nlev        do lev=1,nlev
170  C Note: 2 quantities in Ozone Dataset  C Note: 2 quantities in Ozone Dataset
171        nrec = lev+(time-1)*nlev*2          nrec = lev+(time0-1)*nlev*2  
172        read(ku,rec=nrec) o3        read(ku,rec=nrec) o3
173  #if defined( _BYTESWAPIO )  #if defined( _BYTESWAPIO )
174        call mds_byteswapr4(nlat,o3)        call mds_byteswapr4(nlat,o3)
175  #endif  #endif
176             do lat=1,nlat             do lat=1,nlat
177             oz(lat,nlev-lev+1,time) = o3(lat)             oz(lat,nlev-lev+1,time0) = o3(lat)
178             enddo             enddo
179        enddo        enddo
180        enddo        enddo
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22