New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14534 – NEMO

Changeset 14534


Ignore:
Timestamp:
2021-02-23T12:22:02+01:00 (3 years ago)
Author:
agn
Message:

added extra debugging diagnostics

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/NERC/dev_r11078_OSMOSIS_IMMERSE_Nurser_4.0/src/OCE/ZDF/zdfosm.F90

    r14521 r14534  
    10551055       WRITE(narea+100,'(a,*(g11.3))') ' ghamu[imld-1..ibld+2] =', ( ghamu(ji,jj,jk), jk=jl,jm ) 
    10561056       WRITE(narea+100,'(a,*(g11.3))') ' ghamv[imld-1..ibld+2] =', ( ghamv(ji,jj,jk), jk=jl,jm ) 
     1057       WRITE(narea+100,*) 
    10571058       FLUSH(narea+100) 
    10581059    END IF 
    … …  
    12031204       WRITE(narea+100,'(a,*(g11.3))') ' ghamu[imld-1..ibld+2] =', ( ghamu(ji,jj,jk), jk=jl,jm ) 
    12041205       WRITE(narea+100,'(a,*(g11.3))') ' ghamv[imld-1..ibld+2] =', ( ghamv(ji,jj,jk), jk=jl,jm ) 
     1206       WRITE(narea+100,*) 
    12051207       FLUSH(narea+100) 
    12061208    END IF 
    … …  
    13161318       END DO 
    13171319    END DO 
     1320#ifdef key_osm_debug 
     1321    IF(narea==nn_narea_db) THEN 
     1322       ji=iloc_db; jj=jloc_db 
     1323       jl = imld(ji,jj) - 1; jm = MIN(ibld(ji,jj) + 2, mbkt(ji,jj) ) 
     1324       WRITE(narea+100,'(2(a,g11.3))')'Stokes + buoy + pyc + transport contribs to ghamt/contrib to ghamt/s:  zsc_wth_1=',zsc_wth_1(ji,jj), '  zsc_ws_1=',zsc_ws_1(ji,jj) 
     1325       IF (lpyc(ji,jj)) WRITE(narea+100,'(2(a,g11.3))') 'zsc_wth_pyc=', zsc_wth_pyc(ji,jj), '  zsc_wth_pyc=',zsc_wth_pyc(ji,jj) 
     1326       WRITE(narea+100,'(a,*(g11.3))') ' ghamt[imld-1..ibld+2] =', ( ghamt(ji,jj,jk), jk=jl,jm ) 
     1327       WRITE(narea+100,'(a,*(g11.3))') ' ghams[imld-1..ibld+2] =', ( ghams(ji,jj,jk), jk=jl,jm ) 
     1328       IF( lconv(ji,jj) ) THEN 
     1329          WRITE(narea+100,'(2(a,g11.3))')'Unstable; transport contrib to ghamu/v:  zsc_uw_1=',zsc_uw_1(ji,jj), '  zsc_vw_1=',zsc_vw_1(ji,jj) 
     1330       ELSE 
     1331          WRITE(narea+100,'(3(a,g11.3))')'Stable; transport contrib to ghamu/v:  zsc_uw_1=',zsc_uw_1(ji,jj), '  zsc_vw_1=',zsc_vw_1(ji,jj), & 
     1332               &' zsc_uw_2=',zsc_uw_2(ji,jj) 
     1333       END IF 
     1334       WRITE(narea+100,'(a,*(g11.3))') ' ghamu[imld-1..ibld+2] =', ( ghamu(ji,jj,jk), jk=jl,jm ) 
     1335       WRITE(narea+100,*) 
     1336       FLUSH(narea+100) 
     1337    END IF 
     1338#endif 
    13181339 
    13191340    IF(ln_dia_osm) THEN 
    … …  
    13761397       END DO       ! ji loop 
    13771398    END DO          ! jj loop 
     1399#ifdef key_osm_debug 
     1400    IF(narea==nn_narea_db) THEN 
     1401       ji=iloc_db; jj=jloc_db 
     1402       jl = imld(ji,jj) - 1; jm = MIN(ibld(ji,jj) + 2, mbkt(ji,jj) ) 
     1403       WRITE(narea+100,'(a)')'Tweak gham[uv] to go to zero near surface, add pycnocline viscosity/diffusivity  & set=0 at ibld' 
     1404       WRITE(narea+100,'(a,*(g11.3))') ' ghamt[imld-1..ibld+2] =', ( ghamt(ji,jj,jk), jk=jl,jm ) 
     1405       WRITE(narea+100,'(a,*(g11.3))') ' ghams[imld-1..ibld+2] =', ( ghams(ji,jj,jk), jk=jl,jm ) 
     1406       WRITE(narea+100,'(a,*(g11.3))') ' ghamu[imld-1..ibld+2] =', ( ghamu(ji,jj,jk), jk=jl,jm ) 
     1407       WRITE(narea+100,'(a,*(g11.3))') ' ghamv[imld-1..ibld+2] =', ( ghamv(ji,jj,jk), jk=jl,jm ) 
     1408       WRITE(narea+100,*) 
     1409       FLUSH(narea+100) 
     1410    END IF 
     1411#endif 
    13781412 
    13791413    IF(ln_dia_osm) THEN 
    … …  
    14621496       END DO 
    14631497    END IF ! ln_convmix = .true. 
     1498#ifdef key_osm_debug 
     1499    IF(narea==nn_narea_db) THEN 
     1500       ji=iloc_db; jj=jloc_db 
     1501       jl = imld(ji,jj) - 1; jm = MIN(ibld(ji,jj) + 2, mbkt(ji,jj) ) 
     1502       WRITE(narea+100,'(a)') ' After including KPP Ri# diffusivity & viscosity' 
     1503       WRITE(narea+100,'(a,*(g11.3))') ' zdiffut[imld-1..ibld+2] =', ( zdiffut(ji,jj,jk), jk=jl,jm ) 
     1504       WRITE(narea+100,'(a,*(g11.3))') ' zviscos[imld-1..ibld+2] =', ( zviscos(ji,jj,jk), jk=jl,jm ) 
     1505       WRITE(narea+100,*) 
     1506       FLUSH(narea+100) 
     1507    END IF 
     1508#endif 
    14641509 
    14651510 
    … …  
    14951540          END DO 
    14961541       END DO 
     1542#ifdef key_osm_debug 
     1543       IF(narea==nn_narea_db) THEN 
     1544          ji=iloc_db; jj=jloc_db 
     1545          jl = imld(ji,jj) - 1; jm = MIN(ibld(ji,jj) + 2, mbkt(ji,jj) ) 
     1546          WRITE(narea+100,'(a)') ' After including FK diffusivity & non-local terms' 
     1547          WRITE(narea+100,'(a,*(g11.3))') ' zdiffut[imld-1..ibld+2] =', ( zdiffut(ji,jj,jk), jk=jl,jm ) 
     1548          WRITE(narea+100,'(a,*(g11.3))') ' ghamt[imld-1..ibld+2] =', ( ghamt(ji,jj,jk), jk=jl,jm ) 
     1549          WRITE(narea+100,'(a,*(g11.3))') ' ghams[imld-1..ibld+2] =', ( ghams(ji,jj,jk), jk=jl,jm ) 
     1550          WRITE(narea+100,*) 
     1551          FLUSH(narea+100) 
     1552       END IF 
     1553#endif 
    14971554    ENDIF 
    14981555 
    … …  
    15401597    CALL lbc_lnk_multi( 'zdfosm', ghamt, 'W', 1. , ghams, 'W', 1.,   & 
    15411598         &                  ghamu, 'U', -1. , ghamv, 'V', -1. ) 
     1599#ifdef key_osm_debug 
     1600    IF(narea==nn_narea_db) THEN 
     1601       ji=iloc_db; jj=jloc_db 
     1602       jl = imld(ji,jj) - 1; jm = MIN(ibld(ji,jj) + 2, mbkt(ji,jj) ) 
     1603       WRITE(narea+100,'(a)') ' Final diffusivity & viscosity, & non-local terms' 
     1604       WRITE(narea+100,'(a,*(g11.3))') ' p_avt[imld-1..ibld+2] =', ( p_avt(ji,jj,jk), jk=jl,jm ) 
     1605       WRITE(narea+100,'(a,*(g11.3))') ' p_avm[imld-1..ibld+2] =', ( p_avm(ji,jj,jk), jk=jl,jm ) 
     1606       WRITE(narea+100,'(a,*(g11.3))') ' ghamt[imld-1..ibld+2] =', ( ghamt(ji,jj,jk), jk=jl,jm ) 
     1607       WRITE(narea+100,'(a,*(g11.3))') ' ghams[imld-1..ibld+2] =', ( ghams(ji,jj,jk), jk=jl,jm ) 
     1608       WRITE(narea+100,'(a,*(g11.3))') ' ghamu[imld-1..ibld+2] =', ( ghamu(ji,jj,jk), jk=jl,jm ) 
     1609       WRITE(narea+100,'(a,*(g11.3))') ' ghamv[imld-1..ibld+2] =', ( ghamv(ji,jj,jk), jk=jl,jm ) 
     1610       WRITE(narea+100,*) 
     1611       FLUSH(narea+100) 
     1612    END IF 
     1613#endif 
    15421614 
    15431615    IF(ln_dia_osm) THEN 
Note: See TracChangeset for help on using the changeset viewer.