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

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

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

revision 1.16 by molod, Wed Jul 28 01:25:07 2004 UTC revision 1.19 by molod, Thu Sep 30 16:25:58 2004 UTC
# Line 34  c--------------------------------------- Line 34  c---------------------------------------
34  #include "gridalt_mapping.h"  #include "gridalt_mapping.h"
35  #include "GRID.h"  #include "GRID.h"
36  #include "PARAMS.h"  #include "PARAMS.h"
37    #ifdef ALLOW_MNC
38    #include "MNC_PARAMS.h"
39    #endif
40    
41         integer myThid         integer myThid
42    
# Line 41  c--------------------------------------- Line 44  c---------------------------------------
44         integer im1, im2, jm1, jm2, idim2, jdim2         integer im1, im2, jm1, jm2, idim2, jdim2
45         integer nymdb,nhmsb         integer nymdb,nhmsb
46         character*40 vegdata         character*40 vegdata
47           _RL pressure0(Nrphys)
48         _RL pressure(Nrphys)         _RL pressure(Nrphys)
49         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
50         _RL fracland(sNx,sNy,Nsx,Nsy)         _RL fracland(sNx,sNy,Nsx,Nsy)
# Line 54  c--------------------------------------- Line 58  c---------------------------------------
58         nymdb = nymd0         nymdb = nymd0
59         nhmsb = nhms0         nhmsb = nhms0
60    
61    #ifdef ALLOW_MNC
62           if (useMNC) then
63            call fizhi_mnc_init(myThid)
64           endif
65    #endif
66    
67         call fizhi_alarms(nymdb,nhmsb,deltaTClock)         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
68    
69         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
# Line 88  C And now fill the earth export landtype Line 98  C And now fill the earth export landtype
98    
99  C Compute pressure profile to get methane and n2o values (bottom-up)  C Compute pressure profile to get methane and n2o values (bottom-up)
100    
101         pressure(1)=1000.         pressure0(1)=1000.
102         do L = 2,Nrphys+1         do L = 2,Nrphys+1
103          pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)/100.          pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
104           enddo
105    C Now flip pressure to a top-down array to send to fizhi init chem
106           do L = 1,Nrphys+1
107            pressure(L)=pressure0(Nrphys+2-L)
108         enddo         enddo
109    
110         call fizhi_init_chem(mythid,         call fizhi_init_chem(mythid,
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22