[MITgcm-support] External temperature forcing

[Jonathan Whitefield]

Hi all.

I was hoping to try to set up a temperature forcing of a specific grid
cell, much like the freshwater input forces a cell to include a
certain river runoff. I noticed in the exf_mapfields.F some lines that
alter the heat flux, but I want to change the temperature of the cell
BEFORE fluxes are applied. So, my questions are:

1. Which file(s) would I change to force a temperature?
2. Assuming that my river runoff temperatures are on a 360 x 180 x 12
(i.e. monthly mean) grid, could I just use something like:
    DO j = 1, sNy
      DO i = 1, sNx
          theta(i,j,1,bi,bj) = runoftemp(i,j,bi,bj)
      ENDDO
    ENDDO

If not, could someone please suggest the syntax that I would use? The
above lines are based on exf_mapfields.F (in the #ifdef ALLOW_ATM_TEMP
section), and assuming that where it originally says
theta(i,j,ks,bi,bj), ks is the model grid level.

Thanks!

--
Jonathan Whitefield, M.Sc.
Ph.D. Student, Physical Oceanography

110 O'Neill Building,
P.O. Box 757220
University of Alaska Fairbanks
Fairbanks, AK 99775-7220

Phone: (907) 474-5184
Fax: (907) 474-5863