Changeset 4 for vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC
- Timestamp:
- 02/08/12 17:17:04 (12 years ago)
- Location:
- vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 21 edited
Legend:
- Unmodified
- Added
- Removed
-
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90
r1 r4 48 48 !!---------------------------------------------------------------------- 49 49 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 50 !! $Id: albedo.F90 32 09 2011-12-10 08:04:49Z smasson$50 !! $Id: albedo.F90 3294 2012-01-28 16:44:18Z rblod $ 51 51 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 52 52 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90
r1 r4 68 68 !!---------------------------------------------------------------------- 69 69 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 70 !! $Id: cpl_oasis3.F90 3 116 2011-11-15 20:55:40Z cetlod $70 !! $Id: cpl_oasis3.F90 3294 2012-01-28 16:44:18Z rblod $ 71 71 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 72 72 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis4.F90
r1 r4 69 69 !!---------------------------------------------------------------------- 70 70 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 71 !! $Id: cpl_oasis4.F90 3 186 2011-11-27 08:16:19Z smasson$71 !! $Id: cpl_oasis4.F90 3294 2012-01-28 16:44:18Z rblod $ 72 72 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 73 73 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/fldread.F90
r1 r4 100 100 !!---------------------------------------------------------------------- 101 101 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 102 !! $Id: fldread.F90 3 191 2011-12-05 09:41:41Z davestorkey$102 !! $Id: fldread.F90 3294 2012-01-28 16:44:18Z rblod $ 103 103 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 104 104 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_ice.F90
r1 r4 100 100 !!---------------------------------------------------------------------- 101 101 !! NEMO/OPA 4.0 , NEMO Consortium (2011) 102 !! $Id: sbc_ice.F90 3 116 2011-11-15 20:55:40Z cetlod $102 !! $Id: sbc_ice.F90 3294 2012-01-28 16:44:18Z rblod $ 103 103 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 104 104 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90
r1 r4 89 89 !!---------------------------------------------------------------------- 90 90 !! NEMO/OPA 4.0 , NEMO Consortium (2011) 91 !! $Id: sbc_oce.F90 3 105 2011-11-15 11:33:20Z poddo$91 !! $Id: sbc_oce.F90 3294 2012-01-28 16:44:18Z rblod $ 92 92 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 93 93 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90
r1 r4 41 41 !!---------------------------------------------------------------------- 42 42 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 43 !! $Id: sbcana.F90 2977 2011-10-22 13:46:41Z cetlod $43 !! $Id: sbcana.F90 3294 2012-01-28 16:44:18Z rblod $ 44 44 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 45 45 !!---------------------------------------------------------------------- … … 89 89 nn_tau000 = MAX( nn_tau000, 1 ) ! must be >= 1 90 90 ! 91 qns (:,:) = rn_qns0 92 qsr (:,:) = rn_qsr0 93 emp (:,:) = rn_emp0 94 emps(:,:) = rn_emp0 95 ! 96 utau(:,:) = rn_utau0 97 vtau(:,:) = rn_vtau0 98 taum(:,:) = SQRT ( rn_utau0 * rn_utau0 + rn_vtau0 * rn_vtau0 ) 99 wndm(:,:) = SQRT ( taum(1,1) / ( zrhoa * zcdrag ) ) 100 ! 91 101 ENDIF 92 102 93 qns (:,:) = rn_qns094 qsr (:,:) = rn_qsr095 emp (:,:) = rn_emp096 emps(:,:) = rn_emp097 103 98 104 ! Increase the surface stress to its nominal value during the first nn_tau000 time-steps -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r1 r4 28 28 USE lib_mpp ! distribued memory computing library 29 29 USE wrk_nemo ! work arrays 30 USE timing ! Timing 30 31 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 31 32 … … 84 85 !!---------------------------------------------------------------------- 85 86 !! NEMO/OPA 4.0 , NEMO Consortium (2011) 86 !! $Id: sbcblk_clio.F90 3 186 2011-11-27 08:16:19Z smasson$87 !! $Id: sbcblk_clio.F90 3294 2012-01-28 16:44:18Z rblod $ 87 88 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 88 89 !!---------------------------------------------------------------------- … … 227 228 REAL(wp), POINTER, DIMENSION(:,:) :: zqsb ! sensible heat flux over ocean 228 229 !!--------------------------------------------------------------------- 229 230 ! 231 IF( nn_timing == 1 ) CALL timing_start('blk_oce_clio') 232 ! 230 233 CALL wrk_alloc( jpi,jpj, zqlw, zqla, zqsb ) 231 234 … … 381 384 CALL wrk_dealloc( jpi,jpj, zqlw, zqla, zqsb ) 382 385 ! 386 IF( nn_timing == 1 ) CALL timing_stop('blk_oce_clio') 387 ! 383 388 END SUBROUTINE blk_oce_clio 384 389 … … 444 449 REAL(wp), DIMENSION(:,:,:), POINTER :: z_qlw, z_qsb 445 450 !!--------------------------------------------------------------------- 446 451 ! 452 IF( nn_timing == 1 ) CALL timing_start('blk_ice_clio') 453 ! 447 454 CALL wrk_alloc( jpi,jpj, ztatm, zqatm, zevsqr, zrhoa ) 448 455 CALL wrk_alloc( jpi,jpj,pdim, z_qlw, z_qsb ) … … 625 632 CALL wrk_dealloc( jpi,jpj,pdim, z_qlw, z_qsb ) 626 633 ! 634 IF( nn_timing == 1 ) CALL timing_stop('blk_ice_clio') 635 ! 627 636 END SUBROUTINE blk_ice_clio 628 637 … … 661 670 REAL(wp), DIMENSION(:,:), POINTER :: zps, zpc ! sine (cosine) of latitude per sine (cosine) of solar declination 662 671 !!--------------------------------------------------------------------- 663 672 ! 673 IF( nn_timing == 1 ) CALL timing_start('blk_clio_qsr_oce') 674 ! 664 675 CALL wrk_alloc( jpi,jpj, zev, zdlha, zlsrise, zlsset, zps, zpc ) 665 676 … … 778 789 CALL wrk_dealloc( jpi,jpj, zev, zdlha, zlsrise, zlsset, zps, zpc ) 779 790 ! 791 IF( nn_timing == 1 ) CALL timing_stop('blk_clio_qsr_oce') 792 ! 780 793 END SUBROUTINE blk_clio_qsr_oce 781 794 … … 812 825 REAL(wp), DIMENSION(:,:), POINTER :: zdlha, zlsrise, zlsset ! 2D workspace 813 826 REAL(wp), DIMENSION(:,:), POINTER :: zps, zpc ! sine (cosine) of latitude per sine (cosine) of solar declination 814 !!--------------------------------------------------------------------- 815 827 !!--------------------------------------------------------------------- 828 ! 829 IF( nn_timing == 1 ) CALL timing_start('blk_clio_qsr_ice') 830 ! 816 831 CALL wrk_alloc( jpi,jpj, zev, zdlha, zlsrise, zlsset, zps, zpc ) 817 832 … … 919 934 CALL wrk_dealloc( jpi,jpj, zev, zdlha, zlsrise, zlsset, zps, zpc ) 920 935 ! 936 IF( nn_timing == 1 ) CALL timing_stop('blk_clio_qsr_ice') 937 ! 921 938 END SUBROUTINE blk_clio_qsr_ice 922 939 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r1 r4 34 34 USE lib_mpp ! distribued memory computing library 35 35 USE wrk_nemo ! work arrays 36 USE timing ! Timing 36 37 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 37 38 USE prtctl ! Print control … … 80 81 !!---------------------------------------------------------------------- 81 82 !! NEMO/OPA 3.3 , NEMO-consortium (2010) 82 !! $Id: sbcblk_core.F90 3 186 2011-11-27 08:16:19Z smasson$83 !! $Id: sbcblk_core.F90 3294 2012-01-28 16:44:18Z rblod $ 83 84 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 84 85 !!---------------------------------------------------------------------- … … 242 243 REAL(wp), DIMENSION(:,:), POINTER :: zq_zu ! air spec. hum. at wind speed height 243 244 !!--------------------------------------------------------------------- 244 245 ! 246 IF( nn_timing == 1 ) CALL timing_start('blk_oce_core') 247 ! 245 248 CALL wrk_alloc( jpi,jpj, zwnd_i, zwnd_j, zqsatw, zqlw, zqsb, zqla, zevap ) 246 249 CALL wrk_alloc( jpi,jpj, Cd, Ch, Ce, zst, zt_zu, zq_zu ) … … 396 399 CALL wrk_dealloc( jpi,jpj, zwnd_i, zwnd_j, zqsatw, zqlw, zqsb, zqla, zevap ) 397 400 CALL wrk_dealloc( jpi,jpj, Cd, Ch, Ce, zst, zt_zu, zq_zu ) 401 ! 402 IF( nn_timing == 1 ) CALL timing_stop('blk_oce_core') 398 403 ! 399 404 END SUBROUTINE blk_oce_core … … 450 455 REAL(wp), DIMENSION(:,:,:), POINTER :: z_dqsb ! sensible heat sensitivity over ice 451 456 !!--------------------------------------------------------------------- 452 457 ! 458 IF( nn_timing == 1 ) CALL timing_start('blk_ice_core') 459 ! 453 460 CALL wrk_alloc( jpi,jpj, z_wnds_t ) 454 461 CALL wrk_alloc( jpi,jpj,pdim, z_qlw, z_qsb, z_dqlw, z_dqsb ) … … 610 617 CALL wrk_dealloc( jpi,jpj,pdim, z_qlw, z_qsb, z_dqlw, z_dqsb ) 611 618 ! 619 IF( nn_timing == 1 ) CALL timing_stop('blk_ice_core') 620 ! 612 621 END SUBROUTINE blk_ice_core 613 622 … … 663 672 INTEGER , DIMENSION(:,:), POINTER :: stab ! 1st guess stability test integer 664 673 !!---------------------------------------------------------------------- 665 674 ! 675 IF( nn_timing == 1 ) CALL timing_start('TURB_CORE_1Z') 676 ! 666 677 CALL wrk_alloc( jpi,jpj, stab ) ! integer 667 678 CALL wrk_alloc( jpi,jpj, dU10, dT, dq, Cd_n10, Ce_n10, Ch_n10, sqrt_Cd_n10, sqrt_Cd, L ) … … 739 750 CALL wrk_dealloc( jpi,jpj, dU10, dT, dq, Cd_n10, Ce_n10, Ch_n10, sqrt_Cd_n10, sqrt_Cd, L ) 740 751 CALL wrk_dealloc( jpi,jpj, T_vpot, T_star, q_star, U_star, zeta, U_n10, xlogt, xct, zpsi_h, zpsi_m ) 752 ! 753 IF( nn_timing == 1 ) CALL timing_stop('TURB_CORE_1Z') 741 754 ! 742 755 END SUBROUTINE TURB_CORE_1Z … … 796 809 INTEGER , DIMENSION(:,:), POINTER :: stab ! 1st stability test integer 797 810 !!---------------------------------------------------------------------- 798 811 ! 812 IF( nn_timing == 1 ) CALL timing_start('TURB_CORE_2Z') 813 ! 799 814 CALL wrk_alloc( jpi,jpj, dU10, dT, dq, Cd_n10, Ce_n10, Ch_n10, sqrt_Cd_n10, sqrt_Cd, L ) 800 815 CALL wrk_alloc( jpi,jpj, T_vpot, T_star, q_star, U_star, zeta_u, zeta_t, U_n10 ) … … 887 902 CALL wrk_dealloc( jpi,jpj, stab ) ! interger 888 903 ! 904 IF( nn_timing == 1 ) CALL timing_stop('TURB_CORE_2Z') 905 ! 889 906 END SUBROUTINE TURB_CORE_2Z 890 907 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90
r1 r4 21 21 USE lib_mpp ! distribued memory computing library 22 22 USE wrk_nemo ! work arrays 23 USE timing ! Timing 23 24 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 24 25 USE prtctl ! Print control … … 120 121 & sn_tair , sn_rhm , sn_prec 121 122 !!--------------------------------------------------------------------- 122 123 ! 124 IF( nn_timing == 1 ) CALL timing_start('sbc_blk_mfs') 125 ! 123 126 ! ! ====================== ! 124 127 IF( kt == nit000 ) THEN ! First call kt=nit000 ! … … 263 266 264 267 ENDIF 265 268 ! 269 IF( nn_timing == 1 ) CALL timing_stop('sbc_blk_mfs') 270 ! 266 271 END SUBROUTINE sbc_blk_mfs 267 272 … … 324 329 REAL(wp), DIMENSION(5) :: p_e = (/-0.16,1.0,1.0,1.0,1.0/) 325 330 INTEGER :: kku !index varing with wind speed 326 331 ! 332 IF( nn_timing == 1 ) CALL timing_start('fluxes_mfs') 333 ! 327 334 CALL wrk_alloc( jpi,jpj, rspeed, sh10now, t10now, cdx, ce, shms ) 328 335 CALL wrk_alloc( jpi,jpj, rhom, sstk, ch, rel_windu, rel_windv ) … … 490 497 CALL wrk_dealloc( jpi,jpj, rspeed, sh10now, t10now, cdx, ce, shms ) 491 498 CALL wrk_dealloc( jpi,jpj, rhom, sstk, ch, rel_windu, rel_windv ) 492 499 ! 500 IF( nn_timing == 1 ) CALL timing_stop('fluxes_mfs') 501 ! 493 502 END SUBROUTINE fluxes_mfs 494 503 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r1 r4 48 48 USE lib_mpp ! distribued memory computing library 49 49 USE wrk_nemo ! work arrays 50 USE timing ! Timing 50 51 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 51 52 #if defined key_cpl_carbon_cycle … … 169 170 !!---------------------------------------------------------------------- 170 171 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 171 !! $Id: sbccpl.F90 3 186 2011-11-27 08:16:19Z smasson$172 !! $Id: sbccpl.F90 3294 2012-01-28 16:44:18Z rblod $ 172 173 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 173 174 !!---------------------------------------------------------------------- … … 227 228 & sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf , sn_rcv_cal , sn_rcv_iceflx , sn_rcv_co2 228 229 !!--------------------------------------------------------------------- 229 230 ! 231 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_init') 232 ! 230 233 CALL wrk_alloc( jpi,jpj, zacs, zaos ) 231 234 … … 617 620 CALL wrk_dealloc( jpi,jpj, zacs, zaos ) 618 621 ! 622 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_init') 623 ! 619 624 END SUBROUTINE sbc_cpl_init 620 625 … … 676 681 REAL(wp), POINTER, DIMENSION(:,:) :: ztx, zty 677 682 !!---------------------------------------------------------------------- 678 683 ! 684 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_rcv') 685 ! 679 686 CALL wrk_alloc( jpi,jpj, ztx, zty ) 680 687 … … 862 869 CALL wrk_dealloc( jpi,jpj, ztx, zty ) 863 870 ! 871 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_rcv') 872 ! 864 873 END SUBROUTINE sbc_cpl_rcv 865 874 … … 905 914 REAL(wp), POINTER, DIMENSION(:,:) :: ztx, zty 906 915 !!---------------------------------------------------------------------- 907 916 ! 917 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_ice_tau') 918 ! 908 919 CALL wrk_alloc( jpi,jpj, ztx, zty ) 909 920 … … 1070 1081 CALL wrk_dealloc( jpi,jpj, ztx, zty ) 1071 1082 ! 1083 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_ice_tau') 1084 ! 1072 1085 END SUBROUTINE sbc_cpl_ice_tau 1073 1086 … … 1122 1135 REAL(wp), POINTER, DIMENSION(:,:) :: zcptn, ztmp, zicefr 1123 1136 !!---------------------------------------------------------------------- 1124 1137 ! 1138 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_ice_flx') 1139 ! 1125 1140 CALL wrk_alloc( jpi,jpj, zcptn, ztmp, zicefr ) 1126 1141 … … 1301 1316 CALL wrk_dealloc( jpi,jpj, zcptn, ztmp, zicefr ) 1302 1317 ! 1318 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_ice_flx') 1319 ! 1303 1320 END SUBROUTINE sbc_cpl_ice_flx 1304 1321 … … 1320 1337 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmp3, ztmp4 1321 1338 !!---------------------------------------------------------------------- 1322 1339 ! 1340 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_snd') 1341 ! 1323 1342 CALL wrk_alloc( jpi,jpj, zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1 ) 1324 1343 CALL wrk_alloc( jpi,jpj,jpl, ztmp3, ztmp4 ) … … 1545 1564 CALL wrk_dealloc( jpi,jpj,jpl, ztmp3, ztmp4 ) 1546 1565 ! 1566 IF( nn_timing == 1 ) CALL timing_stop('sbc_cpl_snd') 1567 ! 1547 1568 END SUBROUTINE sbc_cpl_snd 1548 1569 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcdcy.F90
r1 r4 18 18 USE in_out_manager ! I/O manager 19 19 USE lib_mpp ! MPP library 20 USE timing ! Timing 20 21 21 22 IMPLICIT NONE … … 31 32 !!---------------------------------------------------------------------- 32 33 !! NEMO/OPA 3.3 , NEMO-consortium (2010) 33 !! $Id: sbcdcy.F90 2715 2011-03-30 15:58:35Z rblod $34 !! $Id: sbcdcy.F90 3294 2012-01-28 16:44:18Z rblod $ 34 35 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 35 36 !!---------------------------------------------------------------------- … … 76 77 & - paaa * pt1 - zinvtwopi * pbbb * SIN(pccc + ztwopi * pt1) 77 78 !!--------------------------------------------------------------------- 78 79 ! 80 IF( nn_timing == 1 ) CALL timing_start('sbc_dcy') 81 ! 79 82 ! Initialization 80 83 ! -------------- … … 221 224 END DO 222 225 ! 226 IF( nn_timing == 1 ) CALL timing_stop('sbc_dcy') 227 ! 223 228 END FUNCTION sbc_dcy 224 229 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90
r1 r4 23 23 USE lib_mpp ! distribued memory computing library 24 24 USE wrk_nemo ! work arrays 25 USE timing ! Timing 25 26 USE lbclnk ! ocean lateral boundary conditions 26 27 USE lib_fortran … … 41 42 !!---------------------------------------------------------------------- 42 43 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 43 !! $Id: sbcfwb.F90 3 186 2011-11-27 08:16:19Z smasson$44 !! $Id: sbcfwb.F90 3294 2012-01-28 16:44:18Z rblod $ 44 45 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 45 46 !!---------------------------------------------------------------------- … … 68 69 REAL(wp), POINTER, DIMENSION(:,:) :: ztmsk_neg, ztmsk_pos, ztmsk_tospread, z_wgt, zerp_cor 69 70 !!---------------------------------------------------------------------- 71 ! 72 IF( nn_timing == 1 ) CALL timing_start('sbc_fwb') 70 73 ! 71 74 CALL wrk_alloc( jpi,jpj, ztmsk_neg, ztmsk_pos, ztmsk_tospread, z_wgt, zerp_cor ) … … 192 195 CALL wrk_dealloc( jpi,jpj, ztmsk_neg, ztmsk_pos, ztmsk_tospread, z_wgt, zerp_cor ) 193 196 ! 197 IF( nn_timing == 1 ) CALL timing_stop('sbc_fwb') 198 ! 194 199 END SUBROUTINE sbc_fwb 195 200 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r1 r4 14 14 USE oce ! ocean dynamics and tracers 15 15 USE dom_oce ! ocean space and time domain 16 USE domvvl 16 17 USE phycst, only : rcp, rau0 17 18 USE in_out_manager ! I/O manager … … 354 355 355 356 ! Freezing/melting potential 356 ! Calculated over NEMO leapfrog timestep (hence 2*dt), so uses tb in 357 ! frzmlt which is then applied in going from tb to ta. 358 ! May be better using sst_m if not coupling to CICE every time-step 359 360 ! nfrzmlt(:,:)=rau0*rcp*fse3t(:,:,1)*(Tocnfrz-sst_m(:,:))/(2.0*dt) 361 nfrzmlt(:,:)=rau0*rcp*fse3t(:,:,1)*(Tocnfrz-tsb(:,:,1,jp_tem))/(2.0*dt) 357 ! Calculated over NEMO leapfrog timestep (hence 2*dt) 358 nfrzmlt(:,:)=rau0*rcp*fse3t_m(:,:,1)*(Tocnfrz-sst_m(:,:))/(2.0*dt) 362 359 363 360 ztmp(:,:) = nfrzmlt(:,:) … … 499 496 ! Sublimation from the ice is treated in a similar way (included in emp but not emps) 500 497 ! 501 ! This will need re-thinking in the variable volume case and if ice is 'embedded' in the 502 ! ocean rather than floating on top 503 504 emp(:,:) = emp(:,:) - tprecip(:,:)*fr_i(:,:) 498 ! This should not be done in the variable volume case 499 500 IF (.NOT. lk_vvl) THEN 501 502 emp(:,:) = emp(:,:) - tprecip(:,:)*fr_i(:,:) 505 503 506 504 ! Take sublimation into account 507 IF (nsbc == 5 ) THEN 508 emp(:,:) = emp(:,:) + ( emp_ice(:,:) + sprecip(:,:) ) 509 ELSE IF (nsbc == 2 ) THEN 510 emp(:,:) = emp(:,:) - qla_ice(:,:,1) / Lsub 505 IF (nsbc == 5 ) THEN 506 emp(:,:) = emp(:,:) + ( emp_ice(:,:) + sprecip(:,:) ) 507 ELSE IF (nsbc == 2 ) THEN 508 emp(:,:) = emp(:,:) - qla_ice(:,:,1) / Lsub 509 ENDIF 510 511 511 ENDIF 512 512 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_if.F90
r1 r4 33 33 !!---------------------------------------------------------------------- 34 34 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 35 !! $Id: sbcice_if.F90 3 116 2011-11-15 20:55:40Z cetlod $35 !! $Id: sbcice_if.F90 3294 2012-01-28 16:44:18Z rblod $ 36 36 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 37 37 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r1 r4 65 65 !!---------------------------------------------------------------------- 66 66 !! NEMO/OPA 4.0 , UCL NEMO Consortium (2011) 67 !! $Id: sbcice_lim.F90 3 186 2011-11-27 08:16:19Z smasson$67 !! $Id: sbcice_lim.F90 3294 2012-01-28 16:44:18Z rblod $ 68 68 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 69 69 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90
r1 r4 59 59 !!---------------------------------------------------------------------- 60 60 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 61 !! $Id: sbcice_lim_2.F90 3 186 2011-11-27 08:16:19Z smasson$61 !! $Id: sbcice_lim_2.F90 3294 2012-01-28 16:44:18Z rblod $ 62 62 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 63 63 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90
r1 r4 49 49 USE in_out_manager ! I/O manager 50 50 USE lib_mpp ! MPP library 51 USE timing ! Timing 51 52 USE sbcwave ! Wave module 52 53 … … 63 64 !!---------------------------------------------------------------------- 64 65 !! NEMO/OPA 4.0 , NEMO-consortium (2011) 65 !! $Id: sbcmod.F90 3 116 2011-11-15 20:55:40Z cetlod $66 !! $Id: sbcmod.F90 3294 2012-01-28 16:44:18Z rblod $ 66 67 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 67 68 !!---------------------------------------------------------------------- … … 222 223 INTEGER, INTENT(in) :: kt ! ocean time step 223 224 !!--------------------------------------------------------------------- 224 225 ! 226 IF( nn_timing == 1 ) CALL timing_start('sbc') 227 ! 225 228 ! ! ---------------------------------------- ! 226 229 IF( kt /= nit000 ) THEN ! Swap of forcing fields ! … … 363 366 IF( kt == nitend ) CALL sbc_final ! Close down surface module if necessary 364 367 ! 368 IF( nn_timing == 1 ) CALL timing_stop('sbc') 369 ! 365 370 END SUBROUTINE sbc 366 371 -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90
r1 r4 66 66 !!---------------------------------------------------------------------- 67 67 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 68 !! $Id: sbcrnf.F90 3 189 2011-11-29 13:37:29Z charris$68 !! $Id: sbcrnf.F90 3294 2012-01-28 16:44:18Z rblod $ 69 69 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 70 70 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
r1 r4 31 31 !!---------------------------------------------------------------------- 32 32 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 33 !! $Id: sbcssm.F90 2977 2011-10-22 13:46:41Z cetlod $33 !! $Id: sbcssm.F90 3294 2012-01-28 16:44:18Z rblod $ 34 34 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 35 35 !!---------------------------------------------------------------------- -
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90
r1 r4 21 21 USE lib_mpp ! distribued memory computing library 22 22 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 23 USE timing ! Timing 23 24 24 25 IMPLICIT NONE … … 46 47 !!---------------------------------------------------------------------- 47 48 !! NEMO/OPA 4.0 , NEMO Consortium (2011) 48 !! $Id: sbcssr.F90 2715 2011-03-30 15:58:35Z rblod $49 !! $Id: sbcssr.F90 3294 2012-01-28 16:44:18Z rblod $ 49 50 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 50 51 !!---------------------------------------------------------------------- … … 78 79 NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, rn_dqdt, rn_deds, sn_sst, sn_sss, ln_sssr_bnd, rn_sssr_bnd 79 80 !!---------------------------------------------------------------------- 81 ! 82 IF( nn_timing == 1 ) CALL timing_start('sbc_ssr') 80 83 ! 81 84 ! ! -------------------- ! … … 201 204 ENDIF 202 205 ! 206 IF( nn_timing == 1 ) CALL timing_stop('sbc_ssr') 207 ! 203 208 END SUBROUTINE sbc_ssr 204 209
Note: See TracChangeset
for help on using the changeset viewer.