<p><b>laura@ucar.edu</b> 2010-12-21 15:46:06 -0700 (Tue, 21 Dec 2010)</p><p>physics_aquaplanet provides surface boundary conditions for initial testing of surface and PBL parameterizations<br>
</p><hr noshade><pre><font color="gray">Added: branches/atmos_physics/src/core_physics/module_physics_aquaplanet.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_aquaplanet.F                                (rev 0)
+++ branches/atmos_physics/src/core_physics/module_physics_aquaplanet.F        2010-12-21 22:46:06 UTC (rev 654)
@@ -0,0 +1,79 @@
+!=============================================================================================
+ module module_physics_aquaplanet
+ use configure
+ use grid_types
+!use constants
+ 
+ use module_physics_constants
+ use module_physics_vars
+
+ implicit none
+ private
+ public:: physics_aquaplanet_init
+
+
+!These variables will need to be redefined when we have a time manager.
+ integer,parameter,public:: julday = 80!Julian day                               (-).
+ real(kind=RKIND),parameter,public:: gmt    = 0 !Greenwich mean time hour of model start (hr).
+
+ contains
+
+!=============================================================================================
+ subroutine physics_aquaplanet_init(mesh,diag_physics)
+!=============================================================================================
+
+!input and inout arguments:
+!--------------------------
+ type(mesh_type),intent(in):: mesh
+ type(diag_physics_type),intent(in):: diag_physics
+
+!local variables:
+ integer,parameter:: t00_c = 27.
+ integer:: iCell,nCells
+ 
+ real(kind=RKIND),dimension(:),pointer:: latitude
+ real(kind=RKIND),dimension(:),pointer:: longitude
+ real(kind=RKIND),dimension(:),pointer:: xland
+ real(kind=RKIND),dimension(:),pointer:: xice
+ real(kind=RKIND),dimension(:),pointer:: sfc_albedo
+ real(kind=RKIND),dimension(:),pointer:: sfc_temperature
+
+!---------------------------------------------------------------------------------------------
+
+ write(0,*)
+ write(0,*) '--- enter subroutine physics_aquaplanet:'
+
+ nCells = mesh % nCells
+
+ latitude  =&gt; mesh % latCell % array
+ longitude =&gt; mesh % lonCell % array
+ xland           =&gt; diag_physics % xland % array
+ xice            =&gt; diag_physics % xice % array
+ sfc_albedo      =&gt; diag_physics % sfc_albedo % array
+ sfc_temperature =&gt; diag_physics % tsk % array
+
+!set surface conditions to all oceans:
+
+ do iCell = 1, nCells
+    xice(iCell)  = 0.  !no ice
+    xland(iCell) = 2.  !all water
+ enddo
+
+!compute fixed sea-surface temperatures:
+ do iCell = 1, nCells
+    
+    if(latitude(iCell) .gt. -pii/3 .and. latitude(iCell) .lt. pii/3) &amp;
+       sfc_temperature(iCell) = t00_c*(1.-sin(3*latitude(iCell)/2)**2) + t00
+
+ enddo
+
+!set surface albedo:
+ do iCell = 1, nCells
+    sfc_albedo(iCell) = 0.03
+ enddo
+
+ end subroutine physics_aquaplanet_init
+
+!=============================================================================================
+ end module module_physics_aquaplanet
+!=============================================================================================

</font>
</pre>