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

Ignore:
Timestamp:
2019-07-26T11:32:13+02:00 (5 years ago)
Author:
dancopsey
Message:

Add lots of print statements printing out values of files throughout the ocean_ice timestep.. New output will be in files crash_stats_419.out and crash_stats_68.out.

File:
1 edited

Legend:

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

    r10888 r11355  
    13781378         CALL iom_put( 'frq_m', frq_m ) 
    13791379      ENDIF 
     1380 
     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 ) 
     1384      ENDIF 
    13801385       
    13811386      !                                                      ! ========================= ! 
     
    13961401            zemp(:,:) = 0._wp 
    13971402         ENDIF 
     1403 
     1404         IF(narea == 419) THEN 
     1405            WRITE(9419,*) 'max zemp after coupling = ',MAXVAL(  ABS( zemp(:,:) )  ) 
     1406         ENDIF 
     1407 
    13981408         ! 
    13991409         !                                                        ! runoffs and calving (added in emp) 
    14001410         IF( srcv(jpr_rnf)%laction )     rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 
    14011411         IF( srcv(jpr_cal)%laction )     zemp(:,:) = zemp(:,:) - frcv(jpr_cal)%z3(:,:,1) 
     1412 
     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(:,:) )  ) 
     1416         ENDIF 
    14021417  
    14031418         IF( srcv(jpr_icb)%laction )  THEN  
     
    14101425         ELSE                   ;   emp(:,:) =                              zemp(:,:) 
    14111426         ENDIF 
     1427 
     1428         IF(narea == 419) THEN 
     1429            WRITE(9419,*) 'max emp after conversion from zemp = ',MAXVAL(  ABS( emp(:,:) )  ) 
     1430         ENDIF 
     1431 
    14121432         ! 
    14131433         !                                                       ! non solar heat flux over the ocean (qns) 
     
    17331753      ENDIF 
    17341754 
     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,:) 
     1763      ENDIF 
     1764 
    17351765#else 
    17361766      zsnw(:,:) = picefr(:,:) 
     
    18741904      ENDIF 
    18751905 
     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) 
     1914      ENDIF 
     1915 
    18761916#else 
    18771917      zcptsnw (:,:) = zcptn(:,:) 
     
    19141954      CASE( 'oce only' ) 
    19151955         zqsr_tot(:,:  ) = MAX( 0._wp , frcv(jpr_qsroce)%z3(:,:,1) ) 
     1956         zqsr_ice(:,:,:) = 0.0 
    19161957      CASE( 'conservative' ) 
    19171958         zqsr_tot(:,:  ) = frcv(jpr_qsrmix)%z3(:,:,1) 
     
    19391980               zqsr_ice(:,:,jl) = frcv(jpr_qsrice)%z3(:,:,1) 
    19401981            END DO 
     1982 
     1983            IF(narea == 419) THEN 
     1984               WRITE(9419,*) 'zqsr_ice(48,31) oce and ice = ',zqsr_ice(48,31,:) 
     1985            ENDIF 
     1986 
    19411987         ENDIF 
    19421988      CASE( 'mixed oce-ice' ) 
     
    19802026      ENDIF 
    19812027 
     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,:) 
     2031      ENDIF 
     2032 
     2033 
    19822034      !                                                      ! ========================= ! 
    19832035      SELECT CASE( TRIM( sn_rcv_dqnsdt%cldes ) )             !          d(qns)/dt        ! 
     
    19992051         ENDDO 
    20002052      ELSE 
    2001          dqns_ice(:,:,:) = zdqns_ice(:,:,:) 
    2002       ENDIF 
     2053         dqns_ice(:,:,:) = 0.0 
     2054      ENDIF 
     2055 
     2056      IF(narea == 419) THEN 
     2057         WRITE(9419,*) 'dqns_ice(48,31) = ',dqns_ice(48,31,:) 
     2058      ENDIF 
     2059 
    20032060 
    20042061#if defined key_si3       
     
    20102067         qcn_ice(:,:,:) = frcv(jpr_botm)%z3(:,:,:) 
    20112068      END SELECT 
     2069 
     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,:) 
     2073      ENDIF 
     2074 
    20122075      ! 
    20132076      !                                                      ! ========================= ! 
     
    20302093         ! 
    20312094      ENDIF 
     2095 
     2096      IF(narea == 419) THEN 
     2097         WRITE(9419,*) 'qtr_ice_top(48,31) = ',qtr_ice_top(48,31,:) 
     2098      ENDIF 
     2099 
    20322100      ! 
    20332101#endif 
Note: See TracChangeset for help on using the changeset viewer.