[MITgcm-devel] ptracer adjoint

[Jean-Michel Campin]

Hi Matt,

I wonder if this is not related to the default OB values
(i.e., in case you don't provide an OB ptracer value,
it takes the value next to the OB), since, in a FWD run, 
Oliver experienced some exponential growth of tracer 
concentration next to the OB.

Are you in this case (default OB values) ?

Cheers,
Jean-Michel

On Mon, Dec 13, 2010 at 08:07:51PM -0800, Matthew Mazloff wrote:
> Hello,
>
> When the sensitivity to passive tracer concentration hits an open  
> boundary the solution blows up -- anyone seen this before?  Any ideas?
>
> Thanks
> Matt
>
>
>
> On Dec 9, 2010, at 2:00 PM, Matthew Mazloff wrote:
>
>> Hello,
>>
>> I am running a passive tracer adjoint -- all is fine but 2 things I
>> wanted to share
>>
>> 1)  to avoid recomp had to do
>> login4% diff obcs_ad_check_lev1_dir.h ../../../pkg/obcs/.
>> 67,70d66
>> < #ifdef ALLOW_OBCS_PRESCRIBE
>> < CADJ STORE OBNptr0 = comlev1, key = ikey_dynamics
>> < CADJ STORE OBNptr1 = comlev1, key = ikey_dynamics
>> < #endif /* ALLOW_OBCS_PRESCRIBE */
>> 74,77d69
>> < #ifdef ALLOW_OBCS_PRESCRIBE
>> < CADJ STORE OBSptr0 = comlev1, key = ikey_dynamics
>> < CADJ STORE OBSptr1 = comlev1, key = ikey_dynamics
>> < #endif /* ALLOW_OBCS_PRESCRIBE */
>> 81,84d72
>> < #ifdef ALLOW_OBCS_PRESCRIBE
>> < CADJ STORE OBEptr0 = comlev1, key = ikey_dynamics
>> < CADJ STORE OBEptr1 = comlev1, key = ikey_dynamics
>> < #endif /* ALLOW_OBCS_PRESCRIBE */
>> 88,91d75
>> < #ifdef ALLOW_OBCS_PRESCRIBE
>> < CADJ STORE OBWptr0 = comlev1, key = ikey_dynamics
>> < CADJ STORE OBWptr1 = comlev1, key = ikey_dynamics
>> < #endif /* ALLOW_OBCS_PRESCRIBE */
>>
>>
>> 2) ADJptracer is not written by default.  I attempted to add it to
>> adcommon.h and addummy_in_stepping.F, but for whatever reason the
>> output written was all zeros.  So I just hardcoded
>>       IF ( MOD((myIter-1),432).eq.0) THEN
>>         WRITE(suff,'(I10.10)') myIter
>>        CALL WRITE_FLD_XYZ_RL(
>>      &       'ADJptracer.',suff, adptracer, myIter, myThid)
>>       ENDIF
>> to subroutine adthermodynamics from ad_taf_output.f
>> and all seems ok...
>>
>>
>> Thought sharing this experience may help....
>> -Matt
>>
>>
>>
>>
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel