
in release code:  .../models/lnd/clm/tools directory


modfications to mksurfdata code for paleo:


1) mksoicol.F90 (after line 328)

!CAS setting 15 everying in input file
!comment out nsoicol max definition as dependent on
!input file and set explicitly to 20.
  
 ! nsoicol = maxval(soil_color_i)
  nsoicol =  20   
!CAS end mod
  


2) mkvocef.F90
set isoprene emission to zero: mkvocef.F90

example:
!CAS looks like this is hardwired for .5 degree
!CAS default paleo should have all of these values set to zero anyway
  do in = 1,720
     do jn = 1, 360
        if (in < 361) then 
           !ef_btr_i(in,361-jn)=temp_i(in+360,jn)*1.0_r8
        end if
        if (in > 360) then 
           !ef_btr_i(in,361-jn)=temp_i(in-360,jn)*1.0_r8
        end if    
     end do
  end do

notes:
	continue for all ef_* variables (several do loops)
	after half degree stuff has been commented, explicitly
	set to zero.
!CAS 
  ef_btr_i = 0.      ! input grid: EFs for broadleaf trees
  ef_fet_i = 0.      ! input grid: EFs for broadleaf trees
  ef_fdt_i = 0.      ! input grid: EFs for broadleaf trees
  ef_shr_i = 0.      ! input grid: EFs for broadleaf trees
  ef_grs_i = 0.      ! input grid: EFs for broadleaf trees
  ef_crp_i = 0.      ! input grid: EFs for broadleaf trees
!CAS



3) mksrfdat.F90
remove hardwired code for ross ice shelf and don't use dynamic land stuff:   mksrfdat.F90

   a) don't use dynamic land stuff, insert before line 617:

code snippet:

!CAS skip "dynamic land", i.e. for transient runs...not PALEO 
    mksrf_fdynuse = ' '

    if (mksrf_fdynuse /= ' ') then

end code snippet

    b) comment out ross ice shelf hardwire:
          ! Set land values on Ross ice shelf to glacier
!CAS remove for paleo
!
!         if (ldomain%latixy(i,j) < rosslat) then
!
!            pctpft(i,j,:)  = 0._r8
!
!            pctlak(i,j)      = 0._r8
!            pctwet(i,j)      = 0._r8
!            pcturb(i,j)      = 0._r8
!
!            pctgla(i,j)      = 100._r8
!
!         end if
!CAS end mod



