[MITgcm-devel] thsice_get_exf

[Martin Losch]

Hi Jean-Michel,

could you check for me the following parts of thsice_get_exf:
We have:
         IF ( hSnow(i,j).GT.3. _d -1 ) THEN
           iceornot=2
          ELSE
	           iceornot=1
          ENDIF
So that iceornot == 0 never happens. But then we have
                  IF ( iceornot.EQ.0 ) THEN
	           lath  = flamb
	           qsat_fac = cvapor_fac
	           qsat_exp = cvapor_exp
	          ELSE
	           lath  = flamb+flami
	           qsat_fac = cvapor_fac_ice
	           qsat_exp = cvapor_exp_ice
	          ENDIF
which will always be the second case, right? Similarly:
C--- Upward long wave radiation
	          IF ( iceornot.EQ.0 ) THEN
	            emiss = ocean_emissivity
	          ELSEIF (iceornot.EQ.2) THEN
	            emiss = snow_emissivity
	          ELSE
	            emiss = ice_emissivity
	          ENDIF
Where emiss will always be snow or ice_emissivity, but never ocean, correct?

Can we get rid of the iceornot==0 option (otherwise I'll have to introduce more 2d fields for lath and emiss, which I want to avoid)

Martin