[MITgcm-devel] depressed Eta and sIceLoad

[Martin Losch]

What did you actually do? If you just run the ECCO model after  
optimization (that is when the surface fluxes are adjusted so that  
the near surface salinities are not too bad) and suddenly your EmPmR  
changes dramatically (with evap*(1-area) the evap contribution under  
ice is basically zero), I would be surprised if the solution did not  
deteriorate.
If, however, you ran the first guess with the unoptimized surface  
fluxes (which I believe you did), then I am surprised that you have  
an effect that is so much different from mine.

Martin

PS. Just to make sure, the line in question should look like:
> C NOW GET FRESH WATER FLUX
>           EmPmR(I,J,bi,bj)= maskC(I,J,kSurface,bi,bj)*(
>      &         EVAP(I,J,bi,bj)*(ONE-AREA(I,J,2,bi,bj))
>      &         -RUNOFF(I,J,bi,bj)
>      &         +SEAICE_SALT(I,J,bi,bj)*0.92 _d 0/SEAICE_deltaTtherm
>      &         )

On 26 Oct 2006, at 16:58, Patrick Heimbach wrote:

>
> Hi there,
>
> could you hold off for another day with
> putting the EVAP fix.
> I ran a quick test with ECCO using the
> EVAP fix, and things looked really bad there.
> But can't give anything conclusive yet.
>
> -p.
>
>
>
> On Thu, 2006-10-26 at 10:24, Dimitris Menemenlis wrote:
>> Martin, my vote is to add all the fixes, as you describe in your e- 
>> mail, without
>> caring for backward compatibility.  Although I have not yet tried  
>> it, I agree
>> that the evap is potentially the most serious bug, especially for  
>> those
>> integrations that include Southern Ocean.  Regarding too much snow  
>> accumulation,
>> if the problem is with the NCEP or CORE forcing fields, then we  
>> should be fixing
>> the NCEP or CORE forcing fields, not the model.
>>
>> The link:
>> http://ecco2.jpl.nasa.gov/data2/cube/cube38/pickup/HSNOW.jpg
>> is live once again, it was down when you last tried to access it.
>> This particular cubed-sphere integration is driven by NCEP  
>> reanalysis.
>> As Jinlun pointed out the "too-much snow" problem is mostly around  
>> Antarctica.
>> In the Arctic it is less of a problem, although there are some  
>> small regions
>> with unrealistic 25 m snow around Greenland and in the Canadian  
>> Archipelago.
>>
>> The companion figures to above, that show complete contents of  
>> seaice_pickup
>> from same integration and time in Arctic and Antarctic are
>> http://ecco2.jpl.nasa.gov/data2/cube/cube38/pickup/ 
>> pickup_seaiceArctic.jpg
>> http://ecco2.jpl.nasa.gov/data2/cube/cube38/pickup/ 
>> pickup_seaiceAntarctic.jpg
>> A yet different problem is ice that becomes too thick, up to 25 m  
>> in some
>> isolated spots.  Shall we use this same occasion to uncomment
>>
>> #ifdef DO_WE_NEED_THIS
>> c          HEFF(I,J,1,bi,bj)=MIN(MAX_HEFF,HEFF(I,J,1,bi,bj))
>> #endif
>>
>> in growth.F.  We can set a large default MAX_HEFF in  
>> seaice_readparms.F
>>
>> Dimitris
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
> -- 
> Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
> MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
> FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel