- Timestamp:
- 2019-07-26T11:32:13+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/OCE/SBC/sbccpl.F90
r10888 r11355 1378 1378 CALL iom_put( 'frq_m', frq_m ) 1379 1379 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 1380 1385 1381 1386 ! ! ========================= ! … … 1396 1401 zemp(:,:) = 0._wp 1397 1402 ENDIF 1403 1404 IF(narea == 419) THEN 1405 WRITE(9419,*) 'max zemp after coupling = ',MAXVAL( ABS( zemp(:,:) ) ) 1406 ENDIF 1407 1398 1408 ! 1399 1409 ! ! runoffs and calving (added in emp) 1400 1410 IF( srcv(jpr_rnf)%laction ) rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 1401 1411 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 1402 1417 1403 1418 IF( srcv(jpr_icb)%laction ) THEN … … 1410 1425 ELSE ; emp(:,:) = zemp(:,:) 1411 1426 ENDIF 1427 1428 IF(narea == 419) THEN 1429 WRITE(9419,*) 'max emp after conversion from zemp = ',MAXVAL( ABS( emp(:,:) ) ) 1430 ENDIF 1431 1412 1432 ! 1413 1433 ! ! non solar heat flux over the ocean (qns) … … 1733 1753 ENDIF 1734 1754 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 1735 1765 #else 1736 1766 zsnw(:,:) = picefr(:,:) … … 1874 1904 ENDIF 1875 1905 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 1876 1916 #else 1877 1917 zcptsnw (:,:) = zcptn(:,:) … … 1914 1954 CASE( 'oce only' ) 1915 1955 zqsr_tot(:,: ) = MAX( 0._wp , frcv(jpr_qsroce)%z3(:,:,1) ) 1956 zqsr_ice(:,:,:) = 0.0 1916 1957 CASE( 'conservative' ) 1917 1958 zqsr_tot(:,: ) = frcv(jpr_qsrmix)%z3(:,:,1) … … 1939 1980 zqsr_ice(:,:,jl) = frcv(jpr_qsrice)%z3(:,:,1) 1940 1981 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 1941 1987 ENDIF 1942 1988 CASE( 'mixed oce-ice' ) … … 1980 2026 ENDIF 1981 2027 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 1982 2034 ! ! ========================= ! 1983 2035 SELECT CASE( TRIM( sn_rcv_dqnsdt%cldes ) ) ! d(qns)/dt ! … … 1999 2051 ENDDO 2000 2052 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 2003 2060 2004 2061 #if defined key_si3 … … 2010 2067 qcn_ice(:,:,:) = frcv(jpr_botm)%z3(:,:,:) 2011 2068 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 2012 2075 ! 2013 2076 ! ! ========================= ! … … 2030 2093 ! 2031 2094 ENDIF 2095 2096 IF(narea == 419) THEN 2097 WRITE(9419,*) 'qtr_ice_top(48,31) = ',qtr_ice_top(48,31,:) 2098 ENDIF 2099 2032 2100 ! 2033 2101 #endif
Note: See TracChangeset
for help on using the changeset viewer.