/[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.4 by molod, Thu Jun 10 20:53:19 2004 UTC revision 1.5 by molod, Thu Jun 10 21:50:33 2004 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4        subroutine fizhi_init_chem(mythid,nozlats,nozlevs,ozlats,ozlevs,        subroutine fizhi_init_chem(mythid,nozlats,nozlevs,ozlats,ozlevs,
5       . ozone,nwatlats,nwatlevs,watlats,watlevs,water,       . o3,nwatlats,nwatlevs,watlats,watlevs,water,
6       . Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22)       . Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22)
7  C***********************************************************************  C***********************************************************************
8  C Subroutine fizhi_init_chem - routine to read in the ozone and upper  C Subroutine fizhi_init_chem - routine to read in the ozone and upper
# Line 21  C*************************************** Line 21  C***************************************
21        implicit none        implicit none
22  #include "CPP_EEOPTIONS.h"  #include "CPP_EEOPTIONS.h"
23        integer mythid,nozlevs,nozlats,nwatlevs,nwatlats,Nrphys        integer mythid,nozlevs,nozlats,nwatlevs,nwatlats,Nrphys
24        _RL oz(nozlats,nozlevs), water(nwatlats,nwatlevs)        _RL o3(nozlats,nozlevs,12), water(nwatlats,nwatlevs,12)
25        _RL ozlats(nozlats), ozlevs(nozlevs)        _RL ozlats(nozlats), ozlevs(nozlevs)
26        _RL watlats(nwatlats), watlevs(nwatlevs)        _RL watlats(nwatlats), watlevs(nwatlevs)
27        _RL pressure(Nrphys),methane(Nrphys),n2o(Nrphys)        _RL pressure(Nrphys),methane(Nrphys),n2o(Nrphys)
# Line 169  c ------------------------------------- Line 169  c -------------------------------------
169    
170        do time=1,ntime        do time=1,ntime
171        do lev=1,nlev        do lev=1,nlev
172        nrec = lev+(time-1)*nlev*2   ! Note: 2 quantities in Ozone Dataset  C Note: 2 quantities in Ozone Dataset
173          nrec = lev+(time-1)*nlev*2  
174        read(ku,rec=nrec) o3        read(ku,rec=nrec) o3
175             do lat=1,nlat             do lat=1,nlat
176             oz(lat,nlev-lev+1,time) = o3(lat)             oz(lat,nlev-lev+1,time) = o3(lat)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22