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 11367 for NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/OCE/SBC/sbccpl.F90 – NEMO

Ignore:
Timestamp:
2019-07-29T18:10:24+02:00 (5 years ago)
Author:
dancopsey
Message:

More print statements including fix

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/OCE/SBC/sbccpl.F90

    r11355 r11367  
    13791379      ENDIF 
    13801380 
    1381       IF(narea == 419) THEN 
    1382          WRITE(9419,*) 'k_ice, kt, k_fsbc = ',k_ice, kt, k_fsbc 
    1383          WRITE(9419,*) 'MOD( kt-1, k_fsbc ) = ',MOD( kt-1, k_fsbc ) 
     1381      IF(narea == 68) THEN 
     1382         WRITE(968,*) 'k_ice, kt, k_fsbc = ',k_ice, kt, k_fsbc 
     1383         WRITE(968,*) 'MOD( kt-1, k_fsbc ) = ',MOD( kt-1, k_fsbc ) 
    13841384      ENDIF 
    13851385       
    … …  
    14021402         ENDIF 
    14031403 
    1404          IF(narea == 419) THEN 
    1405             WRITE(9419,*) 'max zemp after coupling = ',MAXVAL(  ABS( zemp(:,:) )  ) 
     1404         IF(narea == 68) THEN 
     1405            WRITE(968,*) 'max zemp after coupling = ',MAXVAL(  ABS( zemp(:,:) )  ) 
    14061406         ENDIF 
    14071407 
    … …  
    14111411         IF( srcv(jpr_cal)%laction )     zemp(:,:) = zemp(:,:) - frcv(jpr_cal)%z3(:,:,1) 
    14121412 
    1413          IF(narea == 419) THEN 
    1414             WRITE(9419,*) 'max zemp after runoff = ',MAXVAL(  ABS( zemp(:,:) )  ) 
    1415             WRITE(9419,*) 'max emp after runoff = ',MAXVAL(  ABS( emp(:,:) )  ) 
     1413         IF(narea == 68) THEN 
     1414            WRITE(968,*) 'max zemp after runoff = ',MAXVAL(  ABS( zemp(:,:) )  ) 
     1415            WRITE(968,*) 'max emp after runoff = ',MAXVAL(  ABS( emp(:,:) )  ) 
    14161416         ENDIF 
    14171417  
    … …  
    14261426         ENDIF 
    14271427 
    1428          IF(narea == 419) THEN 
    1429             WRITE(9419,*) 'max emp after conversion from zemp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     1428         IF(narea == 68) THEN 
     1429            WRITE(968,*) 'max emp after conversion from zemp = ',MAXVAL(  ABS( emp(:,:) )  ) 
    14301430         ENDIF 
    14311431 
    … …  
    15961596          
    15971597      ENDIF 
     1598 
     1599      IF(narea == 68) THEN 
     1600         WRITE(968,*) 'Has OASIS sent data = ',(  nrcvinfo(itx) == OASIS_Rcv ) 
     1601         WRITE(968,*) 'srcv(jpr_itx1)%laction = ',srcv(jpr_itx1)%laction 
     1602         WRITE(968,*) 'frcv(jpr_itx1)%z3(54,1,1) = ',frcv(jpr_itx1)%z3(54,1,1) 
     1603         WRITE(968,*) 'frcv(jpr_ity1)%z3(54,1,1) = ',frcv(jpr_ity1)%z3(54,1,1) 
     1604         WRITE(968,*) 'p_taui in sbc_cpl_ice_tau = ',p_taui(54,1) 
     1605         WRITE(968,*) 'p_tauj in sbc_cpl_ice_tau = ',p_tauj(54,1) 
     1606         WRITE(968,*) 'max p_taui in sbc_cpl_ice_tau = ',MAXVAL(ABS(p_taui(:,:))) 
     1607         WRITE(968,*) 'max p_tauj in sbc_cpl_ice_tau = ',MAXVAL(ABS(p_tauj(:,:))) 
     1608      ENDIF 
     1609 
    15981610      ! 
    15991611   END SUBROUTINE sbc_cpl_ice_tau 
    … …  
    17531765      ENDIF 
    17541766 
    1755       IF(narea == 419) THEN 
    1756          WRITE(9419,*) 'emp_tot(48,31) = ',emp_tot(48,31) 
    1757          WRITE(9419,*) 'emp_ice(48,31) = ',emp_ice(48,31) 
    1758          WRITE(9419,*) 'emp_oce(48,31) = ',emp_oce(48,31) 
    1759          WRITE(9419,*) 'sprecip(48,31) = ',sprecip(48,31) 
    1760          WRITE(9419,*) 'tprecip(48,31) = ',tprecip(48,31) 
    1761          WRITE(9419,*) 'evap_ice(48,31) = ',evap_ice(48,31,:) 
    1762          WRITE(9419,*) 'devap_ice(48,31) = ',devap_ice(48,31,:) 
     1767      IF(narea == 68) THEN 
     1768         WRITE(968,*) 'sum emp_tot = ',SUM(emp_tot) 
     1769         WRITE(968,*) 'emp_ice = ',SUM(emp_ice) 
     1770         WRITE(968,*) 'emp_oce = ',SUM(emp_oce) 
     1771         WRITE(968,*) 'sprecip = ',SUM(sprecip) 
     1772         WRITE(968,*) 'tprecip = ',SUM(tprecip) 
     1773         WRITE(968,*) 'evap_ice = ',SUM(evap_ice) 
     1774         WRITE(968,*) 'devap_ice = ',SUM(devap_ice) 
    17631775      ENDIF 
    17641776 
    … …  
    18141826      CASE( 'oce only' )         ! the required field is directly provided 
    18151827         zqns_tot(:,:) = frcv(jpr_qnsoce)%z3(:,:,1) 
     1828         DO jl=1,jpl 
     1829            zqns_ice(:,:,jl) = zqns_tot(:,:) * a_i(:,:,jl) 
     1830         ENDDO 
    18161831      CASE( 'conservative' )     ! the required fields are directly provided 
    18171832         zqns_tot(:,:) = frcv(jpr_qnsmix)%z3(:,:,1) 
    … …  
    18561871      zqns_oce = 0._wp 
    18571872      WHERE( ziceld /= 0._wp )   zqns_oce(:,:) = ( zqns_tot(:,:) - SUM( a_i * zqns_ice, dim=3 ) ) / ziceld(:,:) 
     1873 
     1874      IF(narea == 68) THEN 
     1875         WRITE(968,*) 'sum zqns_oce in sbccpl = ',SUM(zqns_oce) 
     1876         WRITE(968,*) 'sum zqns_tot = ',SUM(zqns_tot) 
     1877         WRITE(968,*) 'sum zqns_ice = ',SUM(zqns_ice) 
     1878         WRITE(968,*) 'sum ziceld = ',SUM(ziceld) 
     1879      ENDIF 
    18581880 
    18591881      ! Heat content per unit mass of snow (J/kg) 
    … …  
    19041926      ENDIF 
    19051927 
    1906       IF(narea == 419) THEN 
    1907          WRITE(9419,*) 'qns_tot(48,31) = ',qns_tot(48,31) 
    1908          WRITE(9419,*) 'qns_oce(48,31) = ',qns_oce(48,31) 
    1909          WRITE(9419,*) 'qns_ice(48,31) = ',qns_ice(48,31,:) 
    1910          WRITE(9419,*) 'qevap_ice(48,31) = ',qevap_ice(48,31,:) 
    1911          WRITE(9419,*) 'qprec_ice(48,31) = ',qprec_ice(48,31) 
    1912          WRITE(9419,*) 'qemp_oce(48,31) = ',qemp_oce(48,31) 
    1913          WRITE(9419,*) 'qemp_ice(48,31) = ',qemp_ice(48,31) 
     1928      IF(narea == 68) THEN 
     1929         WRITE(968,*) 'qns_tot) = ',SUM(qns_tot) 
     1930         WRITE(968,*) 'qns_oce) = ',SUM(qns_oce) 
     1931         WRITE(968,*) 'qns_ice) = ',SUM(qns_ice) 
     1932         WRITE(968,*) 'qevap_ice) = ',SUM(qevap_ice) 
     1933         WRITE(968,*) 'qprec_ice) = ',SUM(qprec_ice) 
     1934         WRITE(968,*) 'qemp_oce) = ',SUM(qemp_oce) 
     1935         WRITE(968,*) 'qemp_ice) = ',SUM(qemp_ice) 
    19141936      ENDIF 
    19151937 
    … …  
    19541976      CASE( 'oce only' ) 
    19551977         zqsr_tot(:,:  ) = MAX( 0._wp , frcv(jpr_qsroce)%z3(:,:,1) ) 
    1956          zqsr_ice(:,:,:) = 0.0 
     1978         DO jl=1,jpl 
     1979            zqsr_ice(:,:,jl) = zqsr_tot(:,:) * a_i(:,:,jl) 
     1980         ENDDO 
    19571981      CASE( 'conservative' ) 
    19581982         zqsr_tot(:,:  ) = frcv(jpr_qsrmix)%z3(:,:,1) 
    … …  
    19812005            END DO 
    19822006 
    1983             IF(narea == 419) THEN 
    1984                WRITE(9419,*) 'zqsr_ice(48,31) oce and ice = ',zqsr_ice(48,31,:) 
     2007            IF(narea == 68) THEN 
     2008               WRITE(968,*) 'zqsr_ice) oce and ice = ',SUM(zqsr_ice) 
    19852009            ENDIF 
    19862010 
    … …  
    20262050      ENDIF 
    20272051 
    2028       IF(narea == 419) THEN 
    2029          WRITE(9419,*) 'qsr_tot(48,31) = ',qsr_tot(48,31) 
    2030          WRITE(9419,*) 'qsr_ice(48,31) = ',qsr_ice(48,31,:) 
     2052      IF(narea == 68) THEN 
     2053         WRITE(968,*) 'qsr_tot) = ',SUM(qsr_tot) 
     2054         WRITE(968,*) 'qsr_ice) = ',SUM(qsr_ice) 
    20312055      ENDIF 
    20322056 
    … …  
    20542078      ENDIF 
    20552079 
    2056       IF(narea == 419) THEN 
    2057          WRITE(9419,*) 'dqns_ice(48,31) = ',dqns_ice(48,31,:) 
     2080      IF(narea == 68) THEN 
     2081         WRITE(968,*) 'dqns_ice) = ',SUM(dqns_ice) 
    20582082      ENDIF 
    20592083 
    … …  
    20682092      END SELECT 
    20692093 
    2070       IF(narea == 419) THEN 
    2071          WRITE(9419,*) 'qml_ice(48,31) = ',qml_ice(48,31,:) 
    2072          WRITE(9419,*) 'qcn_ice(48,31) = ',qcn_ice(48,31,:) 
     2094      IF(narea == 68) THEN 
     2095         WRITE(968,*) 'qml_ice) = ',SUM(qml_ice) 
     2096         WRITE(968,*) 'qcn_ice) = ',SUM(qcn_ice) 
    20732097      ENDIF 
    20742098 
    … …  
    20942118      ENDIF 
    20952119 
    2096       IF(narea == 419) THEN 
    2097          WRITE(9419,*) 'qtr_ice_top(48,31) = ',qtr_ice_top(48,31,:) 
     2120      IF(narea == 68) THEN 
     2121         WRITE(968,*) 'qtr_ice_top) = ',SUM(qtr_ice_top) 
    20982122      ENDIF 
    20992123 
Note: See TracChangeset for help on using the changeset viewer.