[24] | 1 | !> \file prop_th_icetemp.f90 |
---|
| 2 | !! Subroutines for thermal propreties calculation |
---|
| 3 | !< |
---|
| 4 | |
---|
| 5 | !> \namespace prop_th_icetemp |
---|
| 6 | !! Subroutines to calculate thermal prop |
---|
| 7 | |
---|
| 8 | !> SUBROUTINE: Thermal_prop_icetemp |
---|
| 9 | !! |
---|
| 10 | !! Calculate the thermal prop of ice |
---|
| 11 | !! |
---|
| 12 | !! Used modules: |
---|
| 13 | !! - Icetemp_declar |
---|
| 14 | !> |
---|
| 15 | |
---|
| 16 | Subroutine Thermal_prop_icetemp |
---|
[77] | 17 | !$ USE OMP_LIB |
---|
[24] | 18 | Use Icetemp_declar |
---|
| 19 | |
---|
| 20 | Implicit None |
---|
| 21 | |
---|
| 22 | If (Itracebug.Eq.1) Write(Num_tracebug,*)' Entree Dans Routine Thermal_prop_icetemp' |
---|
| 23 | |
---|
[77] | 24 | !$OMP PARALLEL |
---|
| 25 | !$OMP DO COLLAPSE(2) |
---|
[24] | 26 | Do K=1,Nz |
---|
| 27 | Do J=1,Ny |
---|
| 28 | Do I=1,Nx |
---|
| 29 | ! Temperature De Fusion:Tpmp (En Deg C) |
---|
| 30 | ! Capacite Calorifique : Cp En J/Kg/K Puis Multiplie Par Ro En J/M3/K |
---|
| 31 | ! Conductivite : Ct En J/M/K/An |
---|
| 32 | ! Attention Pour La Conductivite C'Est La Formule De Yin-Chao Yen |
---|
| 33 | |
---|
[77] | 34 | ! Ct(I,J,K)=6.62E7 |
---|
| 35 | ! Cp(I,J,K)=2009. |
---|
[24] | 36 | Tpmp(I,J,K)=-0.00087*(K-1)*1./(Nz-1)*H(I,J)!De=1/(Nz_m-1) |
---|
| 37 | Cp(I,J,K)=(2115.3+7.79293*T(I,J,K)) |
---|
| 38 | Cp(I,J,K)=Cp(I,J,K)*Ro ! En J/M3/K |
---|
| 39 | Ct(I,J,K)=3.1014e8*Exp(-0.0057*(T(I,J,K)+273.15)) !En J/M/K/An |
---|
| 40 | End Do |
---|
| 41 | End Do |
---|
| 42 | End Do |
---|
[77] | 43 | !$OMP END DO |
---|
| 44 | !$OMP END PARALLEL |
---|
[24] | 45 | End Subroutine Thermal_prop_icetemp |
---|
| 46 | |
---|
| 47 | !----------------------------------------------------------------------------------------------------------- |
---|
| 48 | |
---|
| 49 | !> SUBROUTINE: Temp_mer |
---|
| 50 | !! Calculate temperature of the sea |
---|
| 51 | !! Used modules: |
---|
| 52 | !! - Icetemp_declar |
---|
| 53 | !> |
---|
| 54 | Subroutine Temp_mer |
---|
| 55 | Use Icetemp_declar |
---|
| 56 | Implicit None |
---|
| 57 | |
---|
| 58 | If (Itracebug.Eq.1) Write(Num_tracebug,*)' Entree Dans Routine Temp_mer_icetemp ' |
---|
| 59 | |
---|
| 60 | Do J=1,Ny |
---|
| 61 | Do I=1,Nx |
---|
| 62 | If (Flot(I,J)) & |
---|
| 63 | Tbmer(I,J) = Acof1*S0mer + Bcof1 + Ccof1*H(I,J)*Ro/Row |
---|
| 64 | End Do |
---|
| 65 | End Do |
---|
| 66 | End Subroutine Temp_mer |
---|