Changeset 3266 for trunk/LMDZ.GENERIC
- Timestamp:
- Mar 13, 2024, 5:13:20 PM (9 months ago)
- Location:
- trunk/LMDZ.GENERIC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/changelog.txt
r3265 r3266 1883 1883 == 13/03/2024 == MT 1884 1884 1885 Very minor change of max sea ice fraction in ocean model 1885 - Very minor change of max sea ice fraction in ocean model 1886 - Preliminary cleanup+commenting in hydrol.F90 -
trunk/LMDZ.GENERIC/libf/phystd/hydrol.F90
r3194 r3266 55 55 real snowlayer 56 56 parameter (snowlayer=33.0) ! 33 kg/m^2 of snow, equal to a layer of 3.3 cm 57 real oceantime 57 58 real oceantime ! this is a relaxation timescale for the oceanbulkavg parameterization 58 59 parameter (oceantime=10*24*3600) 59 60 logical,save :: oceanbulkavg ! relax ocean temperatures to a GLOBAL mean value?60 logical,save :: oceanbulkavg ! simple parameterization to relax ocean temperatures to the global mean value (crude, 0th order parameterization to mimick ocean heat transport) 61 61 62 logical,save :: activerunoff ! enable simple runoff scheme? 62 logical,save :: oceanalbvary ! ocean albedo varies with the diurnal cycle?63 logical,save :: oceanalbvary ! simple parameterization to account for the effect of solar zenith angle on the ocean albedo (for the moment it is not used, but to be included in the future) 63 64 !$OMP THREADPRIVATE(oceanbulkavg,activerunoff,oceanalbvary) 64 65 65 66 ! Arguments 66 67 ! --------- 67 real rnat(ngrid) ! I changed this to integer (RW)68 real rnat(ngrid) ! rnat is terrain type: 0-ocean; 1-continent 68 69 real,dimension(:),allocatable,save :: runoff 69 70 real totalrunoff, tsea, oceanarea … … 175 176 endif 176 177 177 ! call writediagfi(ngrid,"hydrol_rnat"," "," ",2,rnat) 178 !! write (*,*) "oceanarea", oceanarea 178 ! write (*,*) "oceanarea", oceanarea 179 179 180 180 ! add physical tendencies already calculated … … 188 188 enddo 189 189 enddo 190 191 ! call writediagfi(ngrid,"hydrol_zqsurf_iliq_1"," "," ",2,zqsurf(1:ngrid,iliq))192 190 193 191 do ig=1,ngrid … … 350 348 end do ! ig=1,ngrid 351 349 352 ! call writediagfi(ngrid,"hydrol_zqsurf_iliq_2"," "," ",2,zqsurf(1:ngrid,iliq)) 353 354 ! perform crude bulk averaging of temperature in ocean 350 351 ! simple parameterization to perform crude bulk averaging of temperature in ocean 355 352 ! ---------------------------------------------------- 356 if(oceanbulkavg)then 353 if(oceanbulkavg)then 357 354 358 355 oceanarea2=0. … … 433 430 endif !activerunoff 434 431 435 ! call writediagfi(ngrid,"hydrol_zqsurf_iliq_3"," "," ",2,zqsurf(1:ngrid,iliq))436 437 432 ! Re-add the albedo effects of CO2 ice if necessary 438 433 ! ------------------------------------------------- … … 455 450 enddo 456 451 457 ! call writediagfi(ngrid,"hydrol_dqs_hyd_iliq"," "," ",2,dqs_hyd(1:ngrid,iliq))458 459 452 if (activerunoff) then 460 453 call writediagfi(ngrid,'runoff','Runoff amount',' ',2,runoff)
Note: See TracChangeset
for help on using the changeset viewer.