- Timestamp:
- 2011-06-27T13:18:25+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90
r2777 r2789 98 98 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ' 99 99 100 t_bo(:,:) = tfreez( sn(:,:,1) ) * tmask(:,:,1) ! freezing/melting point of sea water [Celcius]100 t_bo(:,:) = tfreez( tsn(:,:,1,jp_sal) ) * tmask(:,:,1) ! freezing/melting point of sea water [Celcius] 101 101 102 102 DO jj = 1, jpj ! ice if sst <= t-freez + ttest 103 103 DO ji = 1, jpi 104 IF( t n(ji,jj,1) - t_bo(ji,jj) >= ttest ) THEN ; zidto(ji,jj) = 0.e0 ! no ice105 ELSE ; zidto(ji,jj) = 1.e0 ! ice104 IF( tsn(ji,jj,1,jp_tem) - t_bo(ji,jj) >= ttest ) THEN ; zidto(ji,jj) = 0.e0 ! no ice 105 ELSE ; zidto(ji,jj) = 1.e0 ! ice 106 106 ENDIF 107 107 END DO
Note: See TracChangeset
for help on using the changeset viewer.