[MITgcm-devel] another bug in growth.F ?

[Dimitris Menemenlis]

Jinlun, the beer/crap comment was in jest.  Everyone who has used pkg/seaice 
appreciates your effort in making this package available to MITgcm and also your 
subsequent help with bug fixes and with other modifications.

Martin, I also find that

> cdm       IF(FICE(I,J,bi,bj).GT.ZERO) THEN
>           IF(atemp(i,j,bi,bj).LE.273.15 _d 0 ) THEN

has very little impact on growth.F both for the forward solution as well as for 
the high forward sensitivity of the model that you and Patrick reported.  What 
does remove the high forward sensitivity is commenting out the snow-melt addition.

> C Now melt snow if there is residual heat left in surface level                 
> C Note that units of YNEG and SEAICE_SALT are m of ice                          
> cdm       IF(RESID_HEAT(I,J,bi,bj).GT.ZERO.AND.                                 
> cdm  &         HSNOW(I,J,bi,bj).GT.ZERO) THEN                                   
> cdm        GHEFF(I,J)=MIN(HSNOW(I,J,bi,bj)/SDF/ICE_DENS,                        
> cdm  &         RESID_HEAT(I,J,bi,bj))                                           
> cdm        YNEG(I,J,bi,bj)=YNEG(I,J,bi,bj)+GHEFF(I,J)                           
> cdm        HSNOW(I,J,bi,bj)=HSNOW(I,J,bi,bj)-GHEFF(I,J)*SDF*ICE_DENS            
> cdm        SEAICE_SALT(I,J,bi,bj)=SEAICE_SALT(I,J,bi,bj)-GHEFF(I,J)             
> cdm       ENDIF                                                                 

So back to where we were before latest exchange.

Dimitris

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720