- Timestamp:
- 2012-10-04T18:31:06+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90
r3402 r3488 106 106 zcoef = z_fwf * rcp 107 107 emp(:,:) = emp(:,:) - z_fwf 108 qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) ! ensure fw correction does not change the heat budget108 qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) ! account for change to the heat budget due to fw correction 109 109 ENDIF 110 110 ! 111 111 CASE ( 2 ) !== fwf budget adjusted from the previous year ==! 112 112 ! 113 IF( kt == nit000 ) THEN ! initialisation114 ! ! Read the corrective factor on precipitations (fwfold)113 IF( kt == nit000 ) THEN ! initialisation 114 ! ! Read the corrective factor on precipitations (fwfold) 115 115 CALL ctl_opn( inum, 'EMPave_old.dat', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 116 116 READ ( inum, "(24X,I8,2ES24.16)" ) iyear, a_fwb_b, a_fwb … … 138 138 zcoef = fwfold * rcp 139 139 emp(:,:) = emp(:,:) + fwfold 140 qns(:,:) = qns(:,:) - zcoef * sst_m(:,:) ! ensure fw correction does not change the heat budget140 qns(:,:) = qns(:,:) - zcoef * sst_m(:,:) ! account for change to the heat budget due to fw correction 141 141 ENDIF 142 142 ! … … 180 180 ! 181 181 emp(:,:) = emp(:,:) + zerp_cor(:,:) 182 qns(:,:) = qns(:,:) - zerp_cor(:,:) * rcp * sst_m(:,:) 182 qns(:,:) = qns(:,:) - zerp_cor(:,:) * rcp * sst_m(:,:) ! account for change to the heat budget due to fw correction 183 183 erp(:,:) = erp(:,:) + zerp_cor(:,:) 184 184 !
Note: See TracChangeset
for help on using the changeset viewer.