Changeset 8970 for branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/TRA
- Timestamp:
- 2017-12-11T10:01:56+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90
r8882 r8970 160 160 ! 161 161 CASE ( 1 ) !* constant flux 162 IF(lwp) WRITE(numout,*) ' ***constant heat flux = ', rn_geoflx_cst162 IF(lwp) WRITE(numout,*) ' ===>> constant heat flux = ', rn_geoflx_cst 163 163 qgh_trd0(:,:) = r1_rau0_rcp * rn_geoflx_cst 164 164 ! 165 165 CASE ( 2 ) !* variable geothermal heat flux : read the geothermal fluxes in mW/m2 166 IF(lwp) WRITE(numout,*) ' ***variable geothermal heat flux'166 IF(lwp) WRITE(numout,*) ' ===>> variable geothermal heat flux' 167 167 ! 168 168 ALLOCATE( sf_qgh(1), STAT=ierror ) … … 172 172 ENDIF 173 173 ALLOCATE( sf_qgh(1)%fnow(jpi,jpj,1) ) 174 IF( sn_qgh%ln_tint ) ALLOCATE( sf_qgh(1)%fdta(jpi,jpj,1,2) )174 IF( sn_qgh%ln_tint ) ALLOCATE( sf_qgh(1)%fdta(jpi,jpj,1,2) ) 175 175 ! fill sf_chl with sn_chl and control print 176 176 CALL fld_fill( sf_qgh, (/ sn_qgh /), cn_dir, 'tra_bbc_init', & … … 186 186 ! 187 187 ELSE 188 IF(lwp) WRITE(numout,*) ' ***no geothermal heat flux'188 IF(lwp) WRITE(numout,*) ' ===>> no geothermal heat flux' 189 189 ENDIF 190 190 !
Note: See TracChangeset
for help on using the changeset viewer.