(u,v)_visc not writing to diagnostics file

[ggerbi]

Using revision 467 I found that the diagnostics were not writing the vertical velocity terms. It looks like this is due to missing code in inp_par.F.

I have added lines at 2767 that seem to have fixed this.

Code: Select all

          ELSE IF (TRIM(KeyWord).eq.'Dout(M3vvis)') THEN
            IF (M3vvis.le.0) THEN
              IF (Master) WRITE (out,280) 'M3vvis'
              exit_flag=5
              RETURN
            END IF
            Npts=load_l(Nval, Cval, Ngrids, Ldiag)
            DO ng=1,Ngrids
              Dout(idDu3d(M3vvis),ng)=Ldiag(ng)
              Dout(idDv3d(M3vvis),ng)=Ldiag(ng)
            END DO

Also, I have not looked at this in detail, but inp_par.F may also be missing a couple lines for M3yvis around 4500 when compared to lines for M3xvis.

[arango]

Yes, the processing of vertical mixing for momentum diagnostics and reporting Dout(M3yvis) were missing. Thank you for reporting this problem. I corrected the repository.