Changeset 834 for trunk/NEMO/LIM_SRC_3/limthd.F90
- Timestamp:
- 2008-03-07T18:11:35+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/limthd.F90
r825 r834 6 6 #if defined key_lim3 7 7 !!---------------------------------------------------------------------- 8 !! 'key_lim3' : LIMsea-ice model8 !! 'key_lim3' LIM3 sea-ice model 9 9 !!---------------------------------------------------------------------- 10 10 !! lim_thd : thermodynamic of sea ice … … 22 22 USE dom_ice ! LIM sea-ice domain 23 23 USE iceini 24 USE limicepoints25 24 USE limthd_dif 26 25 USE limthd_dh … … 84 83 !! * Local variables 85 84 INTEGER :: ji, jj, jk, jl, & 86 zjid , zjjd , &87 85 zji , zjj, & ! dummy loop indices 88 86 nbpb , & ! nb of icy pts for thermo. cal. … … 95 93 REAL(wp) :: & 96 94 zinda , & ! switch for test. the val. of concen. 97 zindb, zindg ,& ! switches for test. the val of arg98 z a , zh, zthsnice, &95 zindb, & ! switches for test. the val of arg 96 zthsnice , & 99 97 zfric_u , & ! friction velocity 100 98 zfnsol , & ! total non solar heat … … 104 102 105 103 REAL(wp) :: & 106 ztmelts , & ! ice layers melting point107 zaaa , & ! numeric value for 2nd degree polynomial108 zbbb , & ! " " " " " "109 zccc , & ! " " " " " "110 zdiscrim , & ! " " " " " "111 zq , & ! sea ice temporary energy of melting112 zqold , & ! test variable for ice heat content113 104 zareamin 114 105 … … 224 215 225 216 ! still need to be updated : fdtcn !!!! 226 ! !-- Lead heat budget (part 1, next one is in limthd_dh , then limthd_lab217 ! !-- Lead heat budget (part 1, next one is in limthd_dh 227 218 ! !-- qldif -- (or qldif_1d in 1d routines) 228 219 zfontn = sprecip(ji,jj) * lfus ! energy of melting … … 510 501 REAL(wp) :: & 511 502 zdes, & ! snow heat content increment (dummy) 512 ztmelts, & ! layer melting point (dummy)513 503 zeps ! very small value (1.e-10) 514 504 … … 560 550 jl !: category number 561 551 562 REAL(wp), DIMENSION(jpl) :: & !: ! goes to trash563 zdummy564 565 552 REAL(wp) :: & !: ! goes to trash 566 553 meance, & !: mean conservation error … … 573 560 INTEGER :: & 574 561 ji,jj,jk, & !: loop indices 575 zji, zjj , zjiindex562 zji, zjj 576 563 !!--------------------------------------------------------------------- 577 564 … … 745 732 jl !: category number 746 733 747 REAL(wp), DIMENSION(jpl) :: & !: ! goes to trash748 zdummy749 750 734 REAL(wp) :: & !: ! goes to trash 751 735 meance, & !: mean conservation error 752 max_cons_err, & !: maximum tolerated conservation error 753 max_surf_err !: maximum tolerated surface error 736 max_cons_err !: maximum tolerated conservation error 754 737 755 738 INTEGER :: & … … 758 741 INTEGER :: & 759 742 ji,jj,jk, & !: loop indices 760 zji, zjj , zjiindex743 zji, zjj 761 744 762 745 !!--------------------------------------------------------------------- … … 974 957 ENDIF 975 958 976 uscomi = 1.0 / ( 1.0 - amax ) ! inverse of minimum lead fraction977 959 rcdsn = hakdif * rcdsn 978 960 rcdic = hakdif * rcdic
Note: See TracChangeset
for help on using the changeset viewer.