[Mitgcm-support] Re: pkg/seaice and pkg/exf

[mitgcm-support at dev.mitgcm.org]

Tom, would you be able to plot FU, FV, EmPmR, Qsw, and Qnet, and see if
there are any differences in Gulf of Mexico, with useSEAICE=.FALSE. and
with useSEAICE=.TRUE..

With SEAICE_EXTERNAL_FLUXES defined, the surface fluxes should be
identical everywhere where there is no sea-ice.

The surface fluxes are dumped automatically by seaice_do_diags.F
when useSEAICE=.TRUE. and SEAICEwriteState=.TRUE. in data.seaice.

For pkg/exf, you will need to add the following lines right after
completion of EXF_GETFORCING in forward_step.F:

      CHARACTER*(MAX_LEN_MBUF) suff

            WRITE(suff,'(I10.10)') myIter
            _BARRIER
            _BEGIN_MASTER( myThid )
            CALL WRITE_FLD_XY_RS( 'FU.',suff,fu,myIter,myThid)
            CALL WRITE_FLD_XY_RS( 'FV.',suff,fv,myIter,myThid)
            CALL WRITE_FLD_XY_RS( 'EmPmR.',suff,EmPmR,myIter,myThid)
            CALL WRITE_FLD_XY_RS( 'Qnet.',suff,Qnet,myIter,myThid)
            CALL WRITE_FLD_XY_RS( 'Qsw.',suff,Qsw,myIter,myThid)
            _END_MASTER( myThid )
            _BARRIER

Dimitris

-- 
Dimitris Menemenlis                      menemenlis at jpl.nasa.gov
Jet Propulsion Lab, MS 300-323           tel: 818-354-1656
4800 Oak Grove Dr, Pasadena CA 91109     fax: 818-393-6720