Class | snowice_frac |
In: |
surface_properties/snowice_frac.f90
|
Subroutine : | |
xy_SurfMajCompIce( 0:imax-1, 1:jmax ) : | real(DP), intent(in ) |
xy_MajCompIceFrac( 0:imax-1, 1:jmax ) : | real(DP), intent(out) |
subroutine CalcMajCompIceFrac( xy_SurfMajCompIce, xy_MajCompIceFrac ) ! 格子点設定 ! Grid points settings ! use gridset, only: imax, jmax, kmax ! 鉛直層数. ! Number of vertical level ! 雪と海氷の定数の設定 ! Setting constants of snow and sea ice ! use constants_snowseaice, only: CO2IceThreshold !!$ logical , intent(in ) :: xy_FlagLand ( 0:imax-1, 1:jmax ) real(DP), intent(in ) :: xy_SurfMajCompIce( 0:imax-1, 1:jmax ) real(DP), intent(out) :: xy_MajCompIceFrac( 0:imax-1, 1:jmax ) ! 作業変数 ! Work variables ! real(DP):: MajCompIceThreshold integer:: i ! 経度方向に回る DO ループ用作業変数 ! Work variables for DO loop in longitude integer:: j ! 緯度方向に回る DO ループ用作業変数 ! Work variables for DO loop in latitude ! 初期化確認 ! Initialization check ! if ( .not. snowice_frac_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if MajCompIceThreshold = CO2IceThreshold do j = 1, jmax do i = 0, imax-1 !!$ if ( xy_FlagLand(i,j) ) then xy_MajCompIceFrac(i,j) = ( xy_SurfMajCompIce(i,j) - 0.0_DP ) / ( MajCompIceThreshold - 0.0_DP ) xy_MajCompIceFrac(i,j) = max( min( xy_MajCompIceFrac(i,j), 1.0_DP ), 0.0_DP ) !!$ else !!$ xy_MajCompIceFrac(i,j) = 0.0_DP !!$ end if end do end do end subroutine CalcMajCompIceFrac
Subroutine : | |
xy_SeaIceConc( 0:imax-1, 1:jmax ) : | real(DP), intent(in ) |
xy_SeaIceFrac( 0:imax-1, 1:jmax ) : | real(DP), intent(out) |
subroutine CalcSeaIceFrac( xy_SeaIceConc, xy_SeaIceFrac ) ! 格子点設定 ! Grid points settings ! use gridset, only: imax, jmax, kmax ! 鉛直層数. ! Number of vertical level ! 雪と海氷の定数の設定 ! Setting constants of snow and sea ice ! use constants_snowseaice, only: SeaIceThreshold !!$ logical , intent(in ) :: xy_FlagOcean ( 0:imax-1, 1:jmax ) real(DP), intent(in ) :: xy_SeaIceConc( 0:imax-1, 1:jmax ) real(DP), intent(out) :: xy_SeaIceFrac( 0:imax-1, 1:jmax ) ! 作業変数 ! Work variables ! integer:: i ! 経度方向に回る DO ループ用作業変数 ! Work variables for DO loop in longitude integer:: j ! 緯度方向に回る DO ループ用作業変数 ! Work variables for DO loop in latitude ! 初期化確認 ! Initialization check ! if ( .not. snowice_frac_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if do j = 1, jmax do i = 0, imax-1 !!$ if ( ( xy_FlagOcean(i,j) ) .and. & !!$ & ( xy_SeaIceConc(i,j) > SeaIceThreshold ) ) then if ( xy_SeaIceConc(i,j) > SeaIceThreshold ) then xy_SeaIceFrac(i,j) = 1.0_DP else xy_SeaIceFrac(i,j) = 0.0_DP end if end do end do end subroutine CalcSeaIceFrac
Subroutine : | |
xy_SurfSnow( 0:imax-1, 1:jmax ) : | real(DP), intent(in ) |
xy_SnowFrac( 0:imax-1, 1:jmax ) : | real(DP), intent(out) |
subroutine CalcSfcLiqSolFrac( xy_SurfSnow, xy_SnowFrac ) ! 格子点設定 ! Grid points settings ! use gridset, only: imax, jmax, kmax ! 鉛直層数. ! Number of vertical level ! 雪と海氷の定数の設定 ! Setting constants of snow and sea ice ! use constants_snowseaice, only: SnowThresholdForFlux !!$ logical , intent(in ) :: xy_FlagLand( 0:imax-1, 1:jmax ) real(DP), intent(in ) :: xy_SurfSnow( 0:imax-1, 1:jmax ) real(DP), intent(out) :: xy_SnowFrac( 0:imax-1, 1:jmax ) ! 作業変数 ! Work variables ! integer:: i ! 経度方向に回る DO ループ用作業変数 ! Work variables for DO loop in longitude integer:: j ! 緯度方向に回る DO ループ用作業変数 ! Work variables for DO loop in latitude ! 初期化確認 ! Initialization check ! if ( .not. snowice_frac_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if if ( SnowThresholdForFlux > 0.0_DP ) then do j = 1, jmax do i = 0, imax-1 !!$ if ( xy_FlagLand(i,j) ) then xy_SnowFrac(i,j) = ( ( xy_SurfSnow(i,j) - 0.0_DP ) / ( SnowThresholdForFlux - 0.0_DP ) )**NOrd xy_SnowFrac(i,j) = max( min( xy_SnowFrac(i,j), 1.0_DP ), 0.0_DP ) !!$ else !!$ xy_SnowFrac(i,j) = 0.0_DP !!$ end if end do end do else xy_SnowFrac = 1.0_DP end if end subroutine CalcSfcLiqSolFrac
Subroutine : | |
xy_SurfSnow( 0:imax-1, 1:jmax ) : | real(DP), intent(in ) |
xy_SnowFrac( 0:imax-1, 1:jmax ) : | real(DP), intent(out) |
subroutine CalcSnowFrac( xy_SurfSnow, xy_SnowFrac ) ! 格子点設定 ! Grid points settings ! use gridset, only: imax, jmax, kmax ! 鉛直層数. ! Number of vertical level ! 雪と海氷の定数の設定 ! Setting constants of snow and sea ice ! use constants_snowseaice, only: SnowThresholdForAlbedo !!$ logical , intent(in ) :: xy_FlagLand( 0:imax-1, 1:jmax ) real(DP), intent(in ) :: xy_SurfSnow( 0:imax-1, 1:jmax ) real(DP), intent(out) :: xy_SnowFrac( 0:imax-1, 1:jmax ) ! 作業変数 ! Work variables ! integer:: i ! 経度方向に回る DO ループ用作業変数 ! Work variables for DO loop in longitude integer:: j ! 緯度方向に回る DO ループ用作業変数 ! Work variables for DO loop in latitude ! 初期化確認 ! Initialization check ! if ( .not. snowice_frac_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if if ( SnowThresholdForAlbedo > 0.0_DP ) then do j = 1, jmax do i = 0, imax-1 xy_SnowFrac(i,j) = ( ( xy_SurfSnow(i,j) - 0.0_DP ) / ( SnowThresholdForAlbedo - 0.0_DP ) )**NOrd xy_SnowFrac(i,j) = max( min( xy_SnowFrac(i,j), 1.0_DP ), 0.0_DP ) !!$ if ( xy_FlagLand(i,j) ) then !!$ xy_SnowFrac(i,j) = & !!$ & ( ( xy_SurfSnow(i,j) - 0.0_DP ) / ( SnowThresholdForAlbedo - 0.0_DP ) )**NOrd !!$ xy_SnowFrac(i,j) = max( min( xy_SnowFrac(i,j), 1.0_DP ), 0.0_DP ) !!$ else !!$ xy_SnowFrac(i,j) = 0.0_DP !!$ end if end do end do else xy_SnowFrac = 1.0_DP end if end subroutine CalcSnowFrac
Function : | |
FlagExist : | logical |
SeaIceConc : | real(DP), intent(in) |
function SeaIceAboveThreshold( SeaIceConc ) result( FlagExist ) ! 雪と海氷の定数の設定 ! Setting constants of snow and sea ice ! use constants_snowseaice, only: SeaIceThreshold real(DP), intent(in) :: SeaIceConc logical :: FlagExist if ( SeaIceConc > SeaIceThreshold ) then FlagExist = .true. else FlagExist = .false. end if end function SeaIceAboveThreshold
Subroutine : |
This procedure input/output NAMELIST#snowice_frac_nml .
subroutine SnowIceFracInit ! モジュール引用 ; USE statements ! ! 種別型パラメタ ! Kind type parameter ! use dc_types, only: STDOUT ! 標準出力の装置番号. Unit number of standard output ! ファイル入出力補助 ! File I/O support ! use dc_iounit, only: FileOpen ! NAMELIST ファイル入力に関するユーティリティ ! Utilities for NAMELIST file input ! use namelist_util, only: namelist_filename, NmlutilMsg ! 作業変数 ! Work variables ! integer:: unit_nml ! NAMELIST ファイルオープン用装置番号. ! Unit number for NAMELIST file open integer:: iostat_nml ! NAMELIST 読み込み時の IOSTAT. ! IOSTAT of NAMELIST read ! NAMELIST 変数群 ! NAMELIST group name ! namelist /snowice_frac_nml/ NOrd ! 実行文 ; Executable statement ! if ( snowice_frac_inited ) return ! デフォルト値の設定 ! Default values settings ! NOrd = 2.0_DP ! NAMELIST の読み込み ! NAMELIST is input ! if ( trim(namelist_filename) /= '' ) then call FileOpen( unit_nml, namelist_filename, mode = 'r' ) ! (in) rewind( unit_nml ) read( unit_nml, nml = snowice_frac_nml, iostat = iostat_nml ) close( unit_nml ) call NmlutilMsg( iostat_nml, module_name ) ! (in) if ( iostat_nml == 0 ) write( STDOUT, nml = snowice_frac_nml ) end if ! 印字 ; Print ! call MessageNotify( 'M', module_name, '----- Initialization Messages -----' ) call MessageNotify( 'M', module_name, 'NOrd = %f', d = (/NOrd/) ) call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) ) snowice_frac_inited = .true. end subroutine SnowIceFracInit
Constant : | |||
module_name = ‘snowice_frac‘ : | character(*), parameter
|
Variable : | |||
snowice_frac_inited = .false. : | logical, save
|
Constant : | |||
version = ’$Name: $’ // ’$Id: snowice_frac.f90,v 1.1 2015/01/29 12:12:21 yot Exp $’ : | character(*), parameter
|