--- MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F	2004/08/20 13:47:36	1.2
+++ MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F	2005/05/24 14:57:00	1.3
@@ -1,6 +1,7 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F,v 1.2 2004/08/20 13:47:36 molod dead $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F,v 1.3 2005/05/24 14:57:00 molod Exp $
 C $Name:  $
 
+#include "FIZHI_OPTIONS.h"
        subroutine fizhi_init_fixed (myThid)
 c-----------------------------------------------------------------------
 c  Routine to initialise the fizhi package.
@@ -17,35 +18,107 @@
 c      the other fixed parameters are:
 c         N2O, Methane (vary with space)
 c         CO2, CFC11, CFC12, CFC22 (set to a global value)
-c   2) This routine will also read the data.fizhi namelist interface
-c      to obtain information about diagnostic selection.
+c   3) For now, the fizhi package contains the alarms and clocks
+c      routines, so this routine will also initialize the alarms.
 c-----------------------------------------------------------------------
        implicit none
-#include "CPP_OPTIONS.h"
 #include "SIZE.h"
 #include "fizhi_SIZE.h"
-#include "land_SIZE.h"
+#include "fizhi_land_SIZE.h"
 #include "EEPARAMS.h"
-#include "chemistry_coms.h"
-#include "earth_coms.h"
+#include "fizhi_chemistry_coms.h"
+#include "fizhi_earth_coms.h"
+#include "fizhi_land_coms.h"
+#include "fizhi_ocean_coms.h"
+#include "chronos.h"
+#include "gridalt_mapping.h"
+#include "GRID.h"
+#include "PARAMS.h"
 
        integer myThid
 
-       integer bi, bj
-       integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
+       integer i,j,L,bi,bj
+       integer im1, im2, jm1, jm2, idim2, jdim2
+       integer nymdb,nhmsb
+       character*40 vegdata
+       _RL pressure0(Nrphys+1)
+       _RL pressure(Nrphys)
+       _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
+       _RL fracland(sNx,sNy,Nsx,Nsy)
+
+       idim2 = sNx+OLx
+       jdim2 = sNy+OLy
+       im1 = 1
+       im2 = sNx
+       jm1 = 1
+       jm2 = sNy
+       nymdb = nymd0
+       nhmsb = nhms0
+
+#ifdef ALLOW_MNC
+       if (useMNC) then
+        call fizhi_mnc_init(myThid)
+       endif
+#endif
 
-       im1 = 1-OLx
-       im2 = sNx+OLx
-       jm1 = 1-OLy
-       jm2 = sNy+OLy
-       idim1 = 1
-       idim2 = sNx
-       jdim1 = 1
-       jdim2 = sNy
+       call fizhi_alarms(nymdb,nhmsb,deltaTClock)
 
        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
-C Dummy routine - nothing to do for now
+       do j = jm1,jm2
+       do i = im1,im2
+        lons(i,j,bi,bj) = xC(i,j,bi,bj)
+        lats(i,j,bi,bj) = yC(i,j,bi,bj)
+       enddo
+       enddo
+       enddo
+       enddo
+       vegdata = 'veg19232.data'
+       call fizhi_init_veg ( mythid, vegdata,im2,jm2,Nsx,Nsy,
+     .         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,
+     .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
+
+C And now fill the earth export landtype
+       do bj = myByLo(myThid), myByHi(myThid)
+       do bi = myBxLo(myThid), myBxHi(myThid)
+        call get_landfrac(im2,jm2,Nsx,Nsy,bi,bj,maxtyp,
+     .        surftype,tilefrac,fracland)
+        do j=jm1,jm2
+        do i=im1,im2
+         landtype(i,j,bi,bj) = surftype(i,j,1,bi,bj)
+         if(fracland(i,j,bi,bj).ge.0.3.and.surftype(i,j,1,bi,bj).ge.100)
+     .            landtype(i,j,bi,bj) = surftype(i,j,2,bi,bj)
+         if(sice(i,j,bi,bj).ne.0.0)landtype(i,j,bi,bj) = 101
+        enddo
+        enddo
+       enddo
+       enddo
+
+C Compute pressure profile to get methane and n2o values 
+C     - First bottom-up (in mb)
+
+       pressure0(1)=1000.
+       do L = 2,Nrphys+1
+        pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
+       enddo
+C Now flip pressure to a top-down array and average to mid levels
+C     to send to fizhi init chem
+       do L = 1,Nrphys
+        pressure(L)=(pressure0(Nrphys+2-L)+pressure0(Nrphys+1-L))/2.
+       enddo
+
+       call fizhi_init_chem(mythid,
+     .      nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone,
+     .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,
+     .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
+
+       do bj = myByLo(myThid), myByHi(myThid)
+       do bi = myBxLo(myThid), myBxHi(myThid)
+        do j=jm1,jm2
+        do i=im1,im2
+         phis_var(i,j,bi,bj) = 0.
+        enddo
+        enddo
        enddo
        enddo