- Timestamp:
- 2010-11-10T08:48:54+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r2291 r2370 75 75 !! NEMO/OPA 3.3 , NEMO-consortium (2010) 76 76 !! $Id$ 77 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)77 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 78 78 !!---------------------------------------------------------------------- 79 79 CONTAINS … … 133 133 ! 134 134 ! (NB: frequency positive => hours, negative => months) 135 ! ! file ! frequency ! variable ! time intep ! clim ! 'yearly' or ! weights ! rotation!136 ! ! name ! (hours) ! name ! (T/F) ! (T/F) ! 'monthly' ! filename ! pairs!137 sn_wndi = FLD_N( 'uwnd10m' , 24 , 'u_10' , .false. , .false. , 'yearly' , '' , '')138 sn_wndj = FLD_N( 'vwnd10m' , 24 , 'v_10' , .false. , .false. , 'yearly' , '' , '')139 sn_qsr = FLD_N( 'qsw' , 24 , 'qsw' , .false. , .false. , 'yearly' , '' , '')140 sn_qlw = FLD_N( 'qlw' , 24 , 'qlw' , .false. , .false. , 'yearly' , '' , '')141 sn_tair = FLD_N( 'tair10m' , 24 , 't_10' , .false. , .false. , 'yearly' , '' , '')142 sn_humi = FLD_N( 'humi10m' , 24 , 'q_10' , .false. , .false. , 'yearly' , '' , '')143 sn_prec = FLD_N( 'precip' , -1 , 'precip' , .true. , .false. , 'yearly' , '' , '')144 sn_snow = FLD_N( 'snow' , -1 , 'snow' , .true. , .false. , 'yearly' , '' , '')145 sn_tdif = FLD_N( 'taudif' , 24 , 'taudif' , .true. , .false. , 'yearly' , '' , '')135 ! ! file ! frequency ! variable ! time intep ! clim ! 'yearly' or ! weights ! rotation ! 136 ! ! name ! (hours) ! name ! (T/F) ! (T/F) ! 'monthly' ! filename ! pairs ! 137 sn_wndi = FLD_N( 'uwnd10m', 24 , 'u_10' , .false. , .false. , 'yearly' , '' , '' ) 138 sn_wndj = FLD_N( 'vwnd10m', 24 , 'v_10' , .false. , .false. , 'yearly' , '' , '' ) 139 sn_qsr = FLD_N( 'qsw' , 24 , 'qsw' , .false. , .false. , 'yearly' , '' , '' ) 140 sn_qlw = FLD_N( 'qlw' , 24 , 'qlw' , .false. , .false. , 'yearly' , '' , '' ) 141 sn_tair = FLD_N( 'tair10m', 24 , 't_10' , .false. , .false. , 'yearly' , '' , '' ) 142 sn_humi = FLD_N( 'humi10m', 24 , 'q_10' , .false. , .false. , 'yearly' , '' , '' ) 143 sn_prec = FLD_N( 'precip' , -1 , 'precip' , .true. , .false. , 'yearly' , '' , '' ) 144 sn_snow = FLD_N( 'snow' , -1 , 'snow' , .true. , .false. , 'yearly' , '' , '' ) 145 sn_tdif = FLD_N( 'taudif' , 24 , 'taudif' , .true. , .false. , 'yearly' , '' , '' ) 146 146 ! 147 147 REWIND( numnam ) ! read in namlist namsbc_core … … 396 396 !! caution : the net upward water flux has with mm/day unit 397 397 !!--------------------------------------------------------------------- 398 REAL(wp), INTENT(in ), DIMENSION(:,:,:) 399 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) 400 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) 401 REAL(wp), INTENT(in ), DIMENSION(:,:,:) 402 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) 403 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) 404 REAL(wp), INTENT( out), DIMENSION(:,:,:) 405 REAL(wp), INTENT( out), DIMENSION(:,:,:) 406 REAL(wp), INTENT( out), DIMENSION(:,:,:) 407 REAL(wp), INTENT( out), DIMENSION(:,:,:) 408 REAL(wp), INTENT( out), DIMENSION(:,:,:) 409 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) 410 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) 411 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) 412 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) 413 CHARACTER(len=1), INTENT(in ) 414 INTEGER, INTENT(in ) 398 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: pst ! ice surface temperature (>0, =rt0 over land) [Kelvin] 399 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pui ! ice surface velocity (i- and i- components [m/s] 400 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pvi ! at I-point (B-grid) or U & V-point (C-grid) 401 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: palb ! ice albedo (clear sky) (alb_ice_cs) [%] 402 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_taui ! i- & j-components of surface ice stress [N/m2] 403 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tauj ! at I-point (B-grid) or U & V-point (C-grid) 404 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_qns ! non solar heat flux over ice (T-point) [W/m2] 405 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_qsr ! solar heat flux over ice (T-point) [W/m2] 406 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_qla ! latent heat flux over ice (T-point) [W/m2] 407 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_dqns ! non solar heat sensistivity (T-point) [W/m2] 408 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_dqla ! latent heat sensistivity (T-point) [W/m2] 409 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tpr ! total precipitation (T-point) [Kg/m2/s] 410 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_spr ! solid precipitation (T-point) [Kg/m2/s] 411 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr1 ! 1sr fraction of qsr penetration in ice (T-point) [%] 412 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr2 ! 2nd fraction of qsr penetration in ice (T-point) [%] 413 CHARACTER(len=1), INTENT(in ) :: cd_grid ! ice grid ( C or B-grid) 414 INTEGER, INTENT(in ) :: pdim ! number of ice categories 415 415 !! 416 416 INTEGER :: ji, jj, jl ! dummy loop indices … … 449 449 ! ----------------------------------------------------------------------------- ! 450 450 SELECT CASE( cd_grid ) 451 CASE( ' B' ) ! B-grid ice dynamics : I-point (i.e. F-point with sea-ice indexation)451 CASE( 'I' ) ! B-grid ice dynamics : I-point (i.e. F-point with sea-ice indexation) 452 452 ! and scalar wind at T-point ( = | U10m - U_ice | ) (masked) 453 453 #if defined key_vectopt_loop … … 605 605 !! 9.0 ! 05-08 (L. Brodeau) Rewriting and optimization 606 606 !!---------------------------------------------------------------------- 607 !! * Arguments608 609 607 REAL(wp), INTENT(in) :: zu ! altitude of wind measurement [m] 610 608 REAL(wp), INTENT(in), DIMENSION(jpi,jpj) :: & … … 646 644 grav = 9.8, & ! gravity 647 645 kappa = 0.4 ! von Karman s constant 648 646 !!---------------------------------------------------------------------- 649 647 !! * Start 650 648 !! Air/sea differences … … 770 768 grav = 9.8, & ! gravity 771 769 kappa = 0.4 ! von Karman's constant 772 770 !!---------------------------------------------------------------------- 773 771 !! * Start 774 772
Note: See TracChangeset
for help on using the changeset viewer.