Changeset 298 for trunk/NEMO/OPA_SRC/SBC/bulk.F90
- Timestamp:
- 2005-09-22T12:58:15+02:00 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/bulk.F90
r258 r298 49 49 50 50 !! * Local declarations 51 REAL(wp), DIMENSION(jpi,jpj) :: & 52 zsst , zsss 51 REAL(wp), DIMENSION(jpi,jpj) :: zsst 53 52 # if ! defined key_ice_lim 54 53 INTEGER :: ji, jj ! dummy loop indices … … 67 66 IF ( .NOT.ln_rstart ) THEN 68 67 gsst(:,:) = ( nfbulk - 1 ) * ( tn(:,:,1) + rt0 ) 69 gsss(:,:) = ( nfbulk - 1 ) * sn(:,:,1)70 68 ENDIF 71 69 ENDIF … … 84 82 85 83 gsst(:,:) = gsst(:,:) + tn(:,:,1) + rt0 86 gsss(:,:) = gsss(:,:) + sn(:,:,1)87 84 88 85 ! Computation of the fluxes 89 86 IF( MOD( kt - 1 , nfbulk ) == 0 ) THEN 90 87 91 !i zsst(:,:) = gsst(:,:) / REAL( nfbulk )92 88 zsst(:,:) = gsst(:,:) / REAL( nfbulk ) * tmask(:,:,1) 93 zsss(:,:) = gsss(:,:) / REAL( nfbulk )94 89 CALL flx_blk( zsst ) 95 90 96 91 gsst(:,:) = 0. 97 gsss(:,:) = 0.98 92 99 93 # if ! defined key_ice_lim … … 110 104 CALL prt_ctl(tab2d_1=tauy , clinfo1=' v-stress : ', mask1=tmask, ovlap=1) 111 105 CALL prt_ctl(tab2d_1=zsst , clinfo1=' sst : ', mask1=tmask, ovlap=1) 112 CALL prt_ctl(tab2d_1=zsss , clinfo1=' sss : ', mask1=tmask, ovlap=1)113 106 ENDIF 114 107 # endif
Note: See TracChangeset
for help on using the changeset viewer.