[MITgcm-support] how to initialize seaice package

[Matthew Mazloff]

OK...sorry for all the emails....problem #2 solved.

I see that permanent checkpoints are only written for seaice (and  
obcs package etc.) at the last time step.  until then the file is  
written into pickup_seaice.ckptA.data.  So all my problems are  
solved...for now.

thanks!
matt


On Sep 13, 2006, at 3:00 PM, Matthew Mazloff wrote:

> 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
>
> _______________________________________________
> 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/decfbd39/attachment.htm>