module m_ECOSM_initialize_sed contains subroutine ECOSM_initialize_sed() use mod_xc use mod_necessary_ecovars use m_ECOSM_bioini implicit none include 'common_blocks.h' include 'biovar.h' include 'ECOSMparam1.h' call ECOSM_bioini if (mnproc==1) print*,'initializing bot_layer: NSED', nsed do j=1-nbdy,jdm+nbdy do i=1-nbdy,idm+nbdy bot_layer(i,j,1,ised1)=20.0*REDF(18)*REDF(1)*REDF(6) !Nitrate bot_layer(i,j,2,ised1)=20.0*REDF(18)*REDF(1)*REDF(6) !Nitrate bot_layer(i,j,1,ised2)=20.0*REDF(20)*REDF(3)*REDF(6) !Silicate bot_layer(i,j,2,ised2)=20.0*REDF(20)*REDF(3)*REDF(6) !Silicate bot_layer(i,j,1,ised3)=2.0*REDF(19)*REDF(2)*REDF(6) !Phosphate bot_layer(i,j,2,ised3)=2.0*REDF(19)*REDF(2)*REDF(6) !Phosphate #if defined (ECOCCO) bot_layer(i,j,1,ised4)=20.0*REDF(18)*REDF(1)*REDF(6) !CaCO3 bot_layer(i,j,2,ised4)=20.0*REDF(18)*REDF(1)*REDF(6) !CaCO3 #endif /* ECOCCO */ enddo enddo end subroutine ECOSM_initialize_sed end module m_ECOSM_initialize_sed