[MITgcm-support] External temperature forcing
Menemenlis, Dimitris (3248)
Dimitris.Menemenlis at jpl.nasa.gov
Mon Jul 1 19:17:40 EDT 2013
Jonathan, you can use pkg/rbcs to do what you describe below in arbitrary locations
http://mitgcm.org/public/r2_manual/latest/online_documents/node237.html
But if you are specifically going to force temperature only where there is runoff,
then use RUNOFTEMP as per example in verification/global_ocean.cs32x15/input.seaice
Dimitris Menemenlis
On Jul 1, 2013, at 3:01 PM, Jonathan Whitefield wrote:
> 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
More information about the MITgcm-support
mailing list