Changeset 2798 for branches/2011
- Timestamp:
- 2011-07-11T16:18:18+02:00 (13 years ago)
- Location:
- branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM
- Files:
-
- 1 added
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_if.F90
r2715 r2798 65 65 !!--------------------------------------------------------------------- 66 66 ! ! ====================== ! 67 IF( kt == nit000 ) THEN ! First call kt=nit000 !68 ! ! ====================== !69 ! set file information70 cn_dir = './' ! directory in which the model is executed71 ! ... default values (NB: frequency positive => hours, negative => months)72 ! ! file ! frequency ! variable ! time intep ! clim ! 'yearly' or ! weights ! rotation !73 ! ! name ! (hours) ! name ! (T/F) ! (T/F) ! 'monthly' ! filename ! pairs !74 sn_ice = FLD_N('ice_cover', -1 , 'ice_cov' , .true. , .true. , 'yearly' , '' , '' )75 76 REWIND ( numnam ) ! ... read in namlist namiif77 READ ( numnam, namsbc_iif )78 79 ALLOCATE( sf_ice(1), STAT=ierror )80 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_ice_if: unable to allocate sf_ice structure' )81 ALLOCATE( sf_ice(1)%fnow(jpi,jpj,1) )82 IF( sn_ice%ln_tint ) ALLOCATE( sf_ice(1)%fdta(jpi,jpj,1,2) )83 84 ! fill sf_ice with sn_ice and control print85 CALL fld_fill( sf_ice, (/ sn_ice /), cn_dir, 'sbc_ice_if', 'ice-if sea-ice model', 'namsbc_iif' )86 !87 ENDIF88 89 CALL fld_read( kt, nn_fsbc, sf_ice ) ! Read input fields and provides the90 ! ! input fields at the current time-step91 67 92 68 IF( MOD( kt-1, nn_fsbc) == 0 ) THEN … … 104 80 ! 105 81 zt_fzp = fr_i(ji,jj) ! freezing point temperature 106 zfr_obs = sf_ice(1)%fnow(ji,jj,1) ! observed ice cover107 82 ! ! ocean ice fraction (0/1) from the freezing point temperature 108 83 IF( sst_m(ji,jj) <= zt_fzp ) THEN ; fr_i(ji,jj) = 1.e0 … … 110 85 ENDIF 111 86 87 zfr_obs = fr_i(ji,jj) 112 88 tn(ji,jj,1) = MAX( tn(ji,jj,1), zt_fzp ) ! avoid over-freezing point temperature 113 89 -
branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90
r2784 r2798 58 58 REAL(wp) :: r1_rau0 ! = 1 / rau0 59 59 60 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_rnf ! structure: river runoff (file information, fields read)60 TYPE(FLD), PUBLIC, ALLOCATABLE, DIMENSION(:) :: sf_rnf ! structure: river runoff (file information, fields read) 61 61 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_s_rnf ! structure: river runoff salinity (file information, fields read) 62 62 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_t_rnf ! structure: river runoff temperature (file information, fields read) -
branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90
r2715 r2798 26 26 REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:,:) :: etot3 !: light absortion coefficient 27 27 28 #if defined key_top && defined key_pisces28 #if defined key_top && (defined key_pisces || defined key_bfm) 29 29 !!---------------------------------------------------------------------- 30 30 !! 'key_top' & 'key_pisces' PISCES bio-model -
branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/par_trc.F90
r2528 r2798 17 17 USE par_cfc ! CFC 11 and 12 tracers 18 18 USE par_my_trc ! user defined passive tracers 19 19 USE par_bfm ! BFM biogeochemical tracers 20 20 21 IMPLICIT NONE 21 22 22 23 ! Passive tracers : Total size 23 24 ! --------------- ! total number of passive tracers, of 2d and 3d output and trend arrays 24 INTEGER, PUBLIC, PARAMETER :: jptra = jp_lobster + jp_pisces + jp_cfc + jp_c14b + jp_my_trc 25 INTEGER, PUBLIC, PARAMETER :: jpdia2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d + jp_c14b_2d + jp_my_trc_2d 26 INTEGER, PUBLIC, PARAMETER :: jpdia3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d + jp_c14b_3d + jp_my_trc_3d 25 INTEGER, PUBLIC, PARAMETER :: jptra = jp_lobster + jp_pisces + jp_cfc + jp_c14b + jp_my_trc + jp_bfm 26 INTEGER, PUBLIC, PARAMETER :: jpdia2d = jp_lobster_2d + jp_pisces_2d + jp_cfc_2d + jp_c14b_2d + jp_my_trc_2d + jp_bfm_2d 27 INTEGER, PUBLIC, PARAMETER :: jpdia3d = jp_lobster_3d + jp_pisces_3d + jp_cfc_3d + jp_c14b_3d + jp_my_trc_3d + jp_bfm_3d 27 28 ! ! total number of sms diagnostic arrays 28 INTEGER, PUBLIC, PARAMETER :: jpdiabio = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_c14b_trd + jp_my_trc_trd 29 INTEGER, PUBLIC, PARAMETER :: jpdiabio = jp_lobster_trd + jp_pisces_trd + jp_cfc_trd + jp_c14b_trd + jp_my_trc_trd +jp_bfm_trd 29 30 30 31 ! 1D configuration ("key_c1d") -
branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/trcini.F90
r2715 r2798 25 25 USE trcini_c14b ! C14 bomb initialisation 26 26 USE trcini_my_trc ! MY_TRC initialisation 27 #if defined key_bfm 28 USE par_bfm 29 #endif 27 30 USE trcdta 28 31 USE daymod … … 94 97 ENDIF 95 98 96 IF( ln_dm2dc .AND. ( lk_pisces .OR. lk_lobster ) ) &97 & CALL ctl_stop( ' The diurnal cycle is not compatible with PISCES or LOBSTER ' )99 IF( ln_dm2dc .AND. ( lk_pisces .OR. lk_lobster .OR. lk_bfm) ) & 100 & CALL ctl_stop( ' The diurnal cycle is not compatible with PISCES or LOBSTER or BFM ' ) 98 101 99 102 IF( nn_cla == 1 ) & … … 118 121 IF( lk_my_trc ) THEN ; CALL trc_ini_my_trc ! MY_TRC tracers 119 122 ELSE ; IF(lwp) WRITE(numout,*) ' MY_TRC not used' 123 ENDIF 124 125 IF( lk_bfm ) THEN ; CALL trc_ini_bfm ! BFM tracers 126 ELSE ; IF(lwp) WRITE(numout,*) ' BFM not used' 120 127 ENDIF 121 128 -
branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/trcnam.F90
r2715 r2798 26 26 USE trcnam_c14b ! C14 SMS namelist 27 27 USE trcnam_my_trc ! MY_TRC SMS namelist 28 USE trcnam_bfm ! BFM namelist 28 29 USE trdmod_trc_oce 29 30 … … 189 190 ELSE ; IF(lwp) WRITE(numout,*) ' MY_TRC not used' 190 191 ENDIF 192 193 IF( lk_bfm ) THEN ; CALL trc_nam_bfm ! BFM tracers 194 ELSE ; IF(lwp) WRITE(numout,*) ' BFM not used' 195 ENDIF 191 196 ! 192 197 END SUBROUTINE trc_nam -
branches/2011/dev_r2784_CMCC1_topbfm/NEMOGCM/NEMO/TOP_SRC/trcstp.F90
r2528 r2798 11 11 !!---------------------------------------------------------------------- 12 12 USE oce_trc ! ocean dynamics and active tracers variables 13 #if ! defined key_bfm 13 14 USE trc 14 15 USE trctrp ! passive tracers transport … … 20 21 USE trdmod_trc_oce 21 22 USE trdmld_trc 23 #else 24 USE api_bfm, ONLY: bio_calc 25 USE par_bfm 26 #ifdef key_obc 27 USE trcobcdta_bfm 28 #endif 29 #endif 22 30 USE iom 23 31 USE in_out_manager … … 49 57 !!------------------------------------------------------------------- 50 58 59 #if defined key_bfm 60 !--------------------------------------------- 61 ! Check the main BFM flag 62 !--------------------------------------------- 63 IF (bio_calc) THEN 64 65 CALL trc_bfm( kt ) ! main call to BFM 66 67 CALL trc_trp_bfm( kt ) ! transport of BFM tracers 68 69 CALL trc_dia_bfm( kt ) ! diagnostic output for BFM 70 END IF 71 #else 72 51 73 IF( MOD( kt - 1 , nn_dttrc ) == 0 ) THEN ! only every nn_dttrc time step 52 74 ! … … 71 93 ! 72 94 ENDIF 95 #endif 73 96 74 97 END SUBROUTINE trc_stp
Note: See TracChangeset
for help on using the changeset viewer.