[MITgcm-devel] ptracer adjoint

[Matthew Mazloff]

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