- Timestamp:
- Feb 22, 2021, 12:44:07 PM (4 years ago)
- Location:
- LMDZ6/branches/LMDZ-tracers
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ-tracers
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ-tracers/libf/phylmd/pbl_surface_mod.F90
r3435 r3851 23 23 USE climb_wind_mod, ONLY : climb_wind_down, climb_wind_up 24 24 USE coef_diff_turb_mod, ONLY : coef_diff_turb 25 USE ioipsl_getin_p_mod, ONLY : getin_p 26 USE cdrag_mod 27 USE stdlevvar_mod 25 28 USE wx_pbl_mod, ONLY : wx_pbl_init, wx_pbl_final, & 26 29 !! wx_pbl_fuse_no_dts, wx_pbl_split_no_dts, & 27 30 !! wx_pbl_fuse, wx_pbl_split 28 31 wx_pbl0_fuse, wx_pbl0_split 32 use config_ocean_skin_m, only: activate_ocean_skin 29 33 30 34 IMPLICIT NONE … … 42 46 INTEGER, SAVE :: iflag_pbl_surface_t2m_bug 43 47 !$OMP THREADPRIVATE(iflag_pbl_surface_t2m_bug) 48 integer, save :: iflag_new_t2mq2m 49 !$OMP THREADPRIVATE(iflag_new_t2mq2m) 50 44 51 !FC 45 52 ! integer, save :: iflag_frein … … 164 171 rlon, rlat, rugoro, rmu0, & 165 172 zsig, lwdown_m, pphi, cldt, & 166 rain_f, snow_f, solsw_m, sol lw_m, &173 rain_f, snow_f, solsw_m, solswfdiff_m, sollw_m, & 167 174 gustiness, & 168 175 t, q, u, v, & … … 177 184 alb_dir_m, alb_dif_m, zxsens, zxevap, & 178 185 alb3_lic, runoff, snowhgt, qsnow, to_ice, sissnow, & 179 zxtsol, zxfluxlat, zt2m, qsat2m, 186 zxtsol, zxfluxlat, zt2m, qsat2m, zn2mout, & 180 187 d_t, d_q, d_u, d_v, d_t_diss, & 181 188 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012 … … 206 213 !jyg< 207 214 !! zxfluxt, zxfluxq, q2m, flux_q, tke, & 208 zxfluxt, zxfluxq, q2m, flux_q, tke_x, 215 zxfluxt, zxfluxq, q2m, flux_q, tke_x, & 209 216 !>jyg 210 217 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012 … … 286 293 USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, carbon_cycle_tr, level_coupling_esm 287 294 USE carbon_cycle_mod, ONLY : co2_send, nbcf_out, fields_out, yfields_out, cfname_out 295 use hbtm_mod, only: hbtm 288 296 USE indice_sol_mod 289 297 USE time_phylmdz_mod, ONLY : day_ini,annee_ref,itau_phy … … 291 299 USE print_control_mod, ONLY : prt_level,lunout 292 300 USE ioipsl_getin_p_mod, ONLY : getin_p 301 use phys_state_var_mod, only: ds_ns, dt_ns, delta_sst, delta_sal 302 use phys_output_var_mod, only: dter, dser, tkt, tks, taur, sss 303 #ifdef CPP_XIOS 304 USE wxios, ONLY: missing_val 305 #else 306 use netcdf, only: missing_val => nf90_fill_real 307 #endif 293 308 294 309 IMPLICIT NONE … … 318 333 REAL, DIMENSION(klon), INTENT(IN) :: snow_f ! snow fall 319 334 REAL, DIMENSION(klon), INTENT(IN) :: solsw_m ! net shortwave radiation at mean surface 335 REAL, DIMENSION(klon), INTENT(IN) :: solswfdiff_m ! diffuse fraction fordownward shortwave radiation at mean surface 320 336 REAL, DIMENSION(klon), INTENT(IN) :: sollw_m ! net longwave radiation at mean surface 321 337 REAL, DIMENSION(klon,klev), INTENT(IN) :: t ! temperature (K) … … 394 410 REAL, DIMENSION(klon), INTENT(OUT) :: zxfluxlat ! latent flux, mean for each grid point 395 411 REAL, DIMENSION(klon), INTENT(OUT) :: zt2m ! temperature at 2m, mean for each grid point 412 INTEGER, DIMENSION(klon, 6), INTENT(OUT) :: zn2mout ! number of times the 2m temperature is out of the [tsol,temp] 396 413 REAL, DIMENSION(klon), INTENT(OUT) :: qsat2m 397 414 REAL, DIMENSION(klon, klev), INTENT(OUT) :: d_t ! change in temperature … … 453 470 REAL, DIMENSION(klon, klev), INTENT(OUT) :: zxfluxu ! u wind tension, mean for each grid point 454 471 REAL, DIMENSION(klon, klev), INTENT(OUT) :: zxfluxv ! v wind tension, mean for each grid point 455 REAL, DIMENSION(klon, nbsrf+1), INTENT(INOUT) 456 REAL, DIMENSION(klon, nbsrf), INTENT(INOUT) 472 REAL, DIMENSION(klon, nbsrf+1), INTENT(INOUT) :: z0m,z0h ! rugosity length (m) 473 REAL, DIMENSION(klon, nbsrf), INTENT(INOUT) :: agesno ! age of snow at surface 457 474 REAL, DIMENSION(klon, nbsrf), INTENT(OUT) :: solsw ! net shortwave radiation at surface 458 475 REAL, DIMENSION(klon, nbsrf), INTENT(OUT) :: sollw ! net longwave radiation at surface 459 476 REAL, DIMENSION(klon, nbsrf), INTENT(OUT) :: d_ts ! change in temperature at surface 460 REAL, DIMENSION(klon, nbsrf), INTENT(INOUT) 477 REAL, DIMENSION(klon, nbsrf), INTENT(INOUT) :: evap ! evaporation at surface 461 478 REAL, DIMENSION(klon, nbsrf), INTENT(OUT) :: fluxlat ! latent flux 462 479 REAL, DIMENSION(klon, nbsrf), INTENT(OUT) :: t2m ! temperature at 2 meter height … … 537 554 REAL, DIMENSION(klon) :: y_flux_u1, y_flux_v1 538 555 REAL, DIMENSION(klon) :: yt2m, yq2m, yu10m 556 INTEGER, DIMENSION(klon, nbsrf, 6) :: yn2mout, yn2mout_x, yn2mout_w 557 INTEGER, DIMENSION(klon, nbsrf, 6) :: n2mout, n2mout_x, n2mout_w 539 558 REAL, DIMENSION(klon) :: yustar 540 559 REAL, DIMENSION(klon) :: ywstar … … 822 841 ! Martin 823 842 843 real, DIMENSION(klon):: ydelta_sst, ydelta_sal, yds_ns, ydt_ns, ydter, ydser, & 844 ytkt, ytks, ytaur, ysss 845 ! compression of delta_sst, delta_sal, ds_ns, dt_ns, dter, dser, tkt, tks, 846 ! taur, sss on ocean points 847 824 848 !**************************************************************************************** 825 849 ! End of declarations … … 838 862 839 863 IF (first_call) THEN 864 865 iflag_new_t2mq2m=1 866 CALL getin_p('iflag_new_t2mq2m',iflag_new_t2mq2m) 867 print*,'pbl_iflag_new_t2mq2m=',iflag_new_t2mq2m 868 840 869 print*,'PBL SURFACE AVEC GUSTINESS' 841 870 first_call=.FALSE. … … 843 872 ! Initialize ok_flux_surf (for 1D model) 844 873 if (klon_glo>1) ok_flux_surf=.FALSE. 874 if (klon_glo>1) ok_forc_tsurf=.FALSE. 845 875 846 876 ! intialize beta_land … … 903 933 zxfluxlat(:)=0. 904 934 zt2m(:)=0. ; zq2m(:)=0. ; qsat2m(:)=0. ; rh2m(:)=0. 935 zn2mout(:,:)=0 ; 905 936 d_t(:,:)=0. ; d_t_diss(:,:)=0. ; d_q(:,:)=0. ; d_u(:,:)=0. ; d_v(:,:)=0. 906 937 zcoefh(:,:,:)=0. ; zcoefm(:,:,:)=0. … … 945 976 !! tke(:,:,is_ave)=0. 946 977 tke_x(:,:,is_ave)=0. 978 947 979 wake_dltke(:,:,is_ave)=0. 948 980 !>jyg … … 955 987 !!! jyg le 10/02/2012 956 988 rh2m_x(:) = 0. ; qsat2m_x(:) = 0. ; rh2m_w(:) = 0. ; qsat2m_w(:) = 0. 957 !!!958 989 959 990 ! 2b) Initialization of all local variables that will be compressed later … … 965 996 yqsurf = 0.0 ; yalb = 0.0 ; yalb_vis = 0.0 966 997 !albedo SB <<< 967 yrain_f = 0.0 ; ysnow_f = 0.0 ; yfder = 0.0 ; ysolsw = 0.0 998 yrain_f = 0.0 ; ysnow_f = 0.0 ; yfder = 0.0 ; ysolsw = 0.0 968 999 ysollw = 0.0 ; yz0m = 0.0 ; yz0h = 0.0 ; yu1 = 0.0 969 1000 yv1 = 0.0 ; ypaprs = 0.0 ; ypplay = 0.0 … … 976 1007 !! d_t_diss= 0.0 ;d_u = 0.0 ; d_v = 0.0 977 1008 yqsol = 0.0 978 ytherm = 0.0 ; ytke=0. 1009 1010 ytke=0. 979 1011 !FC 980 1012 y_treedrg=0. … … 1183 1215 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1184 1216 1217 !--OB this line is not satisfactory because alb is the direct albedo not total albedo 1185 1218 solsw(i,nsrf) = solsw_m(i) * (1.-alb(i,nsrf)) / (1.-alb_m(i)) 1186 1219 ENDDO … … 1205 1238 !>al1 1206 1239 1240 !--OB add diffuse fraction of SW down 1241 DO n=1,nbcf_out 1242 IF (cfname_out(n) == "swdownfdiff" ) fields_out(:,n) = solswfdiff_m(:) 1243 ENDDO 1207 1244 ! >> PC 1208 1245 IF (carbon_cycle_cpl .AND. carbon_cycle_tr .AND. nbcf_out.GT.0 ) THEN … … 1381 1418 ytke_w(j,k) = tke_x(i,k,nsrf)+wake_dltke(i,k,nsrf) 1382 1419 ywake_dltke(j,k) = wake_dltke(i,k,nsrf) 1420 1383 1421 !>jyg 1384 1422 ENDDO … … 1415 1453 ENDDO 1416 1454 ENDIF 1455 1456 if (nsrf == is_oce .and. activate_ocean_skin >= 1) then 1457 if (activate_ocean_skin == 2 .and. type_ocean == "couple") then 1458 ydelta_sal(:knon) = delta_sal(ni(:knon)) 1459 ydelta_sst(:knon) = delta_sst(ni(:knon)) 1460 end if 1461 1462 yds_ns(:knon) = ds_ns(ni(:knon)) 1463 ydt_ns(:knon) = dt_ns(ni(:knon)) 1464 end if 1417 1465 1418 1466 !**************************************************************************************** … … 1455 1503 ENDDO 1456 1504 ENDIF 1505 1457 1506 IF (prt_level >=10) print *,'clcdrag -> ycdragh ', ycdragh 1458 1507 ELSE !(iflag_split .eq.0) … … 1538 1587 print *,' args coef_diff_turb: ycdragh ', ycdragh 1539 1588 print *,' args coef_diff_turb: ytke ', ytke 1589 1540 1590 ENDIF 1541 1591 CALL coef_diff_turb(dtime, nsrf, knon, ni, & … … 1567 1617 print *,' args coef_diff_turb: ycdragh_x ', ycdragh_x 1568 1618 print *,' args coef_diff_turb: ytke_x ', ytke_x 1619 1569 1620 ENDIF 1570 1621 CALL coef_diff_turb(dtime, nsrf, knon, ni, & … … 1821 1872 1822 1873 ! Calculate the temperature et relative humidity at 2m and the wind at 10m 1874 IF (iflag_new_t2mq2m==1) THEN 1875 CALL stdlevvarn(klon, knon, is_ter, zxli, & 1876 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 1877 yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), & 1878 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, & 1879 yn2mout(:, nsrf, :)) 1880 ELSE 1823 1881 CALL stdlevvar(klon, knon, is_ter, zxli, & 1824 1882 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 1825 1883 yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), & 1826 1884 yt2m, yq2m, yt10m, yq10m, yu10m, yustar) 1885 ENDIF 1827 1886 1828 1887 ENDIF … … 1935 1994 yz0m, yz0h, SFRWL,yalb_dir_new, yalb_dif_new, yevap, yfluxsens,yfluxlat,& 1936 1995 ytsurf_new, y_dflux_t, y_dflux_q, slab_wfbils, & 1937 y_flux_u1, y_flux_v1) 1996 y_flux_u1, y_flux_v1, ydelta_sst(:knon), ydelta_sal(:knon), & 1997 yds_ns(:knon), ydt_ns(:knon), ydter(:knon), ydser(:knon), & 1998 ytkt(:knon), ytks(:knon), ytaur(:knon), ysss) 1938 1999 IF (prt_level >=10) THEN 1939 2000 print *,'arg de surf_ocean: ycdragh ',ycdragh … … 2013 2074 ! 2014 2075 !**************************************************************************************** 2015 2016 !!! 2017 !!! jyg le 10/04/2013 2076 !! 2077 !!! 2078 !!! jyg le 10/04/2013 et EV 10/2020 2079 2080 IF (ok_forc_tsurf) THEN 2081 DO j=1,knon 2082 ytsurf_new(j)=tg 2083 y_d_ts(j) = ytsurf_new(j) - yts(j) 2084 ENDDO 2085 ENDIF ! ok_forc_tsurf 2086 2018 2087 !!! 2019 2088 IF (ok_flux_surf) THEN … … 2444 2513 tke_x(i,k,nsrf) = ytke(j,k) 2445 2514 tke_x(i,k,is_ave) = tke_x(i,k,is_ave) + ytke(j,k)*ypct(j) 2515 2446 2516 !>jyg 2447 2517 ENDDO … … 2457 2527 !! tke(i,k,is_ave) = tke(i,k,is_ave) + tke(i,k,nsrf)*ypct(j) 2458 2528 tke_x(i,k,nsrf) = ytke_x(j,k) 2459 tke_x(i,k,is_ave) = tke_x(i,k,is_ave) + tke_x(i,k,nsrf)*ypct(j) 2529 tke_x(i,k,is_ave) = tke_x(i,k,is_ave) + tke_x(i,k,nsrf)*ypct(j) 2460 2530 wake_dltke(i,k,is_ave) = wake_dltke(i,k,is_ave) + wake_dltke(i,k,nsrf)*ypct(j) 2531 2461 2532 2462 2533 !>jyg … … 2538 2609 d_t_w(:,1), d_t_x(:,1), d_t(:,1) 2539 2610 ENDIF 2611 2612 if (nsrf == is_oce .and. activate_ocean_skin >= 1) then 2613 delta_sal = missing_val 2614 ds_ns = missing_val 2615 dt_ns = missing_val 2616 delta_sst = missing_val 2617 dter = missing_val 2618 dser = missing_val 2619 tkt = missing_val 2620 tks = missing_val 2621 taur = missing_val 2622 sss = missing_val 2623 2624 delta_sal(ni(:knon)) = ydelta_sal(:knon) 2625 ds_ns(ni(:knon)) = yds_ns(:knon) 2626 dt_ns(ni(:knon)) = ydt_ns(:knon) 2627 delta_sst(ni(:knon)) = ydelta_sst(:knon) 2628 dter(ni(:knon)) = ydter(:knon) 2629 dser(ni(:knon)) = ydser(:knon) 2630 tkt(ni(:knon)) = ytkt(:knon) 2631 tks(ni(:knon)) = ytks(:knon) 2632 taur(ni(:knon)) = ytaur(:knon) 2633 sss(ni(:knon)) = ysss(:knon) 2634 end if 2540 2635 2541 2636 !**************************************************************************************** … … 2614 2709 !!! jyg le 07/02/2012 2615 2710 IF (iflag_split .eq.0) THEN 2711 IF (iflag_new_t2mq2m==1) THEN 2712 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2713 uzon, vmer, tair1, qair1, zgeo1, & 2714 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2715 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, & 2716 yn2mout(:, nsrf, :)) 2717 ELSE 2616 2718 CALL stdlevvar(klon, knon, nsrf, zxli, & 2617 2719 uzon, vmer, tair1, qair1, zgeo1, & 2618 2720 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2619 2721 yt2m, yq2m, yt10m, yq10m, yu10m, yustar) 2722 ENDIF 2620 2723 ELSE !(iflag_split .eq.0) 2724 IF (iflag_new_t2mq2m==1) THEN 2725 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2726 uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, & 2727 tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2728 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x, & 2729 yn2mout_x(:, nsrf, :)) 2730 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2731 uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, & 2732 tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2733 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w, & 2734 yn2mout_w(:, nsrf, :)) 2735 ELSE 2621 2736 CALL stdlevvar(klon, knon, nsrf, zxli, & 2622 2737 uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, & … … 2627 2742 tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2628 2743 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w) 2744 ENDIF 2629 2745 !!! 2630 2746 ENDIF ! (iflag_split .eq.0) … … 2640 2756 u10m(i,nsrf)=(yu10m(j) * uzon(j))/SQRT(uzon(j)**2+vmer(j)**2) 2641 2757 v10m(i,nsrf)=(yu10m(j) * vmer(j))/SQRT(uzon(j)**2+vmer(j)**2) 2758 ! 2759 DO k = 1, 6 2760 n2mout(i,nsrf,k) = yn2mout(j,nsrf,k) 2761 END DO 2762 ! 2642 2763 ENDDO 2643 2764 ELSE !(iflag_split .eq.0) … … 2650 2771 u10m_x(i,nsrf)=(yu10m_x(j) * uzon_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2) 2651 2772 v10m_x(i,nsrf)=(yu10m_x(j) * vmer_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2) 2773 ! 2774 DO k = 1, 6 2775 n2mout_x(i,nsrf,k) = yn2mout_x(j,nsrf,k) 2776 END DO 2777 ! 2652 2778 ENDDO 2653 2779 DO j=1, knon … … 2663 2789 u10m(i,nsrf) = u10m_x(i,nsrf) + wake_s(i)*(u10m_w(i,nsrf)-u10m_x(i,nsrf)) 2664 2790 v10m(i,nsrf) = v10m_x(i,nsrf) + wake_s(i)*(v10m_w(i,nsrf)-v10m_x(i,nsrf)) 2791 ! 2792 DO k = 1, 6 2793 n2mout_w(i,nsrf,k) = yn2mout_w(j,nsrf,k) 2794 END DO 2795 ! 2665 2796 ENDDO 2666 2797 !!! … … 2919 3050 ! 2920 3051 zxtsol(:) = 0.0 ; zxfluxlat(:) = 0.0 2921 zt2m(:) = 0.0 ; zq2m(:) = 0.0 3052 zt2m(:) = 0.0 ; zq2m(:) = 0.0 ; zn2mout(:,:) = 0 2922 3053 zustar(:)=0.0 ; zu10m(:) = 0.0 ; zv10m(:) = 0.0 2923 3054 s_pblh(:) = 0.0 ; s_plcl(:) = 0.0 … … 2972 3103 zt2m(i) = zt2m(i) + t2m(i,nsrf) * pctsrf(i,nsrf) 2973 3104 zq2m(i) = zq2m(i) + q2m(i,nsrf) * pctsrf(i,nsrf) 3105 ! 3106 DO k = 1, 6 3107 zn2mout(i,k) = zn2mout(i,k) + n2mout(i,nsrf,k) * pctsrf(i,nsrf) 3108 ENDDO 3109 ! 2974 3110 zustar(i) = zustar(i) + ustar(i,nsrf) * pctsrf(i,nsrf) 2975 3111 wstar(i,is_ave)=wstar(i,is_ave)+wstar(i,nsrf)*pctsrf(i,nsrf) … … 3003 3139 zt2m(i) = zt2m(i) + (t2m_x(i,nsrf)+wake_s(i)*(t2m_w(i,nsrf)-t2m_x(i,nsrf))) * pctsrf(i,nsrf) 3004 3140 zq2m(i) = zq2m(i) + q2m_x(i,nsrf) * pctsrf(i,nsrf) 3141 ! 3142 DO k = 1, 6 3143 zn2mout(i,k) = zn2mout(i,k) + n2mout_x(i,nsrf,k) * pctsrf(i,nsrf) 3144 ENDDO 3145 ! 3005 3146 zustar(i) = zustar(i) + ustar_x(i,nsrf) * pctsrf(i,nsrf) 3006 3147 wstar(i,is_ave)=wstar(i,is_ave)+wstar_x(i,nsrf)*pctsrf(i,nsrf) … … 3081 3222 DO nsrf = 1, nbsrf 3082 3223 DO i = 1, klon 3083 zxqsurf(i) = zxqsurf(i) + qsurf(i,nsrf) * pctsrf(i,nsrf)3224 zxqsurf(i) = zxqsurf(i) + MAX(qsurf(i,nsrf),0.0) * pctsrf(i,nsrf) 3084 3225 zxsnow(i) = zxsnow(i) + snow(i,nsrf) * pctsrf(i,nsrf) 3085 3226 ENDDO … … 3142 3283 3143 3284 !albedo SB >>> 3144 SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, &3145 evap, z0m, z0h, agesno, &3146 tsurf,alb_dir,alb_dif, ustar, u10m, v10m, tke)3147 !albedo SB <<<3285 SUBROUTINE pbl_surface_newfrac(itime, pctsrf_new, pctsrf_old, & 3286 evap, z0m, z0h, agesno, & 3287 tsurf,alb_dir,alb_dif, ustar, u10m, v10m, tke) 3288 !albedo SB <<< 3148 3289 ! Give default values where new fraction has appread 3149 3290 3150 3291 USE indice_sol_mod 3292 use phys_state_var_mod, only: delta_sal, ds_ns, dt_ns, delta_sst 3293 use config_ocean_skin_m, only: activate_ocean_skin 3151 3294 3152 3295 INCLUDE "dimsoil.h" … … 3235 3378 alb_dif(i,k,nsrf) = 0.06 3236 3379 ENDDO 3380 if (activate_ocean_skin >= 1) then 3381 if (activate_ocean_skin == 2 & 3382 .and. type_ocean == "couple") then 3383 delta_sal(i) = 0. 3384 delta_sst(i) = 0. 3385 end if 3386 3387 ds_ns(i) = 0. 3388 dt_ns(i) = 0. 3389 end if 3237 3390 ELSE IF (nsrf.EQ.is_sic) THEN 3238 3391 tsurf(i,nsrf) = 271.35
Note: See TracChangeset
for help on using the changeset viewer.