/[MITgcm]/MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F
ViewVC logotype

Contents of /MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F

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


Revision 1.3 - (show annotations) (download)
Tue May 24 14:57:00 2005 UTC (15 years, 1 month ago) by molod
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint58u_post, checkpoint58w_post, checkpoint57m_post, checkpoint57s_post, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint57k_post, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint58r_post, checkpoint57i_post, checkpoint57y_post, checkpoint58n_post, checkpoint58x_post, checkpoint58t_post, checkpoint58h_post, checkpoint57y_pre, checkpoint58q_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint57x_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint58k_post, checkpoint58v_post, checkpoint58s_post, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint58p_post, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint57j_post, checkpoint61z, checkpoint61x, checkpoint61y, checkpoint58b_post, checkpoint58m_post, checkpoint57l_post
Changes since 1.2: +91 -18 lines
Replace read of topog variance with zero-out of array - not needed for this experiment

1 C $Header: /u/gcmpack/MITgcm/verification/fizhi-cs-aqualev10/code/fizhi_init_fixed.F,v 1.1 2005/05/24 14:54:57 molod Exp $
2 C $Name: $
3
4 #include "FIZHI_OPTIONS.h"
5 subroutine fizhi_init_fixed (myThid)
6 c-----------------------------------------------------------------------
7 c Routine to initialise the fizhi package.
8 c
9 c Input: myThid - Process number calling this routine
10 c
11 c Notes:
12 c 1) This routine is the interface to read input datasets and set
13 c other fixed variables for fizhi
14 c the datasets are:
15 c vegetation (data for each tile at every grid point)
16 c ozone (varies with lat, height and time - read it all in
17 c now and interpolate between values later)
18 c the other fixed parameters are:
19 c N2O, Methane (vary with space)
20 c CO2, CFC11, CFC12, CFC22 (set to a global value)
21 c 3) For now, the fizhi package contains the alarms and clocks
22 c routines, so this routine will also initialize the alarms.
23 c-----------------------------------------------------------------------
24 implicit none
25 #include "SIZE.h"
26 #include "fizhi_SIZE.h"
27 #include "fizhi_land_SIZE.h"
28 #include "EEPARAMS.h"
29 #include "fizhi_chemistry_coms.h"
30 #include "fizhi_earth_coms.h"
31 #include "fizhi_land_coms.h"
32 #include "fizhi_ocean_coms.h"
33 #include "chronos.h"
34 #include "gridalt_mapping.h"
35 #include "GRID.h"
36 #include "PARAMS.h"
37
38 integer myThid
39
40 integer i,j,L,bi,bj
41 integer im1, im2, jm1, jm2, idim2, jdim2
42 integer nymdb,nhmsb
43 character*40 vegdata
44 _RL pressure0(Nrphys+1)
45 _RL pressure(Nrphys)
46 _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
47 _RL fracland(sNx,sNy,Nsx,Nsy)
48
49 idim2 = sNx+OLx
50 jdim2 = sNy+OLy
51 im1 = 1
52 im2 = sNx
53 jm1 = 1
54 jm2 = sNy
55 nymdb = nymd0
56 nhmsb = nhms0
57
58 #ifdef ALLOW_MNC
59 if (useMNC) then
60 call fizhi_mnc_init(myThid)
61 endif
62 #endif
63
64 call fizhi_alarms(nymdb,nhmsb,deltaTClock)
65
66 do bj = myByLo(myThid), myByHi(myThid)
67 do bi = myBxLo(myThid), myBxHi(myThid)
68 do j = jm1,jm2
69 do i = im1,im2
70 lons(i,j,bi,bj) = xC(i,j,bi,bj)
71 lats(i,j,bi,bj) = yC(i,j,bi,bj)
72 enddo
73 enddo
74 enddo
75 enddo
76 vegdata = 'veg19232.data'
77 call fizhi_init_veg ( mythid, vegdata,im2,jm2,Nsx,Nsy,
78 . nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,
79 . surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
80
81 C And now fill the earth export landtype
82 do bj = myByLo(myThid), myByHi(myThid)
83 do bi = myBxLo(myThid), myBxHi(myThid)
84 call get_landfrac(im2,jm2,Nsx,Nsy,bi,bj,maxtyp,
85 . surftype,tilefrac,fracland)
86 do j=jm1,jm2
87 do i=im1,im2
88 landtype(i,j,bi,bj) = surftype(i,j,1,bi,bj)
89 if(fracland(i,j,bi,bj).ge.0.3.and.surftype(i,j,1,bi,bj).ge.100)
90 . landtype(i,j,bi,bj) = surftype(i,j,2,bi,bj)
91 if(sice(i,j,bi,bj).ne.0.0)landtype(i,j,bi,bj) = 101
92 enddo
93 enddo
94 enddo
95 enddo
96
97 C Compute pressure profile to get methane and n2o values
98 C - First bottom-up (in mb)
99
100 pressure0(1)=1000.
101 do L = 2,Nrphys+1
102 pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
103 enddo
104 C Now flip pressure to a top-down array and average to mid levels
105 C to send to fizhi init chem
106 do L = 1,Nrphys
107 pressure(L)=(pressure0(Nrphys+2-L)+pressure0(Nrphys+1-L))/2.
108 enddo
109
110 call fizhi_init_chem(mythid,
111 . nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone,
112 . nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,
113 . Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
114
115 do bj = myByLo(myThid), myByHi(myThid)
116 do bi = myBxLo(myThid), myBxHi(myThid)
117 do j=jm1,jm2
118 do i=im1,im2
119 phis_var(i,j,bi,bj) = 0.
120 enddo
121 enddo
122 enddo
123 enddo
124
125 return
126 end

  ViewVC Help
Powered by ViewVC 1.1.22