[MITgcm-support] how to initialize seaice package

[Matthew Mazloff]

Hi again,

1 problem solved...I gave the HeffFile as single precision....oops.

now i just need a pickup_seaice file

-Matt

On Sep 13, 2006, at 2:53 PM, Matthew Mazloff wrote:

> Hi Dimitris,
>
> When I set       SEAICE_initialHEFF = 0.0  AREA is initialized to  
> 0.0, and when I set SEAICE_initialHEFF = 1.0  AREA is initialized  
> to 1.0.   This is good.
>
> However, this is regardless of the fact  that, as STDOUT.0000 shows
> (PID.TID 0000.0001) >      HeffFile  = 'SO6input/SO6_INI_ICECONC.bin',
> so this is bad; Heff should be read in...and AREA should have a  
> pattern that reflects this.  I do not get this...I am not sure  
> why...I will have to work through this issue....it is very  
> odd....luckily it appears to be easy to debug...not too much code
>
> Another bizarre issue with the seaice package.  I am not getting a  
> pickup_seaice.*.data file.  Again the code seems straightforward  
> and yet....
>
> I'll try to work these issues out.  any ideas?
>
> Thanks,
> Matt
>
>
>
>
> On Sep 13, 2006, at 2:22 PM, Dimitris Menemenlis wrote:
>
>> Hi Matt, actually HeffFile is supposed to represent initial  
>> effective ice thickness not area.  Look at seaice_init.F for  
>> initialization logic.
>>
>> First there is
>>
>>             HEFF(I,J,k,bi,bj)=SEAICE_initialHEFF*HEFFM(i,j,bi,bj)
>>
>> where HEFFM is a land mask.  Then there is
>>
>>        IF ( HeffFile .NE. ' ' ) THEN
>>         CALL READ_FLD_XY_RL( HeffFile, ' ', ZETA, 0, myThid )
>>              HEFF(I,J,k,bi,bj) = MAX(ZETA(i,j,bi,bj),ZERO)
>>
>> Finally there is
>>
>>             IF(HEFF(I,J,k,bi,bj).GT.ZERO) AREA(I,J,k,bi,bj)=ONE
>>
>> This has worked for me in the past so I am not sure why it is not  
>> working for you.  How does your initial sea ice area compare with  
>> SO6_INI_ICECONC.bin ?
>>
>> Dimitris
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20060913/36b91678/attachment.htm>