[MITgcm-support] verification_other case global_oce_cs32 fails at runtime in adjoint mode

[Martin Losch]

Hi Dan,

I tried this on my linux box without MPI and I cannot reproduce your problem (I used MITgcm/verification_other.git and not the CVS MITgcm_contrib/verification_other, which appears to be out of date). I grepped the code for “m_boxmean_theta” and only found this:

(base) bkli04l006::build (master)> grep m_boxmean_theta *.f
ad_input_code_ad.f:     $'m_boxmean_theta') then
ad_input_code_ad.f:     $'m_boxmean_theta') then
ad_input_code_ad.f:     $'m_boxmean_theta') then
ad_input_code_ad.f:     $'m_boxmean_theta') then
ad_input_code.f:            if (gencost_barfile(kgen)(1:15).EQ.'m_boxmean_theta') then
ad_taf_output.f:     $'m_boxmean_theta') then
ad_taf_output.f:     $'m_boxmean_theta') then
ad_taf_output.f:     $'m_boxmean_theta') then
ad_taf_output.f:     $'m_boxmean_theta') then
ecco_check.f:     &          (gencost_barfile(k)(1:15).EQ.'m_boxmean_theta').OR.
ecco_phys.f:            if (gencost_barfile(kgen)(1:15).EQ.'m_boxmean_theta') then

(and I made sure that there’s this is really just m_boxmean_theta). Where in your code (which routine) does the model try to read adm_boxmean_theta?

Martin
> On 21. Sep 2020, at 11:00, Dan Jones <dcjones.work at gmail.com> wrote:
> 
> Hello. 
> 
> Apologies for the cross-posting - I've posted this as a GitHub issue, but I thought I should put it here as well.
> 
> I am trying to build and test the global_oce_cs32 verification_other exercise using the code in the input_ad.sens directory. The forward case compiles and runs without error. The adjoint case (built using TAF) compiles without error, but at runtime I receive the following error in STDOUT:
> 
> (PID.TID 0000.0001)  MDS_READVEC_LOC: open file: south30_maskT 
> (PID.TID 0000.0001)  MDS_RD_REC_RL: iRec,Dim =         9          1
> (PID.TID 0000.0001)  MDS_READ_FIELD: filename: adm_boxmean_theta.0000000000.data
> (PID.TID 0000.0001)  MDS_READ_FIELD: File does not exist
> 
> and this error in STDERR:
> 
> (PID.TID 0000.0001) *** ERROR ***  MDS_READ_FIELD: filename: adm_boxmean_theta.0000000000.data
> (PID.TID 0000.0001) *** ERROR ***  MDS_READ_FIELD: File does not exist
> 
> My MITgcm source code is up-to-date with the master. I am running on archer.ac.uk in parallel mode using 24 cores.
> 
> What should I try here? I haven't run into this error before using other adjoint setups, at least not that I can recall. Should I just create an empty "dummy" file to start with? Thanks in advance for any help/guidance.
> 
> Best regards,
> Dan
> 
> --------------------------------------------------------------
> Dr Dan Jones / British Antarctic Survey
> danjonesocean.com / @DanJonesOcean
> --------------------------------------------------------------
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support